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ABSTRACT 


The forests of the Wapiti map area, Alberta (NTS 83L) were studied 
in order to ascertain relationships between forest growth, plant 
community distribution and environmental factors within the western 
boreal and subalpine forest. Quantitative data on tree productivity, 
vegetation and soils were collected from 137 plots. 

The 15 described forest vegetation types of predominantly fire 
origin were dominated by lodgepole pine (Pinus contorta Dougl. var. 
latifolia Engelm.), white spruce (Picea glauca (Moench) Voss) black spruce 
(Picea mariana (Mill.) BSP.) or aspen (Populus tremuloides Michx.), were 
floristically simple (11l vascular plants), had tree productivity rang- 
ing from 0.4 aches ae in the Picea mariana / Ledum groenlandicum / 
Rubus chamaemorus type to 5.0 ae Mee in the Populus balsamifera / 
Rosa acicularis / Thalictrum venulosum type and occurred on soils 
representing Luvisolic, Brunisolic, Gleysolic, Regosolic, Podzolic and 
Organic Orders. Two major successional trends are apparent in the up- 
land boreal forest: lodgepole pine, aspen and balsam poplar to Picea 
glauca / Rubus pubescens - Maianthemum canadense type on moist sites 
at low elevations and lodgepole pine and black spruce forests to Picea 
spp. - Abies lasiocarpa / Hylocomium splendens type on moist to well 
drained sites at higher elevations. Pinus contorta / Menziesia glabella / 
Rubus pedatus type succeeds to Ficea engelmannii - Abies lasiocarpa / 
Menziesia glabella type in closed canopy subalpine forests. 

Multiple linear regression equations were calculated for 83 lodge- 
pole pine-, 30 white spruce-, 15 black spruce-, and 17 aspen - dominated 
plots using mean annual total volume increment, stem density, site index, 


periodic annual total volume increment per tree, age of maximum periodic 
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annual total volume increment and mean annual basal area increment as 
dependent variables. Forty-nine independent variables including five 
external site variables, 21 soil variables, five stand variables, and 
18 vegetation variables accounted for large amounts of the variability 
in the dependent variables: lodgepole pine (44-66%), white spruce (70- 
94%), black spruce (98-99%) and aspen (94 to 99%). Testing the 
equations with an independent sample reinforced their validity. When 
vegetation variables were deleted from the regression, the equations 


lost much of their predictive capability. 
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I. INTRODUCTION 


In recent years a greater appreciation of forested land as a valuable 
natural resource has intensified interest in the evaluation of variation 
in forest site productivity and in inventories of a rapidly changing 
landscape. Forest site classification originated several decades ago in 
Scandinavia (Cajander 1926, Ilvessalo 1954) and other parts of Europe. 
Thorough reviews of forest site classification work have been compiled for 
the United States by Ralston (1964), Jones (1969) and Carmean (1975), 
Carmean (1977) reviewed site classification research on northern forest 
species for the United States and Canada. Daubenmire (1976) comprehensively 
reviewed the use of vegetation in assessment of the productivity of forest 
lands. Only in the past 30 years has any intensive forest-site appraisal 
research been conducted in Canada, particularly in the provinces of Ontario 
(Hills 1952, 1958; Hills and Pierpoint 1960), Quebec (Heimburger 1941, 
Linteau 1955, Lafond 1964, Lemieux 1965), the Maritimes (Loucks 1962, 
Damman 1964) and in British Columbia (Illingworth and Arlidge 1960, 

Lacate 1965). Excellent reviews of forest site classification studies in 
Canada have been prepared by Rennie (1963) and Burger (1972). 

The relatively recent interest in the Canadian forest resources is 
related to two factors. The first is Canada's youth in relation to forest 
resource development. The second is the country's immense area (9,964,000 
ne) Ofewhich over 3,222,700 na is forested land (Forest Econ. Res. Inst. 
1972). This latter factor has in the past lead to belief that there will 
always be enough forested land to provide Canada with lumber, pulp, paper 
and other forest products and benefits. Now there is realization that our 


forested land must be carefully managed if it is to remain productive. 
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Forest-site research in the western boreal forest (Alberta, and 
Mixedwood Section (Bl8a; Rowe 1972) of Saskatchewan and Manitoba has fallen 
behind other parts of Canada and the United States due primarily to the fact 
that only very recently (the last 20 years in Alberta) has any real utiliza- 
tion demand been put upon the western Boreal Forest Region. The growing 
concern that our forest resources are too rapidly being depleted is promp- 
ting here, as elsewhere, the search for means to most effectively utilize 
our forest resource while still maintaining its productivity. 

Some of the earliest work in the western boreal forest was done by 
Brinkman (1931, 1936) in Alberta. Smithers (1956) assessed site productivity 
in dense lodgepole pine stands in the Kananaskis Forest Experiment Station, 
Alberta. Duffy (1964a) used multiple regression techniques to find rela- 
tionships between site factors and growth of lodgepole pine in the Alberta 
foothills. Duffy (1965) developed a forest land classification for the Mixed 
wood section of central Alberta on the basis of differences in soil parent 
material and soil moisture status as they influence white spruce site 
index. Jameson (1963) compared tree growth on two sites in the Riding 
Mountain forest experimental area, Manitoba and in 1965 he related jack pine 
growth to site in the Mixedwood forest Section (B. 18a) of Saskatchewan. 
Dumanski et al. (1973) evaluated the productivity of lodgepole pine forests 
using soil survey maps for the Hinton - Edson area. Lesko and Lindsay (1973) 
related lodgepole pine and white spruce site index within 15 forest types 
to soils in the Chip Lake map area in west-central Alberta. Kabzems et al. 
(1976) described 23 forest ecosystems in the Saskatchewan Mixedwood Section 
(B. 18a) with respect to vegetation, forest productivity and soils. In 
addition, several descriptions of vegetation distribution in northern and 


north-western Alberta have been made by Lewis et al. (1928), Dowding (1929), 
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Raup (1933, 1934, 1946), Moss (1953, 1955), Moss and Pegg (1963), LaRoi 
(1967) and most recently by Achuff and LaRoi (1977). 

In recent years attempts have been made in Biophysical Land Class- 
ification (Jurdant et al. 1975) to differentiate and classify ecologically 
significant segments of the land surface, rapidly, on a small scale 
(reconnaissance survey), to satisfy the need for an initial overview and 
inventory of land and associated wildland resources (Lacate 1969). Bio- 
physical inventories have been conducted in Quebec for several years (e.g. 
Jurdant 1969), and more recently in National Parks (Coen and Holland 1976, 
Wells et al. 1977). 

Biophysical inventories or independent surveys of soils and vegeta- 
tion in conjunction with aerial photography have been used as the basis 
for the production of small scale interpretive maps relative to forestry, 
agriculture, ungulates, waterfowl, sportfish and recreation in the Canada 


Land Inventory program (McCormack 1968). 


Site Evaluation Methods 


Site evaluation studies may take several approaches (Rennie 1963, 
Ralston 1964, Jones 1969, Brickell 1975). A determination of site quality 
may be made by direct methods (measurement) or by indirect methods 


(inference). 


1. Direct methods 


The direct methods of site evaluation usually involve establishing 


plots within mature forest stands in which plot maps may be prepared 


showing the species and position of every tree (Rennie 1963). Various 


‘ Aj 23 , SELL ,CERL) quiet bi 


; 


tnpoe4 seen Cie | 


ay sy ‘sagoes OT 


(oe 


a i) 


itehewt) nod fat z. a 


rial ¢ >. ginemoen dagod tee 
ee “ - 
a ev us eos e LSAnooias 


bate 720En24 Of bis 2 yrasn 
~ 

edairocneveas Ieoreg 
&yive Dts , (bees trie bye 
ror .f9 #0 site 
' 

* ‘ 

oivesress i deobeyieot a 


Nol - 
+6 w mvitoautaoo al RRR 


4 4 - U ' a 


at ey fese Liema 30 BYesa0 O19 oft 302 2, | 
q 


ut 


i 


Iwotaestéw (aedélevns ,o1vd foe 


wi 


ec 


6, welt wiseuleve eff 
30 ¥ oe! sf0L . bee nosealag 
eeen) ebodyen Joe1ih vd ebem oc Yar 


one xzesInt) 


." 


ab or ~~ bgt 
abovidom vetid «4 
pmo ? , z sunib ett 


mensurational data such as height, diameter and age for all trees of a 
stand or a representative sample of the population are taken. More precise 
estimates of increments in growth parameters might involve remeasurement 
of the same trees at five-or ten-year intervals Productivity indices 
obtainable from these types of measurements may be expressed as total 
volume per unit area, basal area increment, annual or periodic volume and 
height increment or site index. Site index or total height at a specified 
convenient age has been the most widely used growth index in site evalua- 
ELON. 

In recent years refinements in direct forest productivity assessment 
have been made through the use of aerial photography. Continuing new 
development (Gimbarzevsky 1975) in false color infra-red remote sensing 
imagery and the assignment of spectral "Signatures" that characterize the 
light reflectance of individual tree species will no doubt prove even 


more useful in this regard. 


2. Indirect Methods 


Indirect methods for estimating site productivity utilize a related 
attribute as a criterion. Five methods may be recognized: climate, 
ground vegetation, soil properties, foliar characteristics and site mapping 
(Rennie 1963). Site mapping tends to be an integration of the former 


approaches. 


Classifications based substantially upon climatic criteria are by 
necessity of a very small scale, eg.Rowe's (1972) Forest Regions of Canada. 
Krajina's (1965) Biogeoclimatic Zonation of British Columbia is a similar 


approach though it gives more detailed climate, vegetation and soils 
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descriptions of the zones than does Rowe's (1972) treatment. 

Studies relating forest productivity to vegetation type have had 
varying degrees of success. Cajander's (1926) work in Finland is a classic 
example of demonstrated productivity-vegetation relationships. Daubenmire 
and his students (1942 et seq.) have successfully related forest productivity 
to habitat type, Canadian examples of studies where vegetation type- 
productivity relationships have been shown include those of Linteau (1955) 
in Quebec, Rowe (1956) for the prairie provinces, Illingworth and Arlidge 
(1960) in British Columbia and Lesko and Lindsay (1973) in Alberta. Duffy 
(1964b) has comprehensively reviewed concepts in plant sociology with respect 
to their use in forest site classification. 

Ordination approaches to site or to vegetation classification treat 
the vegetation as a continuum which has developed in response to environ- 
mental gradients of climate, soil moisture, nutrients and other site 
attributes. An ordination technique often helps elucidate the factors 
responsible for the differences between stand species composition and 
productivity. Examples of ordination approaches include Maycock and Curtis 
(1960) and McIntosh (1967). 

Site evaluation studies relating soil properties to productivity take 
one of two approaches (Ralston 1964): interpretation of forest productivity 
from conventional soil survey data, and derivation of field guides based on 
empirical correlations of site attributes and tree growth. The latter is 
analogous to what Jones (1969) refers to as the factorial approach. Forest 
soils display a whole gamut of characteristics which in general are diffi- 
cult to measure and call for experimental skills in a number of disciplines 
(Rennie 1963). Numerous soil characteristics have been correlated with 


forest growth as noted in the literature. Specific examples are 
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dealt with in the "discussion" section. More conventional criteria studied 
include pH, available nutrients, moisture regime, texture and soil depth. 
The quantitative effect of soil factors and their interaction in relation 
to tree growth (usually expressed as site index) was first tested statisti- 
cally by Coile (1935) utilizing multiple regression analysis. Most studies 
of this type since that time have dealt with physical soil properties, 
though virtually any factor can be included as an independent variable. 

Another indirect method of evaluating one aspect of site quality is 
foliar analysis. It is more useful on nursery stock or young trees, and 
tends to be a reliable criterion when particular nutrients become acutely 
deficient, rather than only slightly so (Rennie 1963). 

Forest site mapping takes several approaches. The map units may be 
separated on the basis of a single attribute such as soil texture or more 
commonly are a synthesis of criteria which reflect differences in the main 
ecological features of the habitat,ie. regional climate, topography, geology, 
soil properties and vegetation cover. The ways in which forest site mapping 
methods differ are in the interpretation of the environmental factors, the 
actual quantitative measurements made to specify them and in the emphasis 
given to different factors (Rennie 1963). The usefulness of a conventional 
soil survey for site evaluation depends upon whether the taxonomic units 
defined are sufficiently uniform ecologically, which in turn depends on the 
degree of coincidence of the classification criteria with the characteristics 
of greatest relevance to problems of forest management (Rennie 1963). Hills 
(1952 et seq.) has used a physiographic approach to site, with landform as 
the basic map unit, with subdivisions according to climate, moisture and 
nutrients. 


The most useful forest productivity maps are likely to be a result of 
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the consideration of as many ecologically significant factors as feasible. 
The present study was initiated in 1972 in conjunction with a 
reconnaissance soil survey of the Wapiti map sheet area (NTS 83L) and had 
the following objectives: (1) To determine the effect of various vegetation 
and site factors upon several measures of forest productivity for the major 
tree species within the study area; (2) To compare the forest productivity 
of the soil mapping units; (3) To describe and classify the major forest 
vegetation types of the Wapiti map area in relation to the soils; (4) To 
describe and aeeenan He the successional relationships of the major forest 
vegetation types of the Wapiti map area. 
It is hypothesized that: (1) the various vegetation and site 
factors are correlated with the forest productivity measures in such a 
way that they can be used in regression analyses for the purposes of 
tree growth prediction; (2) the productivity of the soil mapping units 
differs significantly among the units. In addition it is hypothesized 
that (3) vegetation types can be classified; (4) major successicnal trends 


are discernable. 


-eldiees? 24 esasae? Saesis cre 


& Siiw Ae, aoe tie 


ynan’ de 20 Soe 


etw yous 


S yitanipelots 


+ 
“ 


iz ot ai Bepeiaie 


a 4 ; ’ Aig 
batt BAS (JES UTA) Ree Pesce FaH tare Sts Be yevade Die 
2 i ae LO 
nolessedss Buoktsy. jo Japtze 5 ankn? Bb OF (2) ante pe. 
wetem, ott sc} VIAN LIO8 & jn59c3 to asiciases “1S650Ns6 oe nal 
edin2aphT7q: } S103, Br? 24561 aT > ghee YOots. sa ALAdhw Bae 
Saeco moran os 9 DOA a > oy 3ORS5 oF tt) ee BA nation ug 
> , algae ‘ weg 
oT (h) irien ges or necsselad fs aoe ceo cy wie? aod ie eegy? 
$ae 152 .210Tam ot i eoidanoiieie? Langer s2cd2 ens secon a 
6 gat iaFgen’ Sts ae sagyy 7 
an 
eve ‘brits notistacss moiser et (4) St Se 
& Hogeenh eeavisem vorvitapbesq Seesei Sas avi neh toni § 
Zo. asecttug elt xcs cluak sored ss. 6 Beis" ‘ect caheta cai 
ay 7 
atin pakecen Loox ait Ao wiwitedborg eng (i _ dedste beara Toes 
rn ty ae 
beavers tiie al 24 ne tihes, susie ent aman, 


apne ‘es ia fate 


Avie. FES a oe 4a, 
BaPooos AUS TP. 


II THE WAPITI MAP AREA 


A. LOCATION 


The Wapiti map area (NTS 83L) is located in western Alberta between 


118° and 120° w longitude and between 54° and 55° N latitude covering an 


area of approximately 17,460 aie (5740 me jE LCs Lis at ls: bordered ito 
the west by the British Columbia boundary. The town of Grande Cache lies 
just below the south-western quarter and the city of Grande Prairie is just 


north of the area. The forestry trunk road closely parallels the eastern 


boundary. 


B. SURFICIAL GEOLOGY 


Both Cordilleran and Keewatin ice sheets covered parts of 83L. The 
influence of the Cordilleran ice sheet is restricted to the hills and 
valleys in an area that includes between one-quarter and one-third of the 
map area in the south and west (Bayrock 1972). Virtually all of the 
remainder of the area was covered by the Keewatin (Continental) ice sheet. 
At least two glacial advances occurred. 

Meltwaters flowing by the ice margin deposited thick beds of gravel 
of Pleistocene age. These gravels were subsequently resorted and 
redeposited in river valleys , especially those of the Smoky, Kakwa, Cutbank 
and Simonette rivers. These glaciofluvial deposits are in the form of 
river terraces of coarse gravel, outwash sand and gravel, the proportions 
of each varying locally, and ice-contact deposits such as kames and eskers. 
The latter account for a very small proportion of the Wapiti area. 


A large proglacial lake of Pleistocene age formed in the north and 
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east of the area and gave rise to the large amount of glaciolacustrine 

delta sands, bedded silt and clay, and bedded sand and silt. Aeolian sand, 
derived from glaciolacustrine sand occurs near the glaciolacustrine deposits 
in sheet and dune form. 

A thin deposit of calcareous till derived from the Cordilleran glaciers 
occurs in a fairly large area of ground moraine on both sides of the 
Narraway River; very stony (greater than 30%) cirque and valley glacier 
deposits in the higher foothills; stony till deposits within dissected 
lateral moraine of large valley glaciers, especially in the lower valley of 
the Smoky River; and small pockets of thin (less than 50 cm to bedrock), 
stony (15-30%) pre-Wisconsin Cordilleran till, in which carbonates are 
leached to a depth of 2m, on the level elevated benches. 

Till derived from the Continental glacier includes a few small areas 
in the north and east of the area of hummocky moraine of clayey to sandy 
texture with numerous erratics from the Canadian Shield. Most of the 
Continental till which covers about one-third of 83L is ground moraine 
of clayey to sandy texture, contains numerous erratics of Canadian Shield 
Origin, and is of gentle relief. A fairly large area of thin, very stony 
ground moraine of reworked Tertiary gravel on level topography is present 
on the high plateaus in the southeast corner of the Wapiti area. Large 
areas of this ground moraine of continental origin but lacking erratics 
from the Canadian Shield, are present to the west of the majority of the 
Continental till. This thin till is made almost entirely of local bedrock. 

A few small areas of shale, sandstone, siltstone, coal and conglom- 
erate outcrops are present in the high relief areas in the southwest 
corner of the map area. Shale and sandstone are the predominant bedrock 
types and in general it is these materials that give the topographic 


relief to the landscape to a much greater extent than the surficial deposits. 
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C. BEDROCK GEOLOGY 


Information on bedrock geology in the Wapiti map sheet area is scarce. 
A report by Allan and Carr (1946) describes geology and coal occurrences of 
the Wapiti-Cutbank area. The following general account (1: 1,584,000 map) 
for Alberta is given by Jackson (1975). The Wapiti map area lies within 
two major physiographic regions; the Interior Plains Samide Western 
Cordillera. The Rocky Mountain Foothills section of the Western Cordillera 
is represented in the southwest corner. It is characterized by shales, 
conglomerates, coal, chert and sandstone of Triassic to Lower Cretaceous 
age and Upper Cretaceous shales, concretionary shales and sandstones. 

The Interior plains region is east of the Rocky Mountain Foothills 
and is represented to the south and west of the Wapiti map area by the 
Western Alberta Plains section. This section is characterized by Upper 
Cretaceous shale and Tertiary shales, siltstones and sandstones. The 
Interior Plains Region to the north of the Wapiti map area is represented 
by the Wapiti Plain which is characterized by light grey, Continental and 


deltaic siltstones and shales of Upper Cretaceous age. 
Dee SOTLS 


Representatives of the Luvisolic, Brunisolic, Gleysolic, Regosclic, 
Podzolic and Organic Orders occur within the Wapiti map area. The soils 
were separated on the basis of differences in parent material into 28 
"soil groups" (Twardy 1978). These "soil groups" that form the basis 
for the soils map, were further separated into 66 "soil units” 

distinguished on the.basis of internal drainage, soil subgroup classifica- 


tion and the proportion of other non-mappable inclusions within 
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the units. The ten most important soil groups with their associated parent 
materials in parentheses, listed in order of decreasing area are: (1) 
Edson (Continental till less than 975 m ASL) (2) Donnelly (lacustro- 
till), (3) Mayberne (Continental till greater than 975 m ASL) (4) Copton 
(coarse residual), (5) Torrens (fine residual), (6) Marlboro (medium 
textured Cordilleran till), (7) Robb (moderately coarse to medium textured 
thin Cordilleran till and colluvium), (8) Lodge (thin Edis aelnid ad) 
deposits or aeolian over lacustrine or till), (9) Blackmud (sandy 
alluvium), and (10) Heart (aeolian). 

More detailed descriptions of the soils of the Wapiti map area and 


their extent are discussed by Twardy {1978). 


E. STREAM CHARACTERISTICS 


The main streams of the Wapiti map area are the Smoky, Wapiti, and 
Simonette Rivers that drain toward the north. The Wapiti and the Simonette 
join the Smoky east of Grande Prairie and the Smoky joins the Peace River 
near Peace River Townsite. Important tributaries of the Smoky are the 
Kakwa and Cutbank Rivers which lie almost wholly within the study area. 

The Narraway River is a major tributary of the Wapiti and has its origin 
just to the west of the study area. The Latornell River, a major tributary 
of the Simonette lies wholly within the area. 

With the exception of the larger rivers, most of the tributaries 
flow in steep-sided, narrow valleys and are still down-cutting, indicated 
by the absence of floodplains. The Smoky and Wapiti Rivers have developed 
flood plains at their northerly extent. Several other tributaries in the 
region have developed narrow floodplains but very little deposition has 


occurred elsewhere. Surface water and groundwater is abundant in the 
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Wapiti area as shown by the innumerable perennial streams. The stream 
pattern is dendritic. Floods are common during the mountain snowpack melt 
and flood-related damage can be extensive in terms of washed-out roads, 
bridges and railroads, and severe soil erosion on exposed slopes. The 
heavy rainfall of June 1972, coinciding with peak snowmelt from the Rocky 
Mountains provided many examples of such damage. The higher elevations of 
the Wapiti map area to the south and west yield most of the water with less 
fluctuation than the lower areas as indicated by the greater stream 
density in the south and west. 

Sediment loads in the streams during early summer are large and 
indicate excessive erosional losses on the watersheds. Unfortunately, 
few quantitative data are available, but most of the larger rivers remain 


clouded by suspended sediment until at least mid-August. 


e.) CLIMATE 


The climate of the Wapiti map area is characterized by cool to warm 
summers and cold winters. The area lies within three summer climatic 
Bones (Maciver et al. 19727 Fig. 2). ‘The:northern third of 83h (Fig. 2) 
is characterized for the period including the growing season between May 
1 and September 30 by a mean air temperature of ierah ey an average 
precipitation of 30-41 cm, a potential evapotranspiration of 43-46 cm and 
a frost-free period (greater than 0° C) of 100-120 days (zone 4, MacIver 
et al. 1972). A belt running from west to east through the Wapiti sheet 
and representing approximately one-fourth of the area is characterized for 
the May 1 to September 30 period by a mean air temperature of LO=11ec; 
precipitation of 36-46 cm, a potential evapotranspiration of 41-43 cm 


(ie. a net water surplus), and a frost-free period of 100 to 120 days 


' | ~ 
ree eee pibincmniens 902 
+ bem eeepHine SRE WN pil ruth seme eve annals, | 
eaeas guceniasna Ye eAes? ts ayvienertal a reo 


ai? <epgala late cree® at aoaete 2460, faved. fem 1 
cpap Rete, Guba Beet eh le ann Ge £ 


yong gay guise, ITenor 


des cescoedbeeiemeld at ..eqceet tame: 2o Me gRE Tt FY 


enult ade (eeee oe deer | (TEP E show “Sew: Mee ead oF Baie 


gabe Hodaer} 202 *- t-» fASeeas- Oe 10 ada hs 
Fale TW cnn 


sghiwe “egests oi! a/ ahgok 


Gi: Pca? aah Peewee (226? 


pibnRaaer he ARR Sdap ata Gee A yeeoot Jeac regan 


l ght Zo Teor gud eeltal cerns eis bk BV 


Tr Lay 4 maul a Do of Rha srSwi toy Nets pega 


PER. mn), 


niéeea- SRE. FAX 


q 
’ 


cate os foo yo bes nse oe tan a} nese yen Eee ae oe 7 ; 

; y vy} | | / ke a 
Pale Ween ciahiers eI Hl At pe hat gate ott" cateeedy uo i 
: iy ot 1 

& 


ne 


4 YY 
: mt m io yak Oe a he ee ag teyioem) 
> ie, phan Wea Nn bwitt maentany 8c ie ieee Ns : re) 
; a ed 


Vai 4 Feiilpeiseons noa6 36 CA waz Ls) 5a3 If thy aT ee? | ye4q ht 42: 7. =) cr ee 
| - 5 SS Oe) aoe 
w ; maak) 
oy - ~ ee HPeAove a. UL-GL to Savisieqees 48 1698 4. yd OF. “Wh 
Pe fe | ba i : 
D 


, eas 
ade we cs ch Sg ighisesigegeticqar’ [si ieetoy = ym THO $o-7 aide: a 


i 7 


bi a ; 
Be began) isl Siete. te 19 “o reds goaderyl gaits itt 
uae. nr sS 
: ‘ : J 7 +7 ae : 
- A af Sagceias joes 03 3ug0 omit pataaws Sfed A ‘ oa 
4 a Basta) ah gots 97 to Stud ove ¢istahkivorqye Lessee ave hee 


| Lanerreayes She niet & bonis OE rmeases 1 ye 
oy 
ae a mlaacaess <7apP4 besseay 6, yr ver ak ses om 762. nJigh 
, a 
o hahaan siasiacens a. ail  entaawn a8 eae 
7 
_ 


4.) Lia (ie ae 


~ 


r4e 


ro Maciver af Gl WFD 


é “d a 
4 i) Ks ’ a4 
‘ 
14 { e 1¢—@4 
Fo beri bes Si‘ec4 


Marnieanepicetis  «: 457,46 at 
cm or a 
oe ea aes : a ba ar ay * nt ET Aisi P 


x 


BAOOTHED STATISTICA, CLIMATIC 
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Figure 2: Climatological zonation of the 
Wapiti map area 
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SMOOTHED STATISTICAL CLIMATIC ZONES (redrawn 
from Maclver et al. 1972) 


Selected Seasonal Normals (1 May - 30 September) 


1954 - 1968 
10) ae ; Ox 
Group Temperature ( C) Precipitation Degree - Days ( C 
Number Mean Maximum Minimum cm >-2 0 6 
2 - 10-11 14-17 4-6 ~ 41-46 32-36 24-28 10-14 
4 Ll—12 17-19 4-7 30-41 30-40 32-34 16-20 
10 10-11 16-17 6-7 36-46 30-38 26-32) 12-13 
Potential Frost-Free Kali, ESOS 
Evapotranspiration Day Length Period Free Period 
cm Hours Days Days (7-2°C) _ 
2B 38-43 Wisi6 diese 80-100 100-140 
4 43-46 BES }easiciesi 100-120 140-160 
70 41-43 L5.2=15..4 100-120 140-160 
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(zone 10, MacIver et al. 1972). Approximately the bottom one-third of the 
Wapiti map area has for the same months, a mean air temperature of L0s1i-c, 
precipitation of 41-46 cm, a potential evapotranspiration of 38-43 cm and 
a frost-free period of 80-100 days (zone 2, MacIver et al. 1972). It 

thus becomes apparent that as elevation increases from the NE to the Sw, 
precipitation increases and potential evapotranspiration and air 
temperature decrease (Fig. 2). We thus might expect the most favorable 
climates for tree growth to be in the northeast where temperatures are 
warmer and precipitation is adequate. 

The nearest, full-time weather station is at Grande Prairie to the 
north Gee ier wat. 118° 53' long.). Year-round temperature records are 
made at only one location in the Wapiti map area, The South Wapiti ranger 
station in the north of 83L. The year-round temperature records at the 
Muskeg ranger station below the southeast corner of the area are probably 
representative of much of the southern portion. Weather records from May 
to September inclusive are available for eight forestry fire towers within 
the Wapiti area. The extrapolations made by MacIver et al. (1972) were 
from these data. 

Winters in the Wapiti map area are cold, with January (the coldest 
month) mean air temperatures (mean of daily maximum and minimum) Apa gh ee. 
C at Grande Prairie airport, LG a0n Car Muske J and -15.4° C at South 
Wapiti. Mean monthly air temperatures are below freezing between November 
and March inclusive (Environment Canada 1975). Winter precipitation 
data are available only from the Grande Prairie station. Total precipita- 
tion for the November to March period is 18.5 cm (mostly snow) of which 
3.4 cm fell in January and 3.1 cm each for November and December (Environ- 
ment Canada 1975). The November to March precipitation total represents 


42% of the total annual precipitation at Grande Prairie. In the 
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predominantly higher elevations of the study area south of Grande Prairie, 
winter temperatures would probably be slightly lower and precipitation 


slightly higher, depending upon location. 
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Lit) METHODS 
A. FIELD SAMPLING 


1. Plot Selection 

From an examination of Alberta Forest Service (1: 126,720) forest 
cover maps and aerial photographs supplied by government (Alberta Forest 
Service) and industry (Proctor and Gamble Cellulose ned plot areas 
were chosen to encompass a wide variety of vegetation, soil and landform 
types, within uniform, even-aged and normally-stocked stands ranging from 
45 to over 200 year old. Sampling was concentrated within 70-100 year 
old stands of the various forest types. Most plots were readily accessible 
by road, though four-wheel drive and the Cramer "swamp buggy" were used to 
get to less accessible areas via seismic lines and rough trails. Sampling 
was done primarily within the Upper Foothills (B. 19c) and Lower Foothills 
(B. 19a) Sections of the Boreal Forest Region (Rowe 1972) and to a lesser 
extent in the East Slope Rockies Section (SA. 1) of the Subalpine Forest 
Region. Treeless and alpine areas of 83L were sampled with much less 
intensity than the remainder of the area, and ground work was largeiy of 


a descriptive nature. Sampled plots totaled 137. 


2. Plot Sampling 
Zak SOLS 

In each plot a soil pit was dug in a position with a soil profile 
representative of the plot area. The soil profile was examined and describ- 
ed noting structure, texture, color, thickness of horizons, pH (Lamotte- 
Morgan field kit) as well as plant rooting characteristics and internal 


drainage according to CSSC (1970). Prevalent slope gradient and aspect 


were also recorded. 


19 


sqevor (fSTyesk t+) aai tre seen Lestat rA\> 


Juez 7i (62 2801") Momiemnan aot Pea hager ne wince 
(oa 


seu7s SIG, (+074 sq fubta silanes, ors vas peae) sini “7 
a 


is JE¥ sno Fo yah idee obit 8 eaeqmocnnl Of ¢ 


gicshnsl Hab toe. ye. 
‘ By) a 


s—¥) Leascin S08): Se 6-598 | eaetroM Pit: . 


gost. eriprss ehusze hersgs4- 245 
skew OOTHOV fLARLW DS IR s3 Nene 2sW Snbiemae* Bie Tey, ue vo 


Leni Jséih’...egees san7e abodsey adie to tn 


shabeadous ylibaos Sse" eo928 ~™ * 
a a: 


os Basi 976" “VE ppG che Tom ane 22: Sug Sag ray: bab tdone~so® € 
| ; , Z BAe 
part Apwod sts spats OLPSHSS PATE aos ie 


oats «abet >5 wie 
fi sy a iN 
4 yoaety sao ita iw ee 


3). abbeasea 


aitsnavoo? azwed onF (se: 


wabeelie te She (S° Sh. wor), om COA "Faaio7 lseton, wth p>, Ais 
Scan whachoona Sts 30) 12. Aeh e ree no ites ort wee e 
Rae: foun ihe Soi gie2 eM PG te ee aa fonhats: Bares webs ee 17 
ta yaserse! car Tor anode” br daae $43 390° Sebipignes ond, aie 7 
ha y TiS bg IEIeF athe bedigma Laaonaes cial 
a Be 
oe ot iy, 
a 
izhoxe fies 2 WIiw ait leo? . fot eu eeW 224 Lica oe i“ ae 
“adinvanate ane Beintmers Pow ‘st ridag Lise sit «Esis #649 Basing ame jose saiad 
: ~* . 
g 7 ~eegonaly AY lanostvsn RG qrsnags As saotes co ati a y 
‘ 2 = a 
a Jeoasiri bie saree nit bre Poo + tr pats Fo sane 
a. 1 a ieaal Sasunar’ “aIsit ina tevert orety, ‘Saag itt 
: as | - i m3 nie 
; wt wi - o 7 


2.2 Vegetation 

Within the circular (0.405 ha) plots, the diameters and heights of 
all trees over 1.3 cm diameter at breast height (1.5m) were tallied by 
species. When the tree tally was completed, five to seven healthy 
dominant and co-dominant trees were selected and felled, measured for height, 
and sections removed at 0.3m, 1.5m, and 1.8 or 3.7m intervals thereafter for 
stem analysis. In uneven-aged stands (dominant trees differed by more than 
20 years) additional discs or increment cores were taken from the largest 
trees in the plot for determination cf maximum stand age and evidence of 
stand history (fire, climate, disease, anomalies in growth patterns, etc.). 
Tree canopy cover was visually estimated at 12 random points in the center 
of the quadrats used for the subsidiary vegetation. 

Subsidiary vegetation includes plants other than trees greater than 1.3 
cm DBH, and was sampled by height strata. Estimates of per cent cover by 
species, using a modified Braun-Blanquet scale (Appendix 1), were mace in 
12 randomly placed 1lxlm quadrats for terrestrial bryophytes, herbs and 
dwarf shrub species (up to 0.5 m tall). Shrub cover (over 0.5 m tall) 
estimates were made by species within 5x5m quadrats centered around the 
lxlm quadrats. Tree regeneration density was tallied by species and height 
class within the 5x5m quadrats for individuals less than 1.5m tall. 

A species list for the plot area as a whole was compiled. Unidenti- 
fied specimens were collected for determination in the laboratory. The 
predominant height, vitality and phenology were noted for each species on 
the list. Additional notes were made of predominant plot aspect, slope 
relief, deadfall cover, evidence of disease, animal activity, fire history 
and of any features that would make the area particularly attractivextor 


recreational or natural area uses. 
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B. DATA ANALYSIS 2i 
1. Stem Analysis 

Stem growth data were computed using the stem analysis program of 
Pluth and Cameron (1971). Growth measurements were made at 10-year 
increments (20 years for very suppressed trees) from the outside of the 
stem toward the center on the maximum diameter and at 90° to the maximun. 

The program took the mean diameter as a basis for computation. A tree age 
adjustment was made to account for a one-foot stump height, using 10 years 
for the three Picea species, five for lodgepole pine and one year or no 
correction for aspen. Stem section volumes were calculated as right 
circular cylinders except for the uppermost section or terminus which was 
always calculated as a right circular cone. Computations for the height- 
diameter and height-age curves are derived from procedures outlined for 
stem analysis by Avery (1967). The program calculates periodic annual 
increments for the ten-or 20-year time periods from the most recent period 
backwards in time. Basal areas are calculated as the cross-sectional areas 
at 1.5m from diameters measured inside bark. Using a CalComp plotter, the 
program plotted curves of height-diameter, height-age, volume incrementage 
and basal area increment versus age for each tree (Fig. 3). 

Means and standard error of the mean were calculated for the following 
statistics taken from the individual tree stem analysis data for each species 
sampled in each plot: (1) site index (at 70 years) (2) total volume per tree 
(3) periodic volume increment (70 to 80 years) (4) maximum periodic volume 
increment and (5) age of maximum periodic volume increment. For sample 
trees with ages less than 80 years, linear extrapolations were made for site 
index and periodic volume increment from the previous height and volume 


growth increments. 


2. Plot tally 


A simple computer program facilitated computation of forest stand 


statistics 
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Figure 3: 


Calcomp plots of stem analysis data for one 
tree showing mean annual basal area increment 
vs. age, annual height increment vs. age, 
mean annual total volume increment vs. age, 
height vs. age and height vs. diameter 
relationships. 
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from the plot tally data. Gross volume of trees greater than 1.3 cm DBH 
was computed for each species separately and then summed. Merchantable 
volume of trees with a minimum 7.6 cm diameter top was computed for each 
species separately and then summed. Gross mean annual total volume increment 
(MAI) was determined by dividing the total gross volume by the average 

age of the dominant trees sampled. This tends to give conservative MAI 
values in uneven-aged stands because the volumes of the smaller, younger 
trees are divided by the same age as the larger and older dominants and 
co-dominants. The coefficients used for the volume computations are those 
determined by the Canadian Forestry Service and used by Alberta Forest 
Service. Basal area at 1.5 m was computed by species then summed. Tree 
density (stems/ha) was determined. The program tallied the density of 
dominant, co-dcminant, overtopped and suppressed trees when this information 


was recorded with the tree height and diameter. 


3. Vegetation-Environment Relationships 


A two-dimensional ordination was constructed using the method of Beals 
(1960) and an index of vegetational similarity (Bray and Curtis 1957): 
I=2W/atb, where a= sum of quantitative values of all species in stand A, 
b= the same for stand B, and W= sum of quantitative values the two stands 
have in common. Prominence value (Beals 1960) data for all vascular species 
were used in calculating disimilarity indices (I-l), where PV= mean cover % 
X quadrat frequency %. 

A vegetation synthesis table was constructed with the plot numbers 
arranged horizontally and the species listed vertically in order of 
decreasing constancy (no. of stands in which species occurred) from top 


to bottom. Plots were arrayed according to floristic similarity using the 
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Braun-Blanquet method (Mueller-Dombois and Ellenberg 1974 p. 177). Similarity 
of soils and other environmental factors were also considered in making the 


stand groupings (vegetation types). 


4. Stepwise Multiple Regression 


A stepwise multiple linear regression of the abbreviated Doolittle 
method (Steele and Torrie 1960 p. 289) was computed to relate several 
measures of forest productivity to soil, site and vegetation data taken 
from the stem analysis plots. By definition, a regression equation supplies 
estimates of population means. In practice it may also be used to predict 
events (Steele and Torrie 1960). The multiple regression equation is 
expressed in the form Y = atb, x, +box5+...b x where Y is the predicted 


value of the dependent variable; a is the intercept of the regression 


line on the vertical axis; b b eee ed, are the partial regression coeffi- 


le ipakate 


erents OL VY, On X37" X_, X 


l 2! cee eK are the observed values of the independent 


variables. The program enabled one to select the dependent and forced 
independent variables and to specify which variables are to be omitted 

from the regression. Prior to running the regression the various independent 
variables were plotted against the dependent variables using a computerized 
bivariate plotting program. Where nonlinear relationships between one of 
the independent variables and a dependent variable were apparent, the 
appropriate transformation, linear or nonlinear was applied to the in- 
dependent variable to best approximate a linear relationship between the 
dependent and independent variables after the transformation. The sequence 
addition of the independent variables in the regression is dependent upon 
how much the independent variable contributes to the R? value; the greater 


the contribution, the greater the importance that variable has in accounting 


“_ Ng > Ae a 


YN ESetense  AKTTP 2g sxe hone co ® 


pea 


ody, pAbaws AD ‘Deyebigned hele wee srhitrbs a lee - 


wii) Like? Gate l veages eed. & AOL bas RSE : 


. } 
aed Shyh MOL ves ee ols) CHV eaeeas ae. 


262 gi !'V2 in SoMa iano RIT!) oS yh the EYE WE oli 


ta24 ovSeeu or! oa Be : e355 4 tite si 


5) rie. Leheps iors g3"pa4 eure Lr eas a0 oy eee ‘pth 


py pu PSs '4r 2° S200 wat i 1. ae anh ha Se a 
- a se cy 
oie ‘ hi Sta ae ee , 
WILAS Is \ “epson ge. att, OF E cod Thala , y 
ae > : Ay | 
\ ob 
a4 
as a9 “a os Pea ; ~ y ‘ 3 pa : 
47 , I 
f mney Wish, | Sina ) 259 £0) Pest iene pert D Osta. Ra vad 


Tenes4)5) Baa (ebasgeel an 4 5 3F “4 Sic ieee babuie 
Geidiselo!) Ot uve wejcoliey fataw pia seae NY (fale EL 
‘ eobupetra: ENOM 6 (Gi?  ciWis Cat oad. SAS ce cami er or 
Cee Ses iics.> Ne Poles. > OV geese as per 
~ 2 ea S707 e) BSE NeisITs { me 1st) S3er yy 


‘ a a4) abs CASTES CW ioe TE se ‘eRaTae " o% 


mil. Bay 4O7 is | ; we = ) eual tor 3G 16 24i7 


AT Snaweagd ruby oP: <4 eee & ey Peace Sectis: 7 ia 7 . 
i 3 7 oo - —4 - air ioe 
yo 
TSN Meany. Or Ien sla aeys wh Yet3y vane Py ne - i bn geebe 


Wy : 
; } \ : as, 
i = 3 peli Ste 9 4 ity ORnasudc 7 = id mJ ‘qurua Glas ya old & 
: - = | 
: ‘ 
TERE Te ae: TON a yy Hap! see (al frites Se 
i 7 : : 


_ , t 


- UI mings fi Bht.afs ant ta Seats sone oma 4 


Ss ae es 


26 


for the variation within the data supplied for the dependent variable Re 
(the proportion of the variation in the dependent variable that is accounted 
for by the regression). The regression equations presented for lodgepole 
pine-, white spruce-, black spruce-, and aspen-dominated (according to 
basal area), plots give the dependent variable as a function of the nine 
independent variables accounting for the greatest proportion of the R? 
value. 


The level of significance of the regression was determined from the F 


ratio calculated by the program. 


5. Duncan's Multiple Range Tests 


Initially, soil associations, vegetation type productivity (MAI and SI), 
elevation (nearest 150 m), slope class,soil subgroup and slope aspect were 
cross-stratified against mean annual increment (rotation age) and density. 
Because of the relatively small sample (137 plots), and because some tree 
species are not found on certain soil associations, it was not possible to 
evaluate productivity for all classes listed above using Duncan's New 
Multiple Range test (Duncan 1955). The growth prediction equations thus 
serve as a better means of evaluation in these cases. The stratification 
using Duncan's tests to compare mean annual increment by vegetation type 


proved more satisfactory. 


cigeaney a a St ey sen mcr avibe ot 


trertuawien wR nite ated imales 


agceeankeok yo) Veswoke74 endo line ries goth 
Pud An alr ae Jah CAR! * AUPE rereey coe Peis 


sf One4 76) 5m ) 


< 
Fe : > ‘UReeaT 
- 
oy P 
Df od Alen 
eae tacrss Qo aoes: 
i P, ) 
‘tantics DAN): Ylevres! 
5: SV ke ror] | 
. 7 2ykh 1¢ 
Ona Bis 
4 fc eeSon rf 35 
i { ? i} 
: nie it 
4 
} eon | 


Ae wines Iv ee evse 


Rw sos 5 


pea acpel byt ‘wa cana 


oem aud, 


ion toh Reyne 


ae tert Maine ary" ey solaaneek pe 2 x 
idea a any: vi 


237) IOLoesagey I ag ija 3 
“ue Shoe: ars. ais wth: a vet stain 

33) ee, A — chee eh pat vr 
Hots et shane Lame anasicar ae 
op {chpebes Dion! ciate rier, Bonet: 208 oul 
Bayes t, dseseto is ae ylation 
27g Hiway a ec kh  tevigey aes apoon 
naa? #2 Ags aleve ‘So nitihon, soda s 
7 inci Toe 23 Banoo oF eee _ 


IV RESULTS AND DISCUSSION 


A. VEGETATION ZONATION 

The forest vegetation of the Wapiti map area lies within two Forest 
Regions and four Forest Sections (Rowe 1972). Most of the area lies within 
the Mixedwood (B 18a), Lower Foothills (B 19a) and Upper Foothills (B 19c) 
Sections of the Boreal Forest Region. The remainder of the forest vegeta- 
tion lies within the East Slope Rockies (SAl) section of the Subalpine 
Forest Region. A small portion of the area is alpine tundra. 

Rowe's (1972) forest regions are mapped on a very small scale (1: 
6,336,000) and are delineated primarily by the distribution of the major 
tree species. The descriptions provided for the forest regions, especially 
the elevation criterion, have enabled an estimation of the extent of the 


component forest sections within the Wapiti map area. 


1. Mixedwood (B. 18a) 

The Mixedwood Section of the Boreal Forest Region accounts for 
approximately 34% (5100 1h of the Wapiti map area, primarily in the 
North at elevations less than 900 m and upon predominantly glaciolacustrine 
parent materials. The characteristic forest vegetation of the well drained 
uplands is a mixture of aspen and balsam poplar, white birch, white spruce 
and balsam fir, the last two species most common in old stands (Rowe 1972). 
The mixedwood in the Wapiti map area is broken up by lodgepole pine stand 
outliers of the Lower Foothills Section (B. 19a), though aspen-dominated 


forest is most prevalent. 


2s) Lower Foothills (B. 19a) 


The Lower Foothills Section of the Boreal Forest Region accounts 


for approximately 38% (5700 me) of the Wapiti map area primarily in a 
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band running across the centre of the area at elevations from 900 to 1200 m 
and within predominantly ground moraine of Continental origin. The 
characteristic forest vegetation is dominated by lodgepole pine with lesser 
amounts of aspen and balsam poplar in the wake of fire (Rowe 1972). In 
older stands the white spruce and black spruce are very often present. 


Balsam and subalpine fir show scattered distributions in the upland forest. 
os.) Upper Foothills’ (B.-19c) 


The Upper Foothills Section of the Boreal Forest Region represents 
approximately 21% or 3200 te of the Wapiti map area primarily in the south 
and south-west portions at elevations from 1200 to 1800 m and commonly 
upon colluviated till of both Continental and Cordilleran origin. The 
characteristic forest vegetation is dominated by lodgepole pine with black 
spruce, subalpine fir and white spruce-Engelmann spruce hybrids commonly 
present. The Upper Foothills Section is a transition between the Lower 
Foothills Section and the Subalpine Forest Region. The Upper Foothills is 


characterized by rolling topography and deep valleys. 


4. East Slope Rockies (SA. 1) 


The East Slope Rockies Section of the Subalpine Forest Region represents 
approximately 6% or 825 me of the Wapiti map area primarily in the south- 
west portion at elevations from 1800 to 2000 m (Plate 1) and on the steep 
slope of uplifted Mesozoic shales and sandstones. A thin veneer of 
Cordilleran till is often present. The characteristic forest vegetation is 


dominated by Engelmann spruce and subalpine fir. Lodgepole pine dominates 


younger stands especially after fire (Rowe 1972). The East slope Rockies 
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Plate 1 


Plate 2 


Alpine tundra in southwest corner of Wapiti map area. 
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Section grades into alpine tundra in a few localities. 


5. Alpine Tundra 


Alpine tundra occurs above timberline on the higher mountains at eleva- 
tions above 2000 m in the south-west of the Wapiti map area (Plate 2) 
representing approximately 1% or 90 ae Parent materials of the predom- 
inantly Regosolic soils are Mesozoic shales and sandstones. A thin till 
veneer is sometimes present. Forest vegetation is absent, and consequent- 
ly plots were not established in this zone, Characteristic species of the 
alpine tundra include Dryas hookeriana, Carex nardina, Cassiope tetragona, 


Kobresia bellardii and Silene acaulis. 


B. CLASSIFICATION OF THE FOREST VEGETATION 


The forest vegetation of 83L was classified into 15 forest types on 
the basis of the dominant tree species, floristic composition and by environ- 
ment as inferred from soils. The floristic classification was patterned 
after Braun-BlanquetsS methods as described by Mueller-Dombois and Ellenberg 
(1974) and after a Bray and Curtis (1957) ordination. The concepts of the 
forest types were developed both during field investigations and after the 
plots were sampled and the data analyzed. No attempt was made to restrict 
sampling to certain forest types nor to exclude certain forest types from 
sampling, though certain forest types are not well represented, particularly 
those at high elevations. Table 1 represents 15 concepts of forest 
ecosystems for the Wapiti map area. 

In order to keep a conceptual grouping useful, it should not include 


excessive variability but also should not be so restricted that stands 
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Table 1: 


Percent cover values of plant species with 
greater than 20% presence, arranged in order 


of decreasing presence. 
is in parentheses. 


Vegetation type number 
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fitting the central concept are seldom encountered under the specified 

set of habitat conditions. When vegetation is viewed as a continuum on 

the landscape (Whittaker 1967, 1973), it becomes apparent that there are 
transitions between conceptual units of classification. Because the 
vegetation of the Wapiti map area was sampled within several age classes 
and across gradients of elevation, soil moisture, climate and soils, there 
are many stands which may fit nearly as well with more than one vegetation 
type as recognized here, thus complicating a classification of the 
vegetation. Many authors avoid this dilemma by putting certain restrictions 
upon their sampling, for example by sampling only mature self-perpetuating 
stands. Only very young stands (less than 40 years) and very old stands 
were not sampled in the present study, as the objectives were to characterize 
the forest productivity and floristic composition as well as to ascertain 
successional relationships. The species composition of all the sampled 
plots is presented in Table 1, with transitional plots included with the 
group they resembled most. These groupings appear to be fairly homogeneous 
with respect to soil, other site factors and forest productivity. Most 

of these types characterize a major portion of one of Rowe's (1972) Forest 
Sections represented within the Wapiti map area, but others, notably the 
white spruce types appear to be controlled more by edaphic factors than by 
elevation (climate) and occur over a wide altitudinal range. A brief 
description of each of the forest types follows. Plot locations are given 


in Appendix 3 and plotted on Figure l. 
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Vegetation types of the Wapiti map area 


1. Populus tremuloides / Viburnum edule / Rubus pubescens 


( plots 19,23,26,27,45,46,85,86,87,90,100,106,122,123 ) 


This type represents the aspen forests in the northern sectors of 
83L (Plate 3). They occur at low elevations (610 to 960 m) on gentle 
north slopes (Classes 1 to 5; CSSC 1976). They are young (45 to 125 yr.), 
have lush shrub and herb understorys and often have plentiful white spruce 
regeneration. Constant species (over 80% presence in plots within type) 
which comprise most of the total plant cover, are Viburnum edule, Lonicera 
involucrata, Rosa acicularis, Petasites Palmatus, Rubus pubescens, Epilobium 
angustifolium, Pyrola asarifolia, Aster conspicuus, Maianthemum canadense, 
Cornus canadensis and Linnaea borealis. No seecies have mean cover greater 
than seven percent (Table 1). Aralia nudicaulis and Spiraea lucida are 
locally abundant. The moss layer has generally low cover (Table 1). 

The anticipated climax is the Picea glauca-Rubus pubescens-Maianthemum 
canadense type. The abundant white spruce regeneration and maturing white 
spruce trees in the understory (Table 2) suggest probable eventual dominance 
by white spruce. Forest productivity is variable in the aspen forests 
ranging from CLI class 3 to 5 (Appendix 9) as determined from the plot data. 

Soils are well to imperfectly drained and include in order of 
decreasing abundance, Gleyed Gray Luvisols, Orthic Gray Luvisols, Solonetzic 
Gray Luvisols and Eluviated Eutric Brunisols on fine-textured lacustro- 
till, lacustrine and alluvium overlying lacustro-till parent materials. 


Lodge and Donneliy are the predominant soil groups. 
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Table 2: Regeneration density and standard deviation, regeneration presence, 
ranked total volume of major tre¥species within the 15 vegetation types 
Vegetation l Species . 
Type No. Sw Sb Fa L A B B 
: : + 
1 Regeneration density~. 108-102 14246 2665367) saci65 | 19-59 
(n=14) Regeneration presence 71 14 64 29 14 
Ranked total volume 3 i 6 2 a 4 5 
2 Regeneration density 20 20 
(n=1) Regeneration presence 100 100 
Ranked total volume 2 3 i 4 
‘ : + + 
3 Regeneration density 121-229 189-528 1029-2860 Do-se4ue7=169 13-32 
(n=1l1) Regeneration presence 45 18 36 55 18 18 
Ranked total volume it 4 5 6 2 3 
‘ F + + + + ae 
4 Regeneration density 276-474 180-572 762-1456 34-80 2 3 
(n=l1l) Regeneration presence 64 27 45 18 Bil 
Ranked total volume 1 2 5 6 4 3 7 
: : + + + 
5 Regeneration density 88-110 S625) 978-866 56-95 18747 
(n=7) Regeneration presence 57 14 86 14 14 
Ranked total volume al 5 3 2 4 
: ; + + 
6 Regeneration density 395-596 56-95 
(n=7) Regeneration presence 43 14 
Ranked total volume 2 a 3 
: : = + 
zy Regeneration density 2730-611 25-35 
(n=2) Regeneration presence 100 
Ranked total volume al 
“f - + + +., +. + 
8 Regeneration density 230-463 135-362 399-828 10-21 133-253 97-197 
(n=13) Regeneration presence 38 23 54 vig) eiak Bh3k 
Ranked total volume 2 4 5 1 3 7 6 
; P - + = 
9 Regeneration density 49-67 630-1260 31-47 
(n=4) Regeneration presence HS 25 50 50 50 
Ranked total volume 3 al 2 
= ~ + ~ + 
10 Regeneration density 69622325 73-205 192=301 “23-33 29=—72 15-55 
(n=13) Regeneration presence 54 SB 46 46 15 8 
Ranked total volume = 2 5 iL + 6 
; +. +_ - + + + + 
1k Regeneration density 91-203 524-701 30-98 44-78 33-54 3-8 30-91 
(n=9) Regeneration presence 33 67 22 44 33 ll ll 
Ranked total volume 3 2 5 1 4 
+ + ~ + + 
12 Regeneration density S=33 238-654 422-1280 10-26 1-4 
(n=31) Regeneration presence 10 42 45 16 3 
Ranked total volume 3 2 4 1 5 
; om * + - 
13 Regeneration density 34-69 258-769 1036-1590 47-105 
(n=1l1) Regeneration presence Zi oT 82 18 
Ranked total volume 4 3 2 1 5 
+ 
14 Regeneration density (Sjshlh te) 7 
(n=2) Regeneration presence 100 
Ranked total volume 2 
15 Regeneration density 2S 3707 
(n=l) Regeneration presence 100 100 
Ranked total volume Sel 2 


1 Sw-white spruce, -Sb~biack spruce, Fa-subalpine fir, Lp 


Bp-balsam, poplar, B-white birch, Se-Engelmann spruce 
of stems less than 1.3 cm DBH 


2 Stems ha 


3 Percent of plots of a particular vegetation type in whol 


of tree regeneration occurs. 


-lodgevole pine, A-aspen, 


ch a particular species 
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2. Populus balsamifera / Rosa acicularis / Thalictrum venulosum 

DLO ELS.) 

The balsam poplar forests occur sporadically at low elevations (690 m) 
on alluvial floodplains which are still subject to periodic flooding. The 
herb understory is dense but moss cover is sparse. Constant species are 
Rosa acicularis, Rubus strigosus, Lonicera involucrata, Thalictrum venulosum, 
Equisetum arvense, Aster conspicuus, Mertensia paniculata, Calamagrostis 
Canadensis and Vicia americana. Stands of this type are young (ca. 70 yr.) 
and appear to be succeeding towards the Picea glauca / Rubus pubescens- 
Maianthemum canadense type. This trend is suggested by the many under- 
story species in common to the two types (Table 1) and the presence of 
white spruce seedlings and maturing white spruce trees in the understory 
(Table 2). Balsam poplar productivity is fairly high (class 3), but heart 
-rot in this species occurs at an early age. 


Soils are well drained Orthic Regosols on alluvial parent material. 


3. Picea glauca / Rubus pubescens-Maianthemum canadense 

( plots 24,25,44,47,48,91,94,95,117,118,124 «) 

These white spruce forests (Plate 4) occur at low to medium elevations 
(670 to 1220 m) on generally north-sloping (classes 1 to 6) sites. They 
are generally young (70 to 140 yr.) and have well developed shrub and herb 
understorys. Constant species are Lonicera involucrata, Rosa acicularis, 
Viburnum edule, Rubus pubescens, Maianthemum canadense, Mitella nuda, 
Cornus canadensis, Linnaea borealis and Petasites Ppalmatus. Cornus 
stolonifera, Alnus crispa, Alnus tenuifolia and Aralia nudicaulis are 
often evident in this type but are seldom seen in the other white spruce 


types. Hylocomium splendens is the predominant moss. This type can 
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Plate 3 Populus tremuloides/Viburnum edule/Rubus pubescens 
type (Type l). 


Plate 4 Picea glauca/Rubus pubescens-Maianthemum canadense 
type (Type 3). 
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commonly be seen on depressional sites within aspen forest, suggesting 
that succession advances faster on these sites. The abundance of white and 
black spruce seedlings and fir seedlings in some stands of this type (Table 2) 
should ensure perpetuation of this type as well as increased abundance 
of black spruce and subalpine fir in some stands. Forest productivity is 
good (classes 3 and 4) but appears to be less in the older forests. It is 
distinguished from the wetter Picea glauca/Equisetum eevee aloe nian 
splendens type by the presence of Maianthemum canadense, lower Equisetum 
arvense cover, the absence of Carex capillaris, less moss cover and by 
generally better drained soils. 

Soils are moderately well to imperfectly drained Orthic Gray Luvisols, 
Gleyed Gray Luvisols and Luvic Gleysols on alluvium over lacustro-till, 
lacustro-till, Continental and Cordilleran till. Donnelly, Snipe and Edson 


are the predominant soil groups. 


4. Picea glauca / Eqguisetum arvense / Hylocomium splendens 

(Solots 5,20,49;50,5/,59,61,/8,101,121,136 ) 

The white spruce-horsetail forests (Plate 5) occur at low to moderately 
high elevations (670 to 1450 m) on gentle (classes 1 to 4) generally north- 
facing slopes. They are young to moderately old (80 to 220+ yr.). White 
spruce and subalpine fir regeneration is common in many of the stands 
(Table 2). The understory is herb dominated and a dense Hylocomium splendens 
cover is present with lesser amounts of Ptilium crista-castrensis and 
Pleurozium schreberi (Table 1). Constant species include Rosa acicularis, 
Lonicera involucrata, Equisetum arvense, Petasites palmatus, Mertensia 
Paniculata, Mitella nuda, Cornus canadensis, Linnaea borealis anda Rubus 
pubescens. Carex capillaris,an indicator of the moist conditions of this 
type,is found in approximately one-half of the plots of this type. Forest 


productivity is variable (class 2 to 5) and appears to be less in the older 


forests. 
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Soils are poorly to imperfectly drained peaty Orthic Gleysols, peaty 
Luvic Gleysols and Orthic, Luvic, and Rego Gleysols on Continental till, 
alluvial sand and lacustrotill parent materials. Snipe, Smoky and 


Gunderson are the predominant soil groups. 


5. Picea glauca / Hylocomium splendens 

(Pe pLOts4/ 710,126, 12752267) 1297133° ©) 

The white spruce-feathermoss forests (Plate 6) occur at low to 
moderately high elevations (610 to 1280 m), often along valleys of major 
rivers on old river terraces, on generally north sloping (classes 1 to 6) 
sites. They are young to peat aeel old (70 to 170 yr.), have sparse 
shrub understorys, a herb-dwarf shrub stratum dominated by Cornus canadensis 
and Linnaea borealis and a dense cover of Hyiocomium splendens and Pleurozium 
schreberi.. Lodgepole pine forms a large proportion of the basal area in 
some stands but is being replaced by white spruce. Subalpine fir regenera- 
tion is abundant in most of the stands (Table 2), suggesting eventual 
co-dominance of subalpine fir with white spruce. Constant species are 
Cornus canadensis, Linnaea borealis, Pyrola secunda, Hylocomium splendens 
and Pleurozium schreberi. This type is stable and has moderate forest 
productivity (class 4 to 5). 

Soils are rapidly to well-drained Orthic and Cumulic Regosols, 
Brunisolic Gray Luvisols and Orthic Gray Luvisols on generally coarse- 
textured alluvial sands and gravels, outwash gravels, and coarse-textured 


Cordilleran till. Jarvis and Robb are the predominant soil groups. 
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6. Picea mariana / Ledum groenlandicum / Equisetum arvense 

fe DLocs 15,321,633, 6o7Oleioyis> <) 

The black spruce forests occur at low tc mid-elevations (840 to 
1220 m) on gently sloping (classes 1 to 4) sites with a northerly aspect, 
and may occur in association with the Picea mariana / Ledum grocnlandicum / 
Rubus chamaerorus type. The Slight slope.of stands of this type propably 
allows more oxygen to move through the peaty soil than the black spruce- 
labrador tea-cloudberry type. The stands are young to fairly old (75 to 
180 yr.) and are usually in climax condition, with abundant black spruce 
regeneration (Table 2). The herb understory is often species rich. Constant 
species are Ledum groenlandicum, Vaccinium vitis-idaea, Equisetum arvense, 
Petasites palmatus, Carex capillaris, and Cornus canadensis. The moss 
stratum is dominated by Hylocomium splendens, Pleurozium schreberi and 
lesser amounts of Ptilium crista-castrensis and Sphagnum spp. Though 
maetal trees often attain merchantable size (such as occasional white 
spruce) in this type, overall productivity is rather low (classes 5 and 6). 

Soils are poorly drained Terric Mesisols, peaty Rego Gleysolis and Humic 
Luvic Gleysols on organic and occasionally till parent materials. Kenzie 


and to a lesser extent, Smoky are the predominant soil groups. 


7. Picea mariana / Ledum groenlandicum / Rubus chamaemorus 

(se plots al4;29 2%) 

This type represents the black spruce bog vegetation (Plate 7). The 
bogs occur at low to mid-elevations (915 to 1070 m) in depressions with 
impeded drainage on level sites with hummocky microtopography. These 
open forests are often over 200 years old and can be considered climax. 
The well developed shrub layer is dominated by Ledum groenlandicum. The 


herb-dwarf shrub understory is dominated by Vaccinium vitis-idaea, Rubus 
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chamaemorus and Oxycoccus microcarpus. Sphagnum Spp. are abundant. Tree 
cover is sparse and productivity is very low (class 7) and can be considered 
non-merchantable. 

Soils are poorly drained Typic Mesisols and Fibrisols on MOSS peat 


parent materials. Kenzie is the predominant soil unit. 


8. Pinus contorta / Viburnum edule / Rubus pubescens 

mpstots.o8, 71792793, 97,99 ,9103,104,105, 107,108,116, 12080) 

This is a relatively common forest type in 83L on north-sloping (classes 
1 to 6) sites at mid-elevations (880 to 1040 m). It is characterized by 
young (80 to 90 yr.) lodgepole pine stands of fire Origin and often 
plentiful white spruce, black spruce, subalpine fir and paper birch 
regeneration (Table 2), indicating possible eventual dominance of these 
Species. The occasional occurrence of Oplopanax horridum in this type 
(Appendix 7), and the lush herb understory seem to indicate high precipi- 
tation. Constant species include Viburnum edule, Lonicera involucrata, 
Rosa acicularis, Rubus pubescens, Cornus canadensis, Linnaea borealis, 
Pyrola asarifolia and Epilobium angustifblium. Pleurozium schreberi 
dominates the moss layer. Though not restricted to this type, Sorbus 
scopulina, Spiraea lucida, Maianthemum canadense and Gymnocarpium 
dryopteris may be considered to be indicators of these forests. This 
type is similar to the Pinus contorta / Spiraea lucida / Cornus canadensis 
type but appears to be more moist as indicated by the higher cover of 
Viburnum, Sorbus, Rubus and Gymnocarpium. This type also has affinities 
with the Picea glauca / Rubus pubescens-Maianthemum canadense type towards 
which it might be expected to succeed. Forest productivity is good 


(classes 3 and 4). 
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Pinus contorta/Alnus crispa/Cornus 


canadensis type (Type 10). 


Plate 8 


/Ledum groenlandicum/Rubus 


chamaemorus type (Type 7). 
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Soils are rapidly to imperfectly drained and include Gleyed Gray Luvisols, 
Eluviated and Orthic Eutric Brunisols, "bleached" Gray Luvisols (Brunisolic 
Gray Luvisols with Ael horizon of high color value above a darker Ae2 horizon) 
and Orthic Gray Luvisols on Continental till, sandstone and shale bedrock 
and aeolian sand parent materials. Edson, Torrens and Copton are the 


predominant soil grouos. 


9. Pinus contorta / Spiraea lucida / Cornus canadensis 

(plots 437/20 2,)L09 piss |) 

This type is not common and occurs at low elevations (915 to 1265 m) 
on sloping topography (Classes 4 to 7) in association with the Pinus 
contorta / Viburnum edule / Rubus pubescens type. It is characterized by 
young (less than 125 yr.) lodgepole pine stands of fire origin with a lush 
herb-dominated understory and regenerating aspen, paper birch and some 
black spruce. Constant species include Betula papyrifera, Rosa acicularis, 
Spiraea lucida, Cornus canadensis, Linnaea borealis, Arnica cordifolia, 
Lycopodium annotinum, Lycopodium compianatum, Pyrola asarifolia and 
Pleurozium schreberi, the dominant moss (Table 1). Forest productivity 
is moderate (classes 4 and 5). 

Soils are well to moderately well drained Orthic Gray Luvisols and 
Brunisolic Gray Luvisols on Continental till, Cordilleran till and alluvial 
sand over bedrock parent materials. Soil groups include Mayberne, Edson, 


Robb and Lodge over unconsolidated sandstone bedrock. 


10. Pinus contorta / Alnus crispa / Cornus canadensis 
( “plots °9721722),28;,37,62766)77;83 ,109; 113,114,130; ) 
This type is fairly common and is characterized by young (less than 


100 yr.) lodgepole pine stands of fire origin at moderate elevations 
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(840 to 1280 m). Well developed shrub and herb understorys are present 
(Plate 8), and many stands have abundant white spruce regeneration (Table 2). 
Constant species include Alnus crispa, Vaccinium caespitosum, Cornus 
canadensis, Linnaea borealis, Arnica cordifolia, Petasites palmatus, 
Calamagrostis canadensis, Pyrola asarifolia, Equisetum arvense and 
Pleurozium schreberi,the dominant moss. The shrub Rubus parviflorus has 
sporadic occurrence in this type. 

Adequate tree seedling regeneration following clearcutting might be 
a problem on this type where great increases in competing Alnus and 
Calamagrostis cover such as those observed by Corns and LaRoi (1976) for 
the area north of Edson, Alberta could be anticipated. Forest productivity 
is moderate (classes 4 and 5). 

Soils are moderately well to imperfectly drained Orthic and Gleyed 
Gray Luvisols, Brunisolic Gray Luvisols and "bleached" Gray Luvisols on 
Cordilleran or Continental till and shale and sandstone bedrock. Edson, 


Marlboro and Torrens are the predominant soil groups. 


ll. Pinus contorta / Ledum groenlandicum / Pleurozium schreberi 

( plots 12,13,55,60,70,84,88,89,96 ) 

The lodgepole pine - Labrador tea - feathermoss forests (Plate 9) are 
common at moderate elevations (745 to 1220 m) and are characterized by 
young to fairly old (80 to 180 yr.) lodgepole pine stands of fire origin 
on sloping topography (classes 1 to 5). Black spruce is often present 
(Table 2) and may eventually assume dominance (eg plot 55). Labrador tea 
is usually abundant but dwarf shrub and herb cover is usually sparse 
(Table 1). The moss layer is dominated by a dense cover of Pleurozium 


schreberi, Hylocomium splendens and Ptilium crista-castrensis with lesser 


amounts of the lichen, Peltigera aphthosa which reaches its highest cover 
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inthis: type. Constant species include Ledum groenlandicum, Vaccinium 
vitis-idaea, Cornus canadensis, Linnaea borealis, Pleurozium schreberi 
and Peltigera aphthosa. This type is similar to the lodgepole pine- 
black spruce-tall bilberry type but seems to occur in areas of lower precipi- 
tation at generally lower elevations. It lacks the tall bilberry which is 
often replaced by the blueberry (Vaccinium myrtilloides). Forest productivity 
is moderate (classes 4 and 5). | 

Soils are moderately well to imperfectly drained Orthic Gray Luvisols, 
Brunisolic Gray Luvisols, and "bleached" Gray Luvisols on Continental and 
Cordilleran till, lacustro-till and bedrock parent materials. Soil groups 


include Edson, Snipe, Torrens, Copton, Sneep, and Heart. 


12. Pinus contorta / Picea mariana / Ledum groenlandicum / Vaccinium 
membranaceum 


( plots 1,3,4,6,8,11,17,30,34,36,38,39,51,52,64,66,69,/2, 74-76, 
19-62, 96, 0K 02,115 713757) 


The lodgepole pine-black spruce-Labrador tea-tall bilberry forest type 
(Plate 10) is more extensive than any of the others in the 83L area and 
occurs on gently sloping (classes 1 to 5) sites of variable aspect from 
low to relatively high elevations (840 to 1465 m). It is characterized 
by young to fairly old (65 to 190 yr.) lodgepole pine and black spruce 
stands of fire origin. Black spruce forms a tree understory layer of 
approximately the same age as the pine. Black spruce and subalpine EGE 
regeneration is often abundant (Table 2), indicating probable eventual 
succession to these species. Ledum often forms a dense low shrub understory, 
and herb cover is moderate. Constant species include Ledum groenlandicum, 


Vaccinium membranaceum, Vaccinium vitis-idaea, Cornus Canadensis and Linnaea 


borealis. A dense feathermoss cover of Pleurozium schreberi and Hylocomium 


splendens is usual. This type is transitional to the Pinus contorta / 
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Pinus contorta/Ledum groenlandicum/Pleurozium schreberi 


type (Type ll). 


Plate 9 
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Pinus contorta/ Picea mariana/Ledum groenlandicum 


-Vaccinium membranaceum type (Type 12). 


Plate 10 
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Ledum groenlandicum / Pleurozium schreberi type at lower elevations. At 
the upper limits of type 12, Rubus pedatus is common, and Menziesia glabella, 
Rhododendron albiflorum, Tiarella trifoliata, and Arnica latifolia are 
sporadic in occurrence. Forest productivity is moderate (class 5 with a 
few exceptions). 

Soils are moderately well to imperfectly drained Orthic Gray Luvisols, 
Brunisolic Gray Luvisols and "bleached" Gray Luvisols. Edson, Mayberne and 


Marlboro are the predominant soil units. 


13. Pinus contorta / Menziesia glabella / Rubus pedatus 

(PLO CS) 32,535,907 40-42 75375471297 102,134 ") 

The lodgepole pine-false azalea type (Plate 11) is a subalpine forest 

type restricted to the south and west portions of 83L at high elevations 
(1220 to 1555 m) and is characterized by young to fairly old (70 to 180 yr.) 
lodgepole pine stands of fire origin. Black and white spruce and subalpine 
fir regeneration is commonly abundant in ‘the older stands (Table 2) indicating 
probable eventual succession to these species. A shrub understory of 
Menziesia glabella and often Rhododendron albiflorum, another characteristic 
species of subalpine forests is present. Herb cover is moderate. The 
subalpine herb Veratrum escholtzii occurs sporadically within this type. 
The predominantly shrub understory could provide competition with tree 
seedlings following clearcutting on this type. Constant species include 
Menziesia glabella, Ledum groenlandicum, Vaccinium membranaceum, Vaccinium 
vitis-idaea, Rubus pedatus, Cornus canadensis, Linnaea borealis, Lycopodium 
annotinum, Pleurozium schreberi, Ptilium crista-castrensis and Hylocomium 
Splendens. Forest productivity is low (classes 5 and 6). 

Soils are well to moderately well drained Brunisolic Gray Luvisols 


and "bleached" Gray Luvisols on Continental till and thin Cordilleran till 
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over bedrock parent materials. Mayberne is the predominant soil group 


with lesser occurrences of Marlboro, Torrens, and Robb. 


14. Pinus contorta / Vaccinium myrtilloides / Cladonia spp. 

( plots 16,56 ) 

The lodgepole pine-blueberry-lichen type was not common and was only 
sampled twice. It is restricted to well drained, coarse textured soils. 
The tree canopy is open, and black spruce regeneration is usually present. 
This type is closely related to the lodgepole pine-Labrador tea-feathermoss 
type and is characterized by a similar climate and elevational range. Herb 
and moss cover is sparse but lichen cover may exceed 30%. Constant species 
include Ledum groenlandicum, Vaccinium myrtilloides, Vaccinium vitis-idaea, 
Arctostaphylos uva-ursi, Cornus canadensis, Pleurozium schreberi, and 
Cladonia spp. Forest productivity is low (class 6), probably due to 
limited soil water holding capacity. 

Soils are well drained Orthic Humo-ferric Podzols and Eluviated Eutric 
Brunisols on outwash sand and alluvial sand parent materials. Blackmud is 


the predominant soil group. 


15. Picea engelmannii-Abies lasiocarpa / Menziesia glabella 

Cipiors2e) 

The Engelmann spruce-subalpine fir-false azalea forests (Plate 12) 
form a climax type which occurs on steep (classes 5 to 7) north-facing 
slopes at high elevations (above 1670 m) in the south-west corner of the 
Wapiti map area. Menziesia may form a fairly dense shrub understory but 
herb and low shrub cover is generally sparse. Constant species include 
Menziesia glabella, Phyllodoce empetriformis, Vaccinium membranaceum, 


Rubus pedatus, Pedicularis bracteosa, Cornus canadensis, Lycopodium annotinum 
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and Arnica latifolia. Tree growth is slow (classes 5 and 6) and stands 
are usually not suitable for commercial use. The type is species poor and 
would show a slow recovery after disturbance. 

Soils are moderately well to imperfectly drained Orthic Gray Luvisols 


on Cordilleran till. Robb and Copton are the predominant soil groups. 


C. GEOGRAPHIC RELATIONS AND OTHER STUDIES 


The flora of the Wapiti map area is predominantly Boreal, though it 
contains Cordilleran and Pacific elements (Moss and Pegg 1963), notable 
species with primarily Cordilleran distribution are Picea engelmannii, 
Abies lasiocarpa, Pinus contorta Var. latifolia, Menziesia glabella, 
Rhododendron albifiorum, Sorbus scopulina, Spiraea lucida, Veratrum 
eschscholtzii, Rubus pedatus, Vaccinium membranaceum, Phyllodoce 
empetriformis and Arnica cordifolia. Species with primarily Pacific 
distributions include Oplopanax horridum, Rubus parviflorus, Sambucus 
pubens and Tiarella trifoliata. 

Comparisons between the forest types as described here and those 
described by other workers for the western boreal forest are in large part 
based upon subjective assessments of floristic composition and dominance 
due to a lack of published quantitative data. Early forest classifications 
such as those by Raup (1933, 1934, 1946) and Moss (1953, 1955, Moss and 
Pegg 1963) are general and in a single grouping may include several of the 


vegetation types described here. For example, the "upland spruce forests" 


referred to by Raup (1933, 1934, 1946) would include types 3, 4 and 5 of the 


present study. 


In the following discussion, each of the forest vegetation types or 


the Wapiti map area are compared briefly with similar vegetation described 
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by other authors for the western boreal forest. 

The Populus / Viburnum / Rubus type is included within the aspen 
poplar consociation of the poplar association (Moss 1953) and bears 
strong resemblance to the Populus / Aralia / Linnaea ecosystem of Kabzensc 
et al. (1976) described for central Saskatchewan,and appears to closely 
resemble the aspen forests at Candle Lake Saskatchewan described by Swan 
and Dix (1966; Dix and Swan, 1971). The aspen forests of the Wapiti map 
area bear somewhat less similarity to the Populus tremuloides variant of 
the Ptilium crista-castrensis / Gymnocarpium dryopteris / Abies lasiocarpa- 


Picea glauca association of Wali and Krajina (1973) and to the Mitella- 


(1945) described for north-eastern British Columbia. 


cicularis / Thalictrum venulosum type 


hy 


The Populus balsamifera / Rosa 


, 


corresponds to the balsam poplar consociation of the poplar association 
(Moss 1953) and is similar to the balsam poplar forest vegetation described 
by Dix and Swan (1971) for the Candle Lake area of Saskatchewan. 
Vegetation bearing close resemblance to the Picea / Rubus-Maianthemum 
type has been described by authors in Alberta and Saskatchewan. This 
vegetation corresponds to the white spruce-shrub-herb faciation of the 
white spruce association described for Alberta by Moss (1953); is included 
within the Lonicera / Rubus-Lathyrus group described by LaRoi (1967) for 
western Canada and within the Picea-Abies / Viburnum-Zylocomium community 
type of Achuff and LaRoi (1977) described for the highlands of northern 
Alberta and appears very similar to the white spruce-sarsaparilla associa- 
tion of Lesko and Lindsay (1973). There is also close Similarity to the 
Picea glauca / Populus / Cornus / Rubus ecosystem of Kabzems et al. (1976) 
described in central Saskatchewan and to the white spruce vegetation 


described by Swan and Dix (1966) and Dix and Swan (1971) in the Candle Lake 
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area of Saskatchewan. The most similar described vegetation to the west 

of the Wapiti area is the Ptilium crista-castrensis / Gymnocarpium dryopteris / 
Abies lasiocarpa-Picea glauca association described by Wali and Krajina (1973) 
for north-eastern British Columbia. 

The Picea glauca / Equisetum arvense / Hylocomium splendens type 
corresponds to the feathermoss faciation of the white spruce association 
described for Alberta by Moss (1953); is included within the Picea-Abies / 
Viburnum / Hylocomium community type of Achuff and LaRoi (1977), described 
for the highlands of northern Alberta and appears very similar to the 
white spruce / horsetail association of Lesko and Lindsay (1973) described 
in west-central Alberta and appears comparable to the Picea glauca-Equisetum 
ecosystem described by Kakzems et al. (1976) for central Saskatchewan. 
Somewhat similar vegetation described to the west of the study area includes 
the Gymnocarpium dryopteris / Oplopanax horridum / Abies lasiocarpa-Picea 
glauca association described by Wali and Krajina (1973) for north-eastern 
British Columbia and the Mitella-Equisetum "colony" of Kujala (1945) 
described in eastern British Columbia. 

The Picea glauca / Hylocomium splendens type is also included within 
the feathermoss faciation of the white spruce association described for 
Alberta by Moss (1953) and appears to be analogous to the white spruce- 
feathermoss association of Lesko and Lindsay (1973) and to the Picea glauca / 
Pleurozium type of Kabzems et al. (1976) on well drained sites in central 
Saskatchewan. 

The Picea mariana / Ledum groenlandicum / Equisetum arvense type is 
included within the black spruce bog forest described for Alberta by Lewis 
et al. (1928); within the black spruce-feathermoss association of Moss (1953) 
and within the Vaccinium vitis-idaea-Fragaria-Maianthemum stand group in 


the Picea mariana / Rosa-Ribes triste / Mitella-Mertensia group described 
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for western Canada by LaRoi (1967). The black spruce-Labrador tea-horsetail 
forests in the Wapiti area also seem to correspond closely to the Picea 
mariana / Pleurozium / Hylocomium ecosystem described by Kabzems et al. 
(1976) and to the Sphagnum nemoreum-Pleurozium schreberi-Ptilium crista- 
castrensis-Hylocomium splendens-Cornus canadensis / Picea mariana associa- 
tion described by Wali and Krajina (1973). 

The Picea mariana / Ledum groaenlandicum / Rubus chamaemorus type is 
also included within the black spruce bog forest of Lewis et al. (1928); 
within the black spruce-peat moss association of Moss (1953) and is included 
within the black spruce-peat moss bog association of Lesko and Lindsay (1973). 
This vegetation type corresponds to the Picea mariana / Ledum / Sphagnum 
ecosystem of Kabzems et al. (1976). 

The Pinus contorta / Viburnum edule / Rubus pubescens type and the 
Pinus contorta / Spiraea lucida / Cornus canadensis type do not appear to 
have counterparts described for the western Boreal Forest region but the 
lodgepole pine forests of the Wapiti area have in general, floristic composi- 
tion similar to those described by Newsome and Dix (1968) for the Cypress 
Hills of Alberta and Saskatchewan. The Pinus contorta / Alnus crispa / 
Cornus canadensis type appears most similar to the Pinus / Pleurozium / 
Lycopodium ecosystem of Kabzems et al. (1976) with the major difference 
being that Pinus banksiana seems to occupy soils and aspects in Saskatchewan 
similar to those occupied by Pinus contorta in the Alberta foothills. 

The Pinus contorta / Ledum groenlandicum / Pleurozium schreberi type 
most closely resembles the Pinus / Picea mariana-Pleurozium ecosystem of 
Kabzems et al. (1976) with the exceptions that in Saskatchewan, jack pine 
substitute for lodgepole pine and black spruce appears to be more abundant 
in the comparable community. The Pinus contorta / Vaccinium vitis-idaea 


type of Kujala (1945) also appears to be comparable. 
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The Pinus contorta / Picea mariana / Ledum groenlandicum / Vaccinium 
membranaceum type bears most resemblance to the Vaccinium membranaceum type 
described by Kujala (1945) for eastern British Columbia which is much more 
species rich than its Alberta ccunterpart. The Pinus contorta / Menziesia 
glabella / Rubus pedatus type is primarily subalpine rather than boreal in 
its distribution and is most comparable to the Vaccinium membranaceum type 
described by Kujala (1945) but is not as species rich as the British 
Columbia community. Similar vegetation is common in the mountains of 
Jasper and Banff to the south (Corns and Kojima 1976). In time, the 
lodgepole pine will be replaced by Engelmann spruce and subalpine fir,as it 
has been in some of the high elevation forests of the Wapiti area. The 
Engelmann spruce-Subalpine fir-false azalea type is most similar to the 
subalpine Menziesia-Vaccinium membranaceum "colony" and the subalpine 
Tiarella trifoliata-Rubus pedatus / Rhododendron albiflorum "colony" 
described by Kujala (1945) but is not as species rich as the British 
Columbia representatives. 

The Pinus contorta / Vaccinium myrtilloides / Cladonia spp. type has 
a similar counterpart in the Cladonia gracilis-Arctostaphylos uva-ursi- 
Vaccinium myrtilloides / Pinus contorta association described by Wali and 
Krajina (1973). It is also comparable to the drier stands of the Pinus / 
Vaccinium vitis-idaea / Pleurozium ecosystem described by Kabzems et al. 
(1976) except that in Saskatchewan jack pine substitutes for lodgepole 


pine on rapidly drained sandy sites. 


D. INFLUENCE OF ELEVATION, ASPECT AND INTERNAL SOIL DRAINAGE ON VEGETATION 
TYEE 


A three variate representation of vegetation type relationships in 
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Figure 4: Comparison of the 15 vegetation types 
of the Wapiti map area with respect to 


elevation, slope aspect and internal 
soil drainage. 
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the form of a cube helps to clarify some of the site differences among 
vegetation types. Figure 4 illustrates such relationships where soil 
drainage is represented on the X axis, elevation on the Y axis and slope 
aspect (arranged in order of decreasing insolation received) on the Z 

axis. Warm, low, dry situations are represented in the bottom left corner 
and cool, high, moist situations are represented in the upper right corner. 
The mean elevation and drainage and the modal aspect were plotted for each 
of the 15 vegetation types sampled. 

It is apparent from an examination of Fig. 4 that several of the 
vegetation types occur on predominantly northerly aspects. Others seem 
to have no preference for slope aspect and were thus plotted on a neutral 
(westerly) aspect. Why some vegetation types do not show more pronounced 
affinities for southerly aspects is not apparent. 

As predominant elevation, soil drainage and slope aspect have previously 
been described for each of the vegetation types, only pronounced similarities 
in these factors between floristically distinct vegetation types will be 
discussed in this section. It is expected that the three-dimensional 
presentation should show for example, successional stages of a particular 
type in a similar position on the graph, assuming that the successional 
stages had initially been recognized as separate vegetation types. 

On the basis of elevation alone, the aspen (Type 1) and balsam poplar 
forests (Type 2) are separated at low elevations from the other vegetation 
types. The Engelmann spruce forests (Type 15) and the lodgepole pine-false 
azalea type (Type 13) are separated from the other types at high elevations. 
Most of the other lodgepole pine, white spruce and black spruce types 
cannot be separated on the basis of elevation alone. 

When soil drainage alone is considered, the black spruce bog forests 


(Type 7), the black spruce-horsetail forests (Type 6) and the white spruce- 
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horsetail forests (Type 4) can be separated from the other forest types. 
These three spruce types are similar to each other with respect to elevation 
and slope aspect also. When other factors are considered, it is apparent 
that the black spruce bogs are usually nearly level and that the water 
table is very high. Reducing conditions are present at depth due to the 
anaerobic situation. The white spruce-horsetail type (Type 4) and the black 
spruce-labrador tea-horsetail types (Type 6) both occur on gentle northerly 
slopes but the latter occurs on predominantly organic soils as opposed to 
mineral in the former, and with generally poorer drainage. The floristic 
composition of these sites further separates them. The white spruce-feather- 
moss type (Type 5) is separated from the other white spruce types by the 
higher elevations and better drained soils. The Picea glauca / Rubus 
pubescens-Maianthemum canadense type (Type 3) is moderately well to im- 
perfectly drained and intermediate with respect to internal drainage 
between the other white spruce types. On Fig. 4, Type 3 is similar to the 
Pinus contorta / Viburnum edule / Rubus pubescens type (Type 8) with 
respect to elevation, soil drainage and aspect, type 3 having slightly 
more restricted soil drainage than type 8. The ages of stands of the two 
types overlap. Species composition of the two types is similar and it 
appears likely that type 8 would succeed towards vegetation similar to 
type 3 under the shade of the lodgepole pine canopy. 

Type 8 also has affinities with the Pinus contorta / Spiraea lucida / 
Cornus canadensis type (Type 9), from which it differs by having soils 
with slightly more impeded drainage than type 9. Type 9 might also be 
expected to succeed toward the Picea glauca / Rubus pubescens-Maianthemum 
Canadense type (Type 3) as might the aspen type (Type 1). Successional 
relationships may be hypothesized with the aid of Figure 4 and will be 


discussed in more detail in a later section. 
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Lodgepole types 10, 11 and 12 show differences in soil drainage and 
mean elevation (Fig. 4). Other differences are described in the type 
descriptions. They develop independently and are not related successionally 
to each other. The lodgepole pine-blueberry-lichen type (Type 14) most closely 
resembles the white spruce-feathermoss type (Type 5) in terms of elevation, 
soil drainage and aspect but floristically most closely resembles the 
lodgepole pine-black spruce-Labrador tea-tall bilberry type (Type 12). The 
lodgepole pine-blueberry-lichen type (Type 14) is likely an edaphic climax 
due to the difficulty encountered by spruce and fir regeneration in soil 


with low moisture reserves. 
Bie SUCCESSIONAL RELATIONSHIPS 


The forests of the Wapiti map area are in general very young, and 
with a few exceptions are first generation stands which have arisen in the 
last 200 years. The 15 forest types described have for the most part 
Originated independently and have not succeeded from one of the others. 

In general most of the area appears to be succeeding toward a white spruce 
Or mixed spruce-subalpine fir climax, on the basis of the presence of these 
species in the understory layers of many of the present forests (Table 2). 

In the mixedwood of the north-eastern portion of the map area, white 
MO aominaten stands often occur in depressions within largely aspen- 
dominated forests. It thus appears that the cool, north-facing slopes 
and sites with impeded drainage seem to succeed toward their potential 
"Climax" or more stable condition before those in better drained locations. 

Figure 5 suggests possible successional relationships between the 
The dotted lines indicate extremely slow succession, of the 


forest types. 


order of hundreds of years. The solid lines indicate more rapid succession, 
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of the order of 100 to 300 years from present. In general, the spruce- 
dominated vegetation is more stable than the lodgepole pine-dominated vegeta- 
tion, though it is apparent that succession within some of the lodgepole 

pine types is proceeding very slowly. For example, the nearly 200-year 

Old Pinus contorta/Menziesia glabella/Rubus pedatus forests are very similar 
physiognomically and floristically to much younger (less than 100 yr.) 

stands on similar sites. 

In general it appears that there are at least two successional trends 
from the present forest vegetation to climax: lodgepole pine, aspen and 
balsam poplar to white spruce; and lodgepole pine and upland black spruce 
to mixed spruce-subalpine fir forests (Figure 5). The Pinus contorta / 
Ledum groenlandicum / Pleurozium schreberi type and the Pinus contorta- 
Picea mariana / Ledum groenlandicum / Vaccinium membranaceum type appear 
to be succeeding towards forest vegetation dominated by Picea glauca, P. 
engelmannii (and hybrids of these two), P. mariana and Abies lasiocarpa. 
Lodgepole pine is absent or present only as a relic. The understory is 
floristically poor in the spruce-fir-feathermoss type with the most prominent 
feature being a nearly complete cover of Hylocomium splendens. The mature 
or climax Picea spp. / Abies lasiocarpa-Hylocomium splendens forests are 
sporadic in occurrence and were not sampled due to their "“overmature" 
appearance, being unsuitable for stem analysis sampling. 

The other major successional trend appears to be from types 10, 8 and 
1 to Picea glauca/Rubus pubescens-Maianthemum canadense (Type 3) forests 
based upon understory similarities and often abundant white spruce regenera- 
tion in the pine and aspen forests. Succession appears to be proceeding at 
different rates, being most rapid on the moister sites. Types 9 and 2 
seem to be heading towards type 3 also, but much more slowly. It is 


unlikely that the balsam poplar forests will succeed to white spruce as 
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long as frequent flooding occurs, covering young spruce seedlings with 
fresh alluvium. 

Other successional trends are less evident. The black spruce-Labrador 
tea-horsetail forests (Type 6) might be expected to succeed towards the 
white spruce-horsetail-feathermoss forests (Type 4) if drainage in the 
former was improved by man or on a much longer time scale, by regional 
lowering of the water table. White spruce commonly forms a minor component 
in Type 6 forests (Table 2), Succession from black spruce bog vegetation 
(Type 7) towards the black spruce-Labrador tea-horsetail forests (Type 6) 
and from white spruce-horsetail-feathermoss forest (Type 4) to white spruce- 
dewberry-two leaved Soloman's seal forests (Type 3) might similarly occur 
with increased soil drainage in type 6 in the former case and type 4 in the 
latter. A certain amount of drying of the site might be expected to occur 
naturally since the growing trees of the maturing forest seem likely to 
demand more water for growth, and transpiration, though quantitative data 
are scarce (Satterlund 1972 p. 141). 

The lodgepole pine-green alder-bunchberry type (Type 10) might go in 
one of two directions (Fig. 5); towards Type 3 on imperfectly drained sites 
or towards the better drained white spruce feathermoss forest (Type 5) if 
the site was currently well drained. The former appears more likely since 
most of the pine-alder forests tend to be more imperfectly- than well drained. 
Type 5 in time, with maturation of the already abundant Abies lasiocarpa 
seedlings present in many stands (Table2 ) will probably succeed towards 
the Picea spp.-Abies lasiocarpa / Hylocomium splendens type. 

In the past, due to the frequency of large devastating fires, climax 
forest development was achieved only in isolated areas. Today, with greatly 
improved access and modern fire suppression techniques, fires are much 


smaller and less intense than those of 75 to 85 years ago when most of the 


oe 


Ae ne pet i Yen 

Sobpudat~ssrigs pone opi donkiin pin un staege fa 
iin sbhosug? Deaaaue 62 ba pou ba A houa oo - 
ede ni dpeoiasy aii Gey spent | , | 


i 
= iA: on ee 5 


isvolper yd gaia Sho, Re) veh A wie ee an 9B 
Fipmawos Tote & Keates ‘sees aan Slee wedded > 
4) de Souy, -sriet” Setas sen spepmi tied hut sp aeegReR, 8. ahetwt) 7 
(O ogy?) admeret Tanste a ed raheadaronatn amd eid a 4 
—wiuyse wiide om, sb eeyb) oeerves cheat etaODInf Ladd one ROONgm | 
swanQ Ylaclinie, tein (ft seul), eetaset ison 2'nemoto® bevesd, 2 
Pom + om Lins wk Rombod ate ak eS ways it peedberh 
rojo Ot Hetosgnve we Sale atin ae? Beenive Ro ITI 
Of \hertli dees deethh otbmetied saa ge Baad psiwoae had | 
Stab avlse doabep equa thie shan Ervin. aed deveragy eT oy 
| (hee BP VERE Crd ae: 
al oo Ichbin 19? eget) ea Cp Rete” Hest entg: 
2et1h beciead yioob een week eoee abraeus) we ‘Bid innioacian 
"4 5 coy) J8sg04 eer pate eirgh- gu tie Qankeah aed ied Sas Olay 
quire vlethl o9te @aewe Reeve wee /Hediet trap aieeertes on i" 
‘yea. 01 tio seds~y Ree hiegie sna et Gi bees eeaeees vallansnhy ge 
ee ee eT “a 
sbievny cossaue Qitdéng £64 Selden ated Qode OP Soseene 
Sah Aig MTeoR | eetawod tant “SekeNS pee 

“ecile ,€o.59 erizgeiedasd. apsel Se ehabubes? add bs got (Sead ons Ar : 
yleswog Wiiw , eho. ieeesd Bor) al eine Sovelise seu thdeuo level Sadiod | 
Tob? 206 sec). ~owehiviiogg Aoregewgwe @1i) asebom torias ieee Maen 
| df To tion cette ote Gragg ap ad SS Fe cepans Gn) asietnd oo 


63 


lodgepole pine stands in the Wapiti map area originated (Table 3). How- 
ever, it is still unlikely that climax forest conditions will be achieved 
over large areas even in the absence of fire, due to logging and pulp- 
cutting activities, which are likely to result in successional series 
somewhat different from those after fire. 

If succession was permitted to continue, vegetation similar to that 
depicted in Figure 5 might be expected to develop, with a number of diverse 
early successional stages on different physical sites converging towards a 
small number of climax or mature forest types in the course of hundreds of 


years under the influence of regional climate, sensu Clements (1916 et seq.). 


F. ORDINATION OF STANDS 


The purposes of this section are to describe the vegetational 
interrelations of the 137 sampled plots and to integrate these relation- 
ships with physical habitat data from the plots, using the ordination 
method of continuous classification developed by Bray and Curtis (1957). 

A two-dimensional ordination was constructed using the method of 
Beals (1960) (Fig. 6a) and an index of vegetational similarity (Bray and 
Curtis 1957): I=2W/atb, where a=sum of quantitative values of all species 
in stand A, b= the same for stand B, and W= sum of quantitative values 
the two stands have in common. Prominence value (Beals 1960) data for all 
vascular species (except trees) were used in calculating dissimilarity 
indices (I-l), where PV = mean covers x/ (frequency 3%). To test the validity 
of the ordination, 80 random pairs of stands were selected, and a correla- 
tion coefficient calculated between interstand ordination distances as 
measured with a ruler on the ordination field, and the indices of 


dissimilarity. The r value obtained was significant at the 13% level. 
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Figure 6: 


Ordination of the 137 plots showing 

the distribution of the vegetation types, 
several individual species populations 
and soil groups on the ordination field; 
(a) Plot numbers, (b) Vegetation types, 
(c) Dominant tree species, (d) Pyrola 
asarifolia, (e@) Rubus pubescens, 

(f) Ledum groenlandicum, (g) Hylocomium 
splendens, (h) Pleurozium schreberi, 

(i) Sov group. 
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When the 15 vegetation types are delineated on the ordination field 
(Fig. 6b), it is apparent that several of the types are not distinct from 
adjacent plots of another type. This situation arises as a result of the 
inclusion of species with high cover. Some types have similar understory 
vegetation but different tree dominants. Prominence values for the tree 
stratum were not calculated and were not included because of their tendency 
to determine the ordination due to their generally high cover values. The 
position of the plots on the ordination field is determined largely by the 
Species with the highest cover values; that is, these species contribute 
the most to the a and b values in the Similarity index calculation. The 
most abundant plant species, for example lodgepole pine, are often not the 
best indicators of the environment or physical habitat. As discussed 
previously, the vegetation types were distinguished according to species 
which are by inference, indicative of certain environmental or habitat 
conditions. 

The distribution of certain species within the 137 plots becomes 
apparent when their cover values are plotted on the ordination field. Tree 
cover did not enter into the similarity value calculations, but it became 
apparent that plots with aspen-or white spruce-dominated canopies occurred 
together on the ordination field. These groupings indicate an assemblage 
of understory species that separate the aspen and white spruce plots from 
lodgepole pine-dominated plots for example (Fig. 6c). The white spruce, 
black spruce and lodgepole pine groupings are discontinuous, indicating 
differences in understory vegetation. These differences are discussed in 
the previous section. 

Many other individual species show patterns on the ordination. Five 
of these are illustrated. Pyrola asarifolia (Fig. 6d) is most likely to 


occur in the white spruce, aspen, and lodgepole pine-green alder (Type 10) 
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stands. Rubus pubescens (Fig. 6e) shows a similar distribution. Ledum 
groenlandicium (Fig. 6f) is most abundant in the black spruce-dominated 
Stands and in the lodgepole pine-black spruce-Labrador tea-tall bilberry 
type. The two most common feathermosses also show patterns on the ordina- 
tion. Hylocomium splendens (Fig. 6g) is most abundant in the white spruce 
and in the imperfectly drained black spruce stands. Pleurozium schreberi 
(Fig. 6h) on the other hand, is most abundant in the lodgepole pine- 
dominated stands. 

When the various soil units are plotted on the ordination (Figiai6z) 
some order is apparent. Certain vegetation types and individual plant 
Species have a higher frequency of occurrence on some soils than on others. 
Relationships between vegetation types and soils are discussed in more 


detail in a previous section. 


G. DUNCAN'S MULTIPLE RANGE TESTS 


One of the original objectives of this study was to compare the site 
productivity of the various soil map units. Duncan's (1955) multiple range 
test is an appropriate means of testing for statistically significant 
differences amoung soil map units. However, when MAI for the 137 plots 
was stratified according to map unit, it became apparent that productivity 
varied widely within units and that there were very few plots sampled within 
some units. It is likely that much of the site variation within soil 
units is due to soil and topographic features that are important to tree 
growth, but that are not well described in soil unit definitions. For 
example, surface soil depth, subsoil texture, aspect, slope position and 
Slope gradient are factors often closely related to site quality but may 


vary widely within definitions for certain soil mapping units (Carmean 1977). 
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Thus, much of the variability in MAI within the Wapiti soil units is likely 
due to the exclusion of, or lack of definition of important site factors 
in, the map unit definitions. Duncan's tests on MAI among the soil units 
are therefore not included here. Differences in MAI for all plots, and in 
Site index for pine-dominated plots could be more readily demonstrated. 

The results of Duncan's (1955) multiple range test on means of mean 
annual increment in total volume (MAT) among the 15 vegetation types and 
On means Of lodgepole pine site index in the seven pine-dominated types 
are illustrated in Figs. 7 and 8. Duncan's tests on means of site index 
at 70 years were done only for the lodgepole pine types. 

The most productive type with respect to MAI was the Populus balsam- 
fera/Rosa acicularis/Thalictrum venulosum forest (Type 2). The aspen 
forests (Type 1) ranked third overall in terms of MAI with a mean of 
4.3 fenal yo (ig. Piles 

The most productive white spruce forest type in terms of MAI was the 
Picea glauca/Rubus pubescens-Maianthemum canadense type (Type 3) with a 
mean MAI of 4.7 Heine Saree 7). The least productive white spruce type 
in terms of MAI (3.6 ehaday ee was type 5, Picea glauca / Hylocomium 
splendens. 

Several of the lodgepole pine forest types were more productive than 
some of the white spruce types. The most productive pine type was type 8, 
Pinus contorta/Viburnum edule/Rubus pubescens, which occupies similar sites 


5 es 
to white spruce type 3 (Fig. 4). The average MAI in type 8 was 4.2 m ha 


yr (Fig. 7)z mean site index for pine was 19.8 m at 70 years (Fig. 8). 
The least productive lodgepole pine forest type was type 14, Pinus contorta/ 
Vaccinium myrtilloides/Cladonia spp., with an average MAI of only 1.3 mbar 
roe (Fig. 7) and a mean site index for pine of only 11.9 m at 70 years 


(Fig. 8). Type 15, Picea engelmannii-Abies lasiocarpa / Menziesia glabella 
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Figure 7: Duncan's multiple range test on means of mean annual total volume increment (m ha yr —) for 
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Figure 8: Duncan's multiple range test on means of site index (m) for the seven lodgepole pine- 
dominated vegetation types. Statistically similar means are underlined; significantly 


different means at any confidence level are not underlined, 


19 


had a higher MAI (2.5 mo hee yee) than type 14 but had a site index for 
Engelmann spruce of only 5.5 m at 70 years. 

The black spruce types had the lowest productivity. The Picea mariana- 
Ledum groenlandicum-Equisetum arvense type (Type 6) with an average MAI of 2.1 
Berar yr (Fig. 7) a mean site index of 8.2 m at 70 years was somewhat 
more productive than the Picea mariana/Ledum groenlandicum/Rubus chamaemorus 
bog forests (Type 7) which had an average MAI of only 0.4 Tea re ia (Bigs 7) 
and a mean site index of only 4.3 m at 70 years. 

The fact that several means, representing several vegetation types 
are underlined in any particular grouping in the Duncan's tests, indicates 
that the various vegetation types are not mutually inclusive in terms of 
productivity and that some of them are guite variable in terms of ranges of 
MAI or site index for the lodgepole pine-dominated plots. The Duncan's 
tests do show however, the relationships of the vegetation types to each 
other in terms of MAI, and site index for the pine types and should be 
useful in rapid evaluation of productivity classes for forest management 
purposes. 

Only site index data from plots of vegetation types dominated by lodge- 
pole pine were amenable to Duncan's test analysis because the small number 
of vegetation types dominated by the other tree species (three by white 
spruce, two by black spruce and only one each by aspen, balsam poplar and 
Engelmam spruce) precluded meaningful analyses using Duncan's test. 


H. STEPWISE MULTIPLE REGRESSION 


Reet he Forests of the Wapiti Map Area 


1.1 The Lodgepole Pine Forests 


The sampled lodgepole pine forests have a mean age of 96 years, a 


3, - - : 
mean annual total volume increment (MAI) of 3.0 m ha yr ', a mean site 
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andex of 16.2 m at 70 years and complete their maximum periodic annual total 
volume increment height at a mean age of 77 years. These forests occur 

at a mean elevation of 1133m, have a mean total vascular plant cover of 41% 

and a mean total moss cover of 45%, of which 5% is Bylocomium splendens ard 

32% is Pleurozium schreberi (Table 3). 

The mean soil profile is moderately well drained with a clay inA 
horizon/clay in B horizon ratio of 0.5 indicating Luvisols. The mean depth 
to mottling is 48 cm; mean thickness of moss and litter is 7 cm and mean 
thickness of A horizon (include Ael, Ae2 and Bm if present) is 19 cm. The 
average B horizon (usually Bt) has a thickness of 26 cm, medium subangular 
blocky structure with friable consistence, a hue of 5.8YR, value of 4.5 
and a chroma of 3.6. The average coarse fragment content of the profile is 


ye (lable 3).. 


1.2 The White Spruce Forests 

The sampled white spruce forests have a mean age of 123 years, MAI 
Of 3...43 mha "yoo" a mean site index of 14.0 m at 70 years and complete 
their maximum periodic annual total volume increment at a mean age of 104 
years. These forests occur at a mean elevation of 956 m, have a mean total 
vascular plant cover of 34% and a mean total moss cover of 41% (same as 
lodgepole pine forests), of which 22% is Hylocomium splendens and 11% is 
Pleurozium schreberi (Table 4). 

The mean soil profile is imperfectly drained with a clay in A horizon/ 
clay in B horizon ratio of 0.7 indicating Luvisols. The mean depth to 
mottling is 30.5 cm; mean thickness of moss and litter is 15.2 cm and mean 
thickness of A horizon (includes Ael, Ae2, and Bm if present) is 8.6 cm. 
The average B horizon (usually Bt) has a thickness of 18.8 cm, structure 


grading toward medium angular blocky with friable consistence, a hue of 
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Table 3 Means and standard deviations for lodgepole pine plot data (n=83). 
Dependent variables used in multiple regression analyses 


are indicated by Y; independent variables are indicated by X. 


Variable Mean 
3, -l ; 
Total volume (m ha ") S03. 935 
Plot age- Sa 95.83 
Bie BS =-1 
Mean Annual Increment (m ha yr ) 3.408 
: Sele at 
YIOMAL (rotation ages m sha syr ©) 3.009 
Y2 ,X64 Tree density (stems haa 2036 
Y¥3 Site index at 70yr Poms 
2 
Total volume per tree (m ) 0.347 


Y4 Periodic annual total volume increment per 


tree 0.005 
Maximum PAI per tree 0.007 
Y¥5 Age of maximum PAI ; 76.95 

2 
Total basal area (m ha ) 40.559 


¢ pee de 
Y6 Mean annual basal area increment (m ha yr = 0.423 


Merchantable volume ee 264.420 
Xl Elevation (my LPS 2535 
X2 log Elevation 3.05 
X3 Slope % 8.91 
x4 Aspect* i oelel 
X5 Aspect X Slcpe Emo 
X20 3 Sit+tC in A horizon 55.298 
X21 Thickness A horizon (cm) 18.64 
X22 Thickness A (Si+C A hor.) TO55ia05: 
X23 % Si+cC in B horizon 55.40 
X24 % Clay in B horizon 28.20 
X25 Thickness B horizon (cm) 25 dT 


X26 Thickness B (Si+C B hor.) 1445.15 


+ 


1+ 


Standard deviation 


103.349 


34.93 


0.098 


108.993 


Isai s! 


-0.46 
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X27 


X28 


X29 


X30 


X31 


X32 


X33 


X34 


X35 


X36 


Ey 


X38 


x60 


X61 


X62 


X63 


X39 


x40 


X80 


X81 


X82 


X83 


x84 


x85 


X86 


Table 3 cont. 
Variable 
Clay in A hor/Clay in B hor 
% Stone in profile 
of : I 
Internal profile drainage 


Depth to distinct mottling (cm) 


Hydraulic conductivity (cm day”) 


log Hydraulic conductivity 
: d 
Structure B horizon 
P . d 
Consistence B horizon 
: 2 
Chroma A horizon 
Hue B horizon 
Value B horizon 
Chroma B horizon 
Canopy cover 
Litter cover 
l/log Litter cover 
Deadfall cover 
Thickness organic horizon (cm) 


log Thickness organic 


Regeneration density (stems hao) 


Rosa acicularis cover 

Ledum groenlandicum cover 

log Ledum cover 

Rubus pedatus cover 
Epilobium angustifolium cover 


Cornus canadensis cover 


Ui 


Standard deviation 


0.282 


LS i615 


0.71 


38.99 


192.89 


0.94 


0.65 


ULE cha ae Cea 


Taber ss cont. 


Variable Mean Standard deviation 
X87 Rubus pubescens cover 0.38 z 0.86 
X88 Calamagrostis canadensis cover Oa a 0.41 
X89 Vascular plant cover 41.32 = 18.03 
X90 log Vascular plant cover LS : 0.19 
X91 Hylocomium splendens cover 51.03 = 6.70 
X92 Pleurozium schreberi cover 31.86 E 20575 
X93 log Pleurozium cover 1 al, 4S = 0.54 
X94 Polytrichum commune cover 0.41 a Cass 
X95 Total moss cover 45.01 z 23.42 
X96 log Total moss cover 5S Z 0.36 
X97 Total lichen cover 1.86 = 3.70 
Footnotes 
2 Merchantable volume calculated using minimum top diameter of 7.6cm 
b Xl to XS - site variables; X20 to x40 - soil variables; x60 to x64 - stand 


variables X80 to X97 - vegetation variables 

S Aspect represented by: sine (azimuth + 78.7)+1 

d Internal profile drainage, Structure B horizon and Consistence B horizon 
scales are described in Appendix 2 


e Soil colours are Munsell designation 
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Table 4 Means and standard deviations for white spruce plot data (n=30). 
Dependent variables used in multiple regression analyses 


are indicated by Y¥; independent variables are indicated by X. 


Variable Mean Standard deviation 
Tdtal volome, (mre) 462.194 99.317 
Plot age 122.83 = 35.54 
Be t=le =i + 
Mean Annual Increment (m ha yr ) 4.016 - 1.280 
- - - 
Yl MAI (rotation age; ane we He 3.445 - PASC 
- S 
Y2 Tree density (stems ha * 1507 - 618 
+ 
Y3 Site index at 70 yr 13.96 - 4a; 
3 + 
tal volume per tree (m ) 0.538 = Q2232 
Y4 Periodic annual total volume increment per 2 
tree 0.009 - 0.005 
+ 
Maximum PAI per tree 0.011 - 0.005 
+ 
YS Age of maximum PAI 103.88 - S05 
PRS ie + 
Total basal area (m ha ) 50.803 - 8.908 
F 201s + 
Y6 Mean annual basal area increment (m ha yr ) 0.413 - 0.073 
. 4 
3, -l,a + 
Merchantable volume (m ha ) AUS SON - 94.467 
: b + = 
Xl Elevation (m) 956.06 - 245.74 
+ 
X2 log Elevation 2.98 = -0.41 
+ 
X3 Slope 3% 6.62 = eA. 
(ol + 
X4 Aspect 1.36 = Oss 
X5 Aspect X Slope £6.35 = tees wis) 
_ 
X20 %$ Si+C in A horizon 56.42 - Zoo 
+ 
X21 Thickness A horizon (cm) 8.64 ee 9.37 
aa 
X22 Thickness A (Si+C A hor) 464.49 = 512.02 
X23 % SitC in B horizon 63.20 = 21.70 
+ s 
X24 % Clay in B horizon BS Onc = 16.73 
+ 
X25 Thickness B horizon (cm) 16.87. - LOR 


X26 Thickness B (Si+C hor) 1105.74 - 812.24 
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Table 4 cont. 
Variable 
Clay in A hor/Clay in B hor 
% Stone in profile 
Internal profile Geadaces 


Depth to distinct mottling (cm) 


Hydraulic conductivity (cm day +) 


log Hydraulic conductivity 
: d 
Structure B horizon 
: : d 
Consistence B horizon 
f e 
Chroma A horizon 
Hue B horizon 
Value B horizon 
Chroma B horizon 
Canopy cover 
Litter cover 
lflog Litter cover 
Deadfall cover 
Thickness organic horizon (cm) 


log Thickness organic 


: : ==! 
Regeneration density (stems ha ) 


Rosa acicularis cover 

Ledum groenlandicum cover 

log Ledum cover 

Rubus pedatus cover 
Epilobium angustifolium cover 
Cornus canadensis cover 


Rubus pubescens cover 


Standard deviation 


0.347 


17.49 


1.08 


35.81 
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Table 4, cont. 


86 


Variable Mean Standard deviation 

X88 Calamagrostis canadensis cover 0.09 Z 0.20 
X89 Vascular plant cover Sate 2 22.44 
X90 log Vascular plant cover eae = Q.35 
X91 Hylocomium splendens cover PEI Es 18.80 
X92 Pleurozium schreberi cover 10.60 z 12500 
X93 log Pleurozium cover 0.66 S Olediz 
X94 Polytrichum commune cover 0.61 2 2.84 
X95 Total moss cover 40.75 . 24.01 
X96 log Total moss cover 1.48 0.42 
X97 Total lichen cover 0.42 . 0.42 
Footnotes 
a Merchantable volume calculated using minimum top diameter of 7.6 cm 
5 Xl to X5 - site variables; X20 to x40 - soil variables; X60, to. X63) = 

stand variables x80 to X97 - vegetation variables 
c Aspect represented by: sine (azimuth + 78.7)+l 
d Internal profile drainage, Structure B horizon and Consistence B horizon 


scales are described in Appendix 2 
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5.4YR, value of 4.2 and a chroma of 2.5. The average coarse fragment 


content of the profile is 6% (Table 4). 


1.3 The Black Spruce Forests 

The sampled black spruce forests have a mean age of 138 years, MAI of 
3 ay a mean site index of 11.3 m at 70 years and complete their 
maximum periodic annual total volume increment at a mean age of 128 years. 
These forests occur at a mean elevation of 1044 m, have a mean total vascular 
plant cover of 32% and a mean total moss cover of 63%, (Ole which’ 21%. as 
Hylocomium splendens and 19% is Pleurozium schreberi (Table? 5). 

The mean soil profile is imperfectly drained with a clay in A horizon/ 
clay in B horizon ratio of 0.8 indicating Luvisols but less clay accumula- 
tion in the B than in the white spruce and lodgepole pine forests. The 
mean depth to mottling is 23.9 cm; mean thickness of moss and litter is 
36.1 cm and mean thickness of A horizon (includes Ael, Ae2 and Bm if present) 
235.9 Cm. The average B horizon (usually Bt) also has a thickness of 8.9 
cm, is structureless to medium angular blocky with friable consistence, a 
hue of 4.8YR, value of 4.1 and a chroma of 2.6. The average coarse fragment 


content of the profile is 3% (Table 5). 


1.4 The Aspen Forests 

The sampled aspen forests have a mean age of 79 years, a mean 
Productivity of 4.48 bee a mean site index of 19.2 m at 70 years 
and complete their maximum periodic annual total volume increment at a mean 
age of 56 years. These forests occur at a mean elevation of 777 m, have a 
mean total vascular plant cover of 40% and a mean total moss cover of 10%, 
of which 5% is Hylocomium splendens and 3% is Pleurozium schreberi (Table 6). 


The mean soil profile is moderately well drained with a clay in A horizon/ 


* 
7 


ae eae 


st > ert 


i3 & 


3 bw 


tees e82  -8+%. 30 suo mis a be t.s 36 


knsot acem c han 80) So zOVe5 Jnalq 1 


% eplecuge!s go: €£ Bow earreinio 4k. 42, ff  moened wa E 


“ 


= 


as 
oy 76  £' So wetted 6928 faee a 


arew tiew taetecaban #i sl iSong 2 ; 


ml 7) 


- 
ip olde?) 20 ei silttome ¢ 


j 
7 


stea2ot eovIde a 


a # Sve sented? ssuvge soetd 


i ngicavelse geen 6 35 38% 
sso ese fetot hesm # Das ane 
—owets ef 7el tos enshaa 
>} «f{iios@ ife 
efusll’ eniten tard 6.0 Qo Cheers 
Fyry i> aenlvee adie ae aw nang . 


Soom es wited OS 220k ok ptt 


<5 3exl A 2e eanitioms 

® 
on yihesey) castro? 3 egareve 

ew : ay) slmucé cei sped oF &e6 sleg 


7 
a«' mele 4s fap 1.) Io sola 


ma i Se sf fiewg 


etaer9@ 


SOS NGI0 wus aiestod neges 


” 


f~ | 

| wie aoew oo. Oy ed we Ob. Be 

iV 7 

ietes Launnmk nfSel reg mom 4m shay 
' ] 


® neo 4 34 SHOGe etestc3 oned? 


i 


Table 5 Means and standard deviations for black 
Dependent variables used in multiple reg 
are indicated by Y¥; independent variable 

Variable Mean 
7 cape! 
Total volume (m ha ) 253.494 
Plot age 138.25 
3. 5—]) =} 
Mean Annual Increment (m ha ‘yr ~) 2.003 
P , 3, .-1 -1 
Yl MAI (rotation age; m ha yr ) g Me 9 Bes 
oe ~ =k 
Y2 Tree density (stems ha ) 2858 
Y¥3 Site index at 70 yr 8.25 
3 ‘ 
Total volume per tree (m ) 0.191 
¥4 Pericdic annual total volume increment per 
tree 0.002 
Maximum PAI per tree 0.093 
Y5 Age of maximum PAI 128.16 
A -1 
Total basal area (m ha )} 39232 


. : oe Rowe. 1 an 
Y6 Mean annual basal area increment {m ha yr ) 0.284 


=p 
Merchantable volume fae os 185.825 
= b 

Xl Elevation (m) 1043.69 
X2 log Elevation 3.02 
X3 Slope ¢% 3-50 

c ee 
X4 Aspect i ey 
X5 Aspect X Slope 4.56 
X20 % Si+C in A horizon 68.10 
X21 Thickness A horizon (cm) 8.37 
X22 Thickness A (Si+C A hor) 556.77 
X23 % Si+C in B horizon 58.80 
X24 % Clay in B horizon 30.43 
X25 Thickness B horizon (cm) 8.97 


X26 Thickness B (Si+C B hor) 625.85 
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X27 


X28 


X29 


X30 


X31 


X32 


X33 


X34 


X35 


X36 


X37 


X38 


X60 


x61 


X62 


X63 


X39 


X40 


X80 


X81 


X82 


X83 


X84 


X85 


X86 


X87 


Table 5 cont. 
Variable 
Clay in A hor/Clay in B hor 
$ Stone in profile 


Internal profile Boa races 


Depth to distinct mottling (cm) 


Hydraulic conductivity (cm gaye) 


log Hydraulic conductivity 
: d 
Structure B horizon 
: ; d 
Consistence B horizon 
E e 
Chroma A horizon 
Hue B horizon 
Value B horizon 
Chroma B horizon 
Canopy cover 
Litter cover 
l/log Litter cover 


Deadfall cover 


Thickness organic horizon (cm) 


log Thickness organic 


F : -1 
Regeneration density (stems ha ) 


Rosa acicularis cover 
Ledum groenlandicum cover 
log Ledum cover 


Rubus pedatus cover 


Epilobium angustifolium cover 


Cornus canadensis cover 


Rubus pubescens cover 


Standard deviation 
0.269 


8.02 
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Table 5) cont. 


Variable Mean Standard deviation 
X88 Calamagrostis canadensis cover On5 Z One 
X89 Vascular plant cover eal es ph fl 
X90 log Vascular piant cover 1.46 ie 0.23 
X91 Hylocomium splendens cover 20.93 2 lish (a)? 
X92 Pleurozium schreberi cover 18.64 = PAO) A 
X93 log Pleurozium cover 0.92 = 0.67 
X94 Polytrichum commune cover O13 . 0.29 
X95 Total moss cover o3ng 2 ss 14.86 
X96 log Total moss cover is IS) = Q.11 
X97 Total lichen cover 0.82 = Pe 6 
Footnotes 
a Merchantable volume calculated uSing minimum top diameter of 7.6 cm 
b Xl to X5 - site variables; X20 to x40 - soil variables; X50 to 63- 


“stand variables x80 to X97 - vegetation variables 

c Aspect represented by: sine (azimuth + 78.7)+1 

d Internal profile drainage, Structure B horizon and Consistence B horizon 
scales are described in Appendix 2 


e Soil colours are Munsell designation 
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Table 6 Means and standard deviations for aspen plot data (n=17). 
Dependent variables used in multiple regression analyses 


are indicated by Y; independent variables are indicated by X. 


Variable 
Total volume jane 
Plot age 
32, vale =) 
Mean Annual Increment (m ha yr  ) 
ab 


Yl MAI (rotation age; m hal iy ) 
Y2 Tree density (stems pen 

Y3 Site index at 70 yr 

Total volume per tree fm) 


Y4 Periodic annual total volume increment 
per tree 


Maximum PAI per tree 
Y5 Age of maximum PAI 


Total basal area (ene 


: ah = 
Y6 Mean annual basal area increment (m ha yr 


; ; ey Sate 
Merchantable volume (m ha ) 


Xl Elevation ie 

X2 log Elevation 

X3 Slope % 

x4 Aspect® 

X5 Aspect X Slope 

X20 % Si+C in A horizon 

X21 Thickness A horizon (cm) 
X22 Thickness A (Si+C A hor) 
X23 % Si+C in B horizon 


horizon 


Ww 


X24 % Clay in 


X25 Thickness B horizon (cm) 


Ww 


X26 Thickness (Si+C B hor) 


Mean 
362.552 
78.85 

4.605 


4.511 


56.53 


42.608 


1) 0.540 


308.075 


776.34 


aie Oe 


24.35 


1708.22 


Standard deviation 


144.496 


23.16 


1.205 


1.433 


0.005 


0.005 


31.45 


WA xi) 


O5157 


136.829 


88.45 


-0.44 
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Table 6 cont. 


Variable Mean Standard deviation 
X27 Clay in A hor/Clay in B hor 0.406 z 0.295 
X28 % Stone in profile 0.62 z 1.08 
: : d - 
X29 Internal profile drainage PLS Tal - Os75 
X30 Depth to distinct mottling (cm) 29.26 ox S\5) 3p 
<3) Eydrsulie conductivity (cm day) 49.78 f 132.06 
; Wess - 
X32 log Hydraulic conductivity -0.06 - al eds} 
; d fe 
X33 Structure B horizon 4.29 = 1.36 
: : d + 
X34 Consistence B horizon sty = 0.64 
: e + 
X35 Chroma A horizon 2.24 - 0.83 
+ 
X36 Hue B horizon 5.74 - Bie 
+ 
X37 Value B horizon A ABIS = 0.60 
+ 
X38 Chroma B horizon Jy SE) = Lea 
he 
X60 Canopy cover 53-95 - 8.93 
2 + 
X61 Litter cover TAAL = Ly s08 
+ 
X62 l/log Litter cover -1.33 - ON Ly, 
+ 
X63 Deadfall cover S222 - SHALOM 
+ 
X39 Thickness organic horizon (cm) 7.69 - 6.47 
+ 
X40 log Thickness organic 0.89 - 0.54 
: : -1 + 
X80 Regeneration density (stems ha ) 319 - 361 
+ 
X81 Rosa acicularis cover 2.94 - 3102 
+ 
X82 Ledum groenlandicum cover 0.09 - Ov27 
+ 
X83 log Ledum cover -0.01 ~ 0.03 
+ 
X84 Rubus pedatus cover 0.00 - 0.00 
+ 
X85 Epilobium angustifolium cover 3.08 = Beno 
+ 
X86 Cornus canadensis cover 5.56 - 3.98 


X87 Rubus pubescens cover Qa - aly gees) 
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Table 6 cont. 


Variable Mean Standard deviation 

X88 Calamagrostis canadensis cover 0.42 as Oe 
X89 Vascular plant cover 40.29 : 22.97 
X90 log Vascular plant cover 1G SS) af On21 
X91 Hylocomium splendens cover sys akit 4 967 
X92 Pleurozium schreberi cover OZ a4 4.02 
X93 log Pleurozium cover OnS5 g 0.40 
X94 Polytrichum commune cover 0.10 M 0.41 
X95 Total moss cover Ses) 2 12.74 
X96 log Total moss cover 0.57 . 0.68 
X97 Total lichen cover 0.12 ‘ 0.17 
Footnotes 
a Merchantable volume calculated using minimum top a taner ae OL oem 
b XL to XS - site variables; X20 to x40 - soil variables; X60 to X63 - stand 

variables; X80 to X97 - vegetation variables 
S Aspect represented by: sine (azimuth + 78.7)+1 
d Internal profile drainage, Structure B horizon and Consistence B horizon 


scales are described in Appendix 2 
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elay~in B: horizon-ratio of 0.4 indicating Luvisols. The mean depth to 
mottling is 30.5 cm; mean thickness of moss and litter is 7.6 cm and mean 
thickness of A horizon (includes Ael, Ae2 and Bm if present) is 11.7 cm. 
The average B horizon (usually Bt) has a thickness of 24.4 cm, medium 
subangular blocky structure with firm consistence, a hue of 6.7YR, value 
of 4.4 and a chroma of 3.6. The average coarse fragment content of the 
Profile is 1%. (Table 6). 

The semi-quantitative scales for the soil attributes entered in the 


regression analyses are presented in Appendix 4. 


2. Transformation of Independent Variables 


The assumption behind the use of the stepwise multiple regression is 
that the relationship between the dependent variables and the independent 
variables is linear. To test whether or not this was in fact the case, 
each of the dependent variables was plotted against each of the independent 
variables using a computer scattergram program. A visually estimated 
best-fit line was drawn through the points on each of the resulting scatter- 
grams. In most cases a straight line seemed to well approximate the 
relationship between the variable pairs. In other cases a nonlinear 
relationship was evident. In the latter cases the log of the values 
(in the case of litter cover, the reciprocal of the log) of the independent 
variables was plotted against the dependent variables. The log transform- 
ation was regarded as best approximating a linear relationship where 
originally a nonlinear relationship existed. The effectiveness of the log 
transformation was evaluated by comparing the simple correlation coefficients 
between the dependent and independent variables before and after transform- 


ation. The simple correlation between the dependent variable and the 
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transformed independent variable should be greater than that between the 
dependent and non-transformed independent variable. The log-transformed 
variables were subsequently entered into the regression as separate independent 
variables (Tables 3 to 6). The untransformed independent variables were 
allowed to remain in the regression because some dependent variables showed 

a more linear relationship with the untransformed independent variable. 

The stepwise multiple regression program permitted the selection of either 

the transformed or untransformed independent variable depending on which 
exhibited the better linear relationship with the dependent variable. 

In two other cases a different transformation was applied. It seemed 
likely that an interaction between two related independent variables (ie. 
product of the two) might contribute more to the R? value in the regression 
than might the two of them alone. The independent variables involved were 
Si+C content of the A and B horizons and the thicknesses of those horizons. 
If the interaction variable entered the regression before one or the other 
components, that is to say, contributed more to the R? value, then the 
two components (thickness and Si+C content) were deleted. Conversely, if 
Si+C content of the A or B horizon or thickness of the A or B horizon 
entered the regression before the interaction variable, then the interaction 
variable was deleted. It became apparent that the interaction of Si+C 
content and horizon thickness did not contribute more to the R? value in 


the black spruce and aspen equations (Tables 12 and 13), than did Si+c 


content alone. 


3. Selection of the Dependent Variables 


For many years, the usual dependent variable used in forest-site 


studies has been site index, or height at a certain index age, usually the 


ae Aaeessed fic, nem: 


ae sbsehianiextetiot ail 


ai afi of ce ne 
‘ a ne lain xe ‘pAsoam Te ag c * 
‘TS OMeiacts dak ati ease GS OF pare 7 ! 


dy yaene” peat Chapt ay natn 
ee A sash wince 


Stk eens coe b ee ae ba vial once relia " , 
; ae ks 


ee aa cided baie vue! Adie 4 


of 
ae wl - 


as eae | 


yt ted ex’ oe Se a RAO NT apianeayee | ab 


Potaw CMTE fata 


sh adiaeiey sdotriaaibact bometanering 20:8 
tow Seunbmged vow ai MA heanatoeten redo 


i. eB i 


Hameae 23 


gu. het. Rm Ve gyi, ed serbsi sco 


‘a oe 
i» SS) eet A ae 3, et 


oie sedi oth) xcae aud? ates rapier font. “At 
laa aneteeegeberd itt orn lope oa senay 


eyew beebowns Fmdtes 


VAY dO apace rods ond toni enone a ites, 8 da Dot 


ate dt alec et iucabeateys' ate bes ete aicakym's 
CSS es aw Lee oa ‘eth ne 2a bosuele rene Kee Perer 
ae ae pcaabat why gee arya: Cp fain: or bit AY ed Ov a me 


or 
were) ed. wea ‘gat hess iguiones es a) is wo “fh: Ghd Bo Jaws 

jooratned off: abdd yskGat ray) onto Yoda t: meets ‘excihond aotasoseer ena 4 

be ed neettomeid Gait) gaa ral stint adh | bedenst anh 


nk Suter, *reioee sont stpik easton son beh bonatoiss sos inet Grit 2 


re 


TS, 


im bb ed ae Nene “ ey lao) wnolawipe qogme Bris. pore a 


sAie-teexdh ol bake. whee alla ld Javed ott) yaeesy oghett see = 
edt leven oe sess pie trae Se einen, 0 cael NA rs 


rotation age for that species (Alban 1972, 1974; Carmean L965), 42.967:,, 19704 
1971; Graney and Ferguson 1971, 1972; Holmes and Tackle 1962). In other 
studies (Duffy and England 1967, Duffy and Knight 1967, Kabzems et al. 1972), 
mean annual total volume increment (MAI) has been used as an indication of 
forest site productivity. This measure is the basis of the Canada Land 
Inventory land capability rating for forestry. Periodic annual total 

volume increment (PAI), or the volume increment per unit area at a specified 
period of years divided by the number of years, has sometimes been used as 

a measure of forest productivity (e.g. Dumanski et al. 1973). In this thesis 
several dependent variables including MAI (rotation age), density, site 
index, PAI/tree, age of maximum PAI and mean annual basal area increment 

are used. 

Periodic annual total volume increment per tree (PAI) will be discussed 
first. It represents the periodic annual total volume increment for the ten 
years prior to rotation age (70 years for lodgepole pine and aspen; 100 
years for the three spruce species) as determined from the stem analyses. 

This value could not be multiplied by density to find PAI per unit area 
because bias in selection of the stem analysis sample trees from which PAI 
was estimated. This bias was towards trees of superior form and apparent 
lack of mechanical or disease damage such as might be expected of trees in 
a well-managed forest stand. More elaborate corrections to estimate PAI 
On an areal basis were not explored due to the opportunity to determine a 
good estimate of MAI, generally regarded as a superior expression of site 
productivity, from the stem analysis data. 

Mean annual increment is basically the total tree volume of a given 
area of land divided by the age of the forest in years, MAI estimates are 
most useful when the forest stand in question is at or near the commercially- 


accepted rotation age for that species in the area in question (i.e. when 
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both PAI and MAI has peaked). Tree growth (volume or height) exhibits the 
sigmoid growth curve characteristic of many Organisms and populations, 

with slow initial growth, rapid intermediate growth and finally much 

reduced growth as the organism approaches maturity. This trend is illustrated 
in the stem analysis computer plot of height and volume growth of a lodgepole 
pine tree in plot 12 (Fig. 3). Thus, mean annual increment estimates tend 

to underestimate potential productivity in old mature stands because volume 
increment decreases after a certain age. Conversely, if a MAI estimate is 
made in a stand much younger than the rotation age, MAI may Be overestimated, 
assuming that the population growth will normally have tapered off by the 
time it reaches the rotation age. It is especially desirable that prior 

to harvesting pulpwood species such as lodgepole pine, in order to maximize 
monetary return, the age of maximum PAI be at least attained, although it 

is generally regarded as desirable to wait until the MAI growth curve also 
has peaked, with the rotation age taken as the age when the PAI and MAI 
curves cross (Davis 1954 p. 228). For species usually used for lumber, 

such as white spruce, it is often desirable to let the stand grow for a 
longer time to improve the quality of the saw timber harvested. It is 
apparent from Davis' (1954) discussion, that rotation age can not be 
determined on the basis of gross productivity alone. Different utilization 


Standards and sites will cause some deviation of rotation age from average 


(Kirby 1962). 


3.1 Volume Ratio Calculations 


A method to standardize MAI estimates to a base age (70 years for 
lodgepole pine and aspen and 100 years for the three spruce species) is used 
here. These base ages approximate those currently viewed as optimum 


rotation ages in the western Alberta forest industry, and not those ages 
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suggested to be optimum from data of this study. The volume ratio corrects 
MAI's forward in stands less than the base age and backwards in stands 
greater than the base age and is specific to the different tree species 
sampled for the area concerned. For each plot, the mean volume at rotation 
age as determined from the stem analysis, was divided by the mean volume 

of the sample trees at time of sampling to obtain the "volume ratio." 

These volume ratios were plotted against the respective age for all field 
plots, with separate graphs for each of the tree species concerned. A 


curvilinear relationship was evident. Volume ratio plotted against log 


age produced a nearly linear relationship. Subsequently a linear regression 


of volume ratio on log age was calculated for each tree species (Fig. 4.9) 
Next the volume ratio at rotation age was calculated for each tree species 
using the regression equations. This value was then multiplied by the 
gross MAI to give MAI at rotation age. It is this value that is used as 

a dependent variable in the stepwise multiple regression equations 
calculated. Itis likely that this method probably gives somewhat conserva- 
tive estimates of MAI for the stands much older than rotation age because 


it does not take into account mortality which has occurred between the 
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rotation age and the age at time of sampling. However, it is believed that in 


most cases, the decrease in the MAI estimate due to decreased productivity 
after annual increment has peaked, is greater than loss due to mortality 
during that period of time. The volume ratio regressions are illustrated 
ahoFig. 9. 

From the volume ratio graphs it became apparent for the white spruce 
plots that there were two very distinct white spruce populations in terms 
of volume growth (Fig. 9b). A closer examination of the sample plot data 
revealed that the population represented by the lower regression line 


(Fig. 9b) was composed of sample trees which occupied understory position 
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Figure 9: Volume ratio equations for the four 


tree species sampled: (a) lodgepole 
pine, (b) white spruce, (c) Black 
spruce, (d) aspen. 
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Figure 9a: lodgepole pine 
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Figure 9b: White spruce 
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Figure 9c: Black spruce 
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Figure 9d: Aspen 
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in either aspen or lodgepole pine forests and were thus somewhat suppressed. 

Tree density is not usually used as a dependent variable as it is in 
this thesis. It seemed desirable to determine if factors other than stand 
history and age were important in determining stand density. 

Site index is used in this study as it has been in many others. It 
is still the most commonly used expression of forest productivity. 

The age of maximum periodic annual total volume increment (PAI) is 
used here as a dependent variable, and is proposed as a means of estimating 
optimum rotation age for a particular site independently of direct measure- 
ments of the trees, enabling this determination to be made on cleared or 
burned sites. 

Mean annual basal area increment is another dependent variable used 
in this study but is not commonly used elsewhere. It is expected that it 
should show similar relationships to the independent variables as shown by 


MAT e 


4. The Multiple Regression Equations 

The equations which are listed in Tables 7 to 14 represent data from 
83 lodgepole pine stands, 30 white spruce stands, 15 black spruce stands 
and 17 aspen stands. For each of the 36 equations the independent variables 
are listed in order of importance, relative to the amount they contribute 
to the R? value. A brief description taken from Tables 3 to 6 listing 
the variables that entered the regression, their means, variances and 
standard deviations for each of the sampled stand groups, follows. Component 
vegetation types within each of the stand groups are described in more 


detail in the previous sections. 
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With one exception, the black Spruce equations accounted for the largest 
amount of variability in the six dependent variables: MAI, tree density, 
Site index, PAI per tree, age of maximum PAI, and mean annual basal area 
increment, with lesser amounts of variability accounted for by the aspen, 
white spruce and lodgepole pine equations (Tables 7 to 14). The only 
exception was that the tree density prediction equation for aspen (R7=0.99) 
was Slightly higher than that for black spruce (R7=0.98). The prediction 
equations were on the whole very satisfactory with a large proportion of 
the variation in the dependent variables being accounted for by the equations. 
For black spruce and aspen the prediction equations were nearly equally 
effective with Ro values ranging from 0.94 to 0.99 (Tables 9 and 10). For 
the white spruce stands, the strongest prediction equation was for mean 
annual basal area increment (R7=0.94) and the weakest was for stem density 
(R7=0.70). For the lodgepole pine stands the strongest prediction equation 
was for PAI per tree (R7=0.76) and the weakest was for mean annual basal 
area increment (R7=0.44). Site index and MAI are probably the most widely 
used expressions of site productivity. For all four species, site index 
was predicted with slightly greater precision than was MAI. For all the 
equations presented, the R? values are significant at the 1% probability 
level. 

To determine the effect of the physical site factors alone upon the 
selected dependent variables,MAI (rotation age), site index and age of 
maximum PAI, the vegetation cover values and vegetation-related variables 
were deleted from the regression. Equations derived in this manner should 
allow estimation of the dependent variables upon sites where trees, under- 


story and litter have been obliterated by major disturbances such as 
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logging or fire, and should provide estimates for soil survey sites where 
plot data can not be obtained. Deletion of the vegetation-related indepen- 
dent variables resulted in a decrease in the R? and F values for the 
regression equations. These decreases were greatest in the lodgepole pine 
equations, though all F ratios were still significant at the 1% level (Table 
11). The white spruce equations showed the next greatest decrease in R? 
and F values but F ratios were still significant at the 1% level except for 
age of maximum PAI (significant at 5% level; Table 12). The aspen equations 
showed a smaller reduction in the R? value when vegetation variables were 
deleted than did those for lodgepole pine and white spruce, but showed very 
substantial decreases in the F ratio values for MAI and site index equations 
(Table 13), which were not significant at even the 5% level. The black 
spruce equations showed a tendency similar to those for aspen with relatively 
small decreases in R? but large decreases in F ratio. The F ratio for the 
black spruce MAI equation was significant only at the 10% level while that 
for site index was significant at the 5% level. The black spruce age of 
maximum PAI equation remained significant at the 13 level (Table 14). 
Hagglund and Lundmark (1977) discovered that site index equations 
calculated for Scots pine and Norway spruce in Sweden by vegetation type 
were more reliable than equations calculated for the whole tree species 
group as was done in the present study. The small number of sampled stands 
within several of the vegetation types prohibited this approach for the 
Wapiti map area. Lowry (1976) reported that when both individual species 
and species groups were used in multiple regression equations for black 
spruce in eastern Canada, the species groups showed higher predictive 
precision than did individual species. Species groups were not used in 
the present study, and it seems unlikely this approach would contribute 
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substantially in most cases to the already high R’ values. Hofman (1976) 
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also recently discovered the utility of vegetation variables in estimating 
beech productivity in East Germany. Twenty site factors including soil 
physical and chemical properties plus the ground cover percent of nine 
indicator species were used. 

The significance of the individual independent variables in the 
regression equations was determined by t-tests. The level of Significance 
(5 or 1% level) is indicated for each independent variable in the equations 
of Tables 7 to 14. Some variables with non-significant t-tests appear in 
the regressions. These variables however still contribute at least 1% to 
the R? value up to a maximum of nine independent variables per equation 
(the criterion for order of occurrence in regression equation being the 
contribution to Bae Discussion on the ecological significance of the 
independent variables is limited to those established to be statistically 
Significant as determined by the t-tests in addition to being significantly 


correlated with the dependent variable (r value significant to at least p=.10). 


5. Independent Testing of Regression Equations 


In order to determine the utility of the regression equations presented 
in Tables 6 to 9 it seemed desirable to attempt to test their validity. 
Multiple regression equations, even though they have a high R? value may 
not necessarily be reliable when applied to data independent of that from 
which the equations were derived, due to inadequate sampling, incomplete 
data analysis or chance correlation. Errors due to the above factors will 
not be detected until the equation is tested with independent data 
(McQuilkin 1976). This test was accomplished by randomly selecting 70% of 
the plots for the tree species in question and developing growth predic- 


tion equations by stepwise multiple regression techniques, from the data 
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from those plots, using MAI and site index as dependent variables. The 
equations developed in this manner were used to predict MAI and site index 

in the remaining 30% of the plots. A simple correlation coefficient between 
actual and predicted values for MAI and site index was subsequently 
calculated. This procedure was repeated four times, with vegetation variables 
included and deleted for both lodgepole pine and white spruce plots. The 
equations developed from the random selection were terminated after nine 

steps like the originals and after as many steps as were necessary to 

include all the independent variables present in the Original equations. 

The best of the four equations, ie. those which corresponded most closely 
to the equations developed from all plots for the particular tree species 
and dependent variables concerned are presented in Tables 15 and 16. 

In general, there appeared to be good correspondence between the 
Original equations and those derived from the random selection, though the 
agreement was better for the lodgepole pine equations than for those develop- 
ed for the white spruce plots. The best equation developed from the random 
selection for lodgepole pine MAI (vegetation variables included) contained 
the same nine variables as the original equation in nine steps though the 
variables were arranged in slightly different order (Tables 7 vs. 15). 

The MAI (vegetation variables deleted) equation calculated from the random 
plot selection for pine contained five of the original six variables in 

the first six steps (Tables 11 vs. 15). Twelve steps were required to 
include all six of the original variables. The results for the pine site 
index equations developed from the random sample were similar, with eight 
of the original nine variables being present in the equation including 
vegetation variables and all seven of the variables of the original equation 
with vegetation variables deleted (Tables 7 vs. 15). With the exception of 


MAI (vegetation variables deleted), the F-ratios calculated for the lodgepole 
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pine equations derived from the random plot selection and all the simple 


correlation ccefficients between observed and predicted values for MAI 


and site index were significant to at least p=.05. 


The best equations developed from the random selection of white spruce 
plots showed less resemblance to the original equations than did the pine 
equations. The equations for white spruce MAI developed from the random 
plot selection contained four out of nine and four out of seven of the 
variables present in the original equations with vegetation variables 
included and deleted, respectively (Tables 8 vs. 16 and 12 vs. L6)i She 
white spruce site index equations contained six out of nine and five out 
of seven of the variables present in the original equations with vegetation 
variables included and deleted, respectively (Tables 8 vs. 16 and 12 vs. LG). 
Only the equations for white spruce MAI and site index (both with vegetation 
variables included) had F-ratios significant at p=.05. None of the Simple 
correlation coefficients calculated between observed and predicted values 
for the four white spruce equations in Table 16 were significant at p=.05. 

In summary, the lodgepole pine regression equations are in general, 
reliable and may be used with caution within the Wapiti map area. The 
white spruce equations are less reliable. A larger and broader data base 
encompassing a wide range of soil and site characteristics is desirable for 
white spruce, black spruce and aspen forests within the Wapiti map area, 


if the equations for these species are to be widely used. 
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I. DEPENDENT - INDEPENDENT VARIABLE RELATIONSHIPS 


The following discussion is an attempt to discern ecological rela- 
tionships among the various dependent variables and independent variables 
that had statistically significant t-tests (to at least p = .05) within 
the growth prediction equations (Tables 17 to 20). Attempts to explain 
relationships between the dependent and independent variables in question 
were made only where significant simple correlations to at least the 
p = .10 level occurred between the dependent and independent variables, 
as indicated by asterisks on the simple correlation coefficients in 
Tables 17 to 20. Thus in the following discussion when references are 
made to "Significant" independent variables, these variables had 
statistically significant t-tests (to at least p = .05), and when 
references are made to "correlations" between dependent and independent 
variables, it can be assumed that the simple correlation coefficients 
in question are statistically significant to at least the p = .10 level. 

Variables discussed represent physical habitat attributes and 
vegetation attributes (species cover). The physical site attributes are 
those of the external environment of the site such as slope gradient, 
Slope aspect and elevation and also readily field-determinable physical 
properties of the soil, especially as they may influence soil moisture 


regime and tree rooting characteristics. 


1. External Site Factors 


External site factors include elevation, slope gradient, slope 


aspect and slope gradient X aspect interaction. 
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Table 17 


Elevation (m) 


log Elevation 
Slope percent 


Slope aspect 
Ranect X slope 

% Si+C in A horizon 
Thickness A horizon (cm) 
Thickness A 


% Si+C in B horizon 


% Clay in B horizon 


Thickness B horizon (cm) 


Thickness B 
Clay in A 
% Stone in profile 


Internal profile drainage 


Depth to distinct 
mottling (cm) 


Hydraulic conductivity 
(cm day °) 


log Hydraulic conductivity 
Structure of B horizon 
Consistence B horizon 
Chroma A horizon 

Hue B horizon 


Value B horizon 


log Total moss cover 


Total lichen cover 


for the lodgepole pine plots (n=83) 


Mean annual 
total volume 


(Si+C in A hor) 


(Si+C in B hor) 


hor/Clay in B hor 


wae 


xs 
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Simple correlation coefficients between dependent 
variables and significant independent variables as 
determined by t-tests in multiple regression equations 


PAT 
per 


Age of 
maximum 


Mean annual 


Stem basal area 


Site 


increment density index tree PAT increment 
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Significant at p = .05 
Significant at p = .10 


Independent variable occurring only in equations with 
vegetation variables deleted 
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Table 17 cont. 


Mean annual 
total volume 


Stem 


Site 


PAI 
per 


Age of 
maximum 


increment density index tree PAL 


Chroma B horizon 

Canopy cover + 
Litter cover 

l/log Litter cover - 
Deadfall cover = 
Stem density 


Thickness organic 
horizon (cm) 


log Thickness organic 


Regeneration density 
(stems ha -) 


Rosa acicularis cover 

Ledum groenlandicum cover 

log Ledum cover 

Rubus pedatus cover 

Epilobium angustifolium cover 
Cornus canadensis cover aay eas 
Rubus pubescens cover 
Calamagrostis canadensis cover 
Vascular plant cover 

log Vascular plant cover 
Hylocomium splendens cover 
Pleurozium schreberi cover 
log Pleurozium cover 


Polytrichum commune cover 


Total moss cover 
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Mean annual 
basal area 
increment 
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Table 18 


Elevation (m) 

log Elevation 

Slope percent 

Slope aspect 

Aspect X slope 

% Si+C in A horizon 


Thickness A horizon (cm) 


Thickness A (Si+C in A hor) 


%¢ Si+C in B horizon 
% Clay in B horizon 


Thickness B horizon (cm) 


Thickness B (Si+C in B hor) 


“Clay in A hor/Clay in B hor 


% Stone in profile 
Internal profile drainage 


Depth to distinct 
mottling (cm) 


Hydraulic conductivity 
(cm day ) 


log Hydraulic conductivity 
Structure of B horizon 
Consistence B horizon 
Chroma A horizon 

Hue B horizon 

Value B horizon 


Chroma B horizon 


Simple correlation coefficients between dependent 
variables and significant independent variables as 
determined by t-tests in multiple regression equations 


for the white spruce plots (n = 30) 


502 
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547) 


42 


-O5 
Vv 


Mean annual 
total volume 
increment 


ae 


41 


ex 


Stem 
density 


PAIL 
Site per 
index ree 
Sede ae. 
- .53 
- .02 
w 
+09 
x* 
- .41 
Vv 
= S28 


Age of 
maximum 


26 
v 


+24 


pee 


Mean annual 
basal area 
increment 


- ~44 
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Table 18 cont. 
Mean annual PAI Age of Mean annual 
fotal volume Stem Site per maximum basal area 
increment density index tree PAI increment 
Canopy cover ae gras, = ABs 
Litter cover 
kx xk 
l/log Litter cover - .40 + .49 
kk x* 
Deadfall cover - .43 = SA Ae Sete 
Stem density 
* 
Thickness organic horizon (cm) 1 Oe tie > 
log Thickness organic 
Regeneration density 
(stems ha 7) + .22 
Rosa acicularis cover => Als! eee Oe so Sty 
“ * wae 
Ledum groenlandicum cover - .26 — eas} aS 2 One) ei 80 
log Ledum cover 
Rubus pedatus cover 
Epilobium angustifolium cover =) 6 
Cornus canadensis cover = 235 a Oo) 
2 wa 
Rubus pubescens cover Te 
Calamagrostis canadensis « — 
cover - .10 + .16 - .04 
ak 
Vascular plant cover - .44 
xx 
log Vascular plant cover + .44 
Hylocomium splendens cover 
Pleurozium schreberi cover 
log Pleurozium cover 
Rex 
Polytrichum commune cover + .46 
Total moss cover 
log Total moss cover ee? 
Total lichen cover 
**k Significant at p = .01l 
** Significant at p = .05 
* Significant at p = .10 
v Independent variable occurring only in equations with 


vegetation variables deleted 
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Table 19 


Elevation (m) 

log Elevation 

Slope percent 

Slope aspect 

Aspect X slope 

% Si+C in A horizon 


Thickness A horizon (cm) 


Thickness A (Si+C in A hor) 


$ Si+C in B horizon 


% Clay in B horizon 


Thickness 


Thickness 


B 


B 


horizon 


(cm) 


(Si+C in B hor) 


Clay in A hor/Clay in B hor 


$ Stone in profile 


Internal profile drainage 


Depth to distinct 
mottling (cm) 


Hydraulic conductivity 
(cm day -) 


log Hydraulic conductivity 


Structure of B horizon 
Consistence B horizon 
Chroma A horizon 

Hue B horizon 


Value B horizon 


Simple correlation coefficients between dependent 
variables and significant independent variables as 
determined by t-tests in multiple regression equations 


for the black spruce plots (n = 15). 


Mean annual 
total volume 
increment 


Stem 


density 


PAI 
Site per 
index tree 
oe AGW) 
zt 
+ 653 
Vv 
+ .50 
* 
+7560 - .19 
wg Te 
Vv 
gee Ae; 
Vv 
= 29 
Vv 
$e 2O 
= sm ep Les 
a 
ee 


Age of 
maximum 
PAI 


Mean annual 
basal area 
increment 
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Table 19 cont. 


Mean annual PAI Age of 
total volume Stem Site per maximum 
increment density index tree PAI 
; xx 

Chroma B horizon + .82 —2554 aera) 

Canopy cover + 220 ao lke! a aes) 

Litter cover 

. xk* 

l/log Litter cover + 560 

Deadfall cover 27 

Stem density - 04, 


Thickness organic horizon (cm) 
log Thickness organic 


Regeneration density 
(stems ha ~) 


Rosa acicularis cover + .42 = .24 
Ledum groenlandicum cover - .46 

log Ledum cover 

Rubus pedatus cover =~ 14 sb? 
Epilobium angustifolium e Tea Bees 

cover aps dhs} =o acts) 
Cornus canadensis cover ye se) 
Rubus pubescens cover ao Oke! 
Calamagrostis canadensis cover =e, 
Vascular plant cover 

log Vascular plant cover 

Hylocomium splendens cover = ow = 5s 
Pleurozium schreberi cover 
log Pleurozium cover 
Polytrichum commune cover 
Total moss cover 


log Total moss cover a BA =o TA 


Total lichen cover + .39 


#** Significant at p = .01 


** Significant at p = .05 


" 
. 
~ 
OQ 


* Significant at p 
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Mean annual 
basal area 
increment 


Sal 


v Independent variable occurring only in equations with 


vegetation variables deleted 
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Table 20 Simple correlation coefficients between dependent 
variables and significant independent variables as 
determined by t-tests in multiple regression equations 


a for the aspen plots (n = 17). 


Mean annual PAIL Age of Mean annual 
total volume Stem Site per maximum basal area 


increment density index . tree PAI increment 


Elevation (m) 
log Elevation + .28 

Slope percent tae ae 

Slope aspect = 20 
Aspect X slope 

% Si+C in A horizon +35 
Thickness A horizon (cm) + .40 
Thickness A (Sit+tC in A hor) + ,43 

% Si+C in B horizon =e 

% Clay in B horizon 

Thickness B horizon (cm) 

Thickness B (Si+C in B hor) 

Clayean A hor/Clay, in) B hor =e LD + .24 

% Stone in profile =e Gee 
Internal profile drainage = 18 


Depth to distinct 
mottling (cm) - .29 


Hydrauli¢ conductivity 
(cm day ) 


log Hydraulic conductivity 

Structure of B horizon - .34 m5 2 + .29 = s0 
Consistence B horizon - ae ae Elly 

Chroma A horizon = 215 - .Ol1 

Hue B horizon = 55) 


Value B horizon Pee LE + oe mete SET, 


parity fa iinet 


ub Rs et So he 
* sa Ger, rh ci 
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Table 120 cont. 


Mean annual PAIL Age of Mean annual 
total volume Stem Site per maximum basal area 


increment GenSBtye se nde Xen Omer A increment 


Chroma B horizon =e Sill 

Canopy cover 

Litter cover 

l/log Litter cover 

Deadfall cover - .20 + .06 
Stem density - 19, 


Thickness organic horizon 
(cm) + .61 


log Thickness organic eS ~- .10 
Regeneration density (stems tee = Oe) 
Rosa acicularis cover 
Ledum groenlandicum cover =O 
log Ledum cover 
Rubus pedatus cover 
Epilobium angustifolium cover ihe 2. cee 
Cornus canadensis cover + eel +2 0 Piero) 
Rubus pubescens cover - .42 43 
Calamagrostis canadensis cover =.25 
Vascular plant cover a2 Cos! 

xx xxx 
log Vascular plant cover - .60 +525 a Ars ii/ 
Hylocomium splendens cover set Skt 
Pleurozium schreberi cover te il = .14 
log Pleurozium cover ’ tae 
Polytrichum commune cover 
Total moss cover 
log Total moss cover 


Total lichen cover 


*** Significant at p = .Ol 
** Significant at p = .05 
* Significant at p = .10 


v Independent variable occurring only in equations with 
vegetation variables deleted 
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1.1 Elevation (Xl, X2) 


Elevation or its log was a significant independent variable in seven 
equations: white spruce MAI (vegetation variables included and deleted), 
Site index (vegetation variables included and deleted) and mean annual basal 
area increment; black spruce site index; and aspen age of maximum PAI 
(vegetation variables deleted) (Tables 8, 9, 12, 14). All the correlations 
between the dependent variables and elevation were negative with the 
exception of white spruce mean annual basal area increment (Tables 18-20). 

The negative correlation between growth and elevation is expected and. 
is an indication of the cooler climate and shorter growing season at high 
elevations. It is notable that where the elevations of the sampled stand 
groups are within fairly narrow ranges, such as for black spruce and aspen, 
factors other than elevation contributed more to the Ro value in the 
regression equations (Tables 9, 10, 13, 14). The positive correlation of 
elevation with white spruce mean annual basal area increment may be a 
reflection of a tendency at high elevations, for volume growth to be added 
to the lower part of the tree (greater taper) as MAI is negatively correlated 
with elevation (Table 18). Other authors have shown similar relationships 
between productivity and elevation. For example, Mader (1963) demonstrated 
that five-year volume increment for red pine (Pinus resinosa Ait.) was 
negatively correlated with increasing elevation. Graney and Ferguson (1971) 
showed that shortleaf pine (Pinus echinata Mill.) site index in Arkansas 


was also negatively correlated with increasing elevation. 


1.2 Slope gradient (X3) 


Slope gradient was a significant independent variable in five equations: 
lodgepole pine density; white spruce MAI and site index, both with vegetation 


variables deleted; black spruce age of maximum PAI (vegetation variables 
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nye27 
deleted) and aspen age of maximum PAI (vegetation variables deleted) (Tables 
io. 2, 14). The simple correlation coefficients between slope gradient 
and the growth expressions were positive, with the exception of lodgepole 
pine density (Tables 17-20). 

The positive correlation of slope gradient with several of the dependent 
variables is in agreement with other workers. Dumanski et al. (1973) found 
that increasing topographic steepness, implying better soil drainage, was 
associated with increased PAI for lodgepole pine, with optimal slopes being 
approximately 15-30%, for the Hinton-Edson area of Alberta. Duffy (1964a) 
found that lodgepole pine height growth was positively correlated with slope 
gradient on till parent materials in the Rocky Mountain House area, Alberta. 
Della-Bianca and Olson (1961) also observed that height growth was positive- 
ly correlated with percent slope in Piedmont hardwood and pine-hardwood 
upland forests. Myers and Van Deusen (1960) noted similar relationships 
for ponderosa Sys eee index in the Black Hills area of South Dakota. 

The negative correlation of percent slope with lodgepole pine density 
is perhaps an indication of the difficulty of seedling establishment on 
slopes where steepness results in increasing erosion and where young seedlings 
Or pine seed may also be washed away. 

The positive correlation of slope angle with aspen age of maximum 
PAI (Table 20) is believed to be due to aspen’ high frequency of occurrence 
on nearly level to gently sloping, imperfectly drained sites, and subsequent 
better growth on these sites. In the lodgepole pine plots Slope is 
negatively correlated with the thickness of the Ae horizon, (r =-.25) 
reflecting the lesser amount of leaching possible on steep slopes as well 
as probable greater erosion on these slopes, particularly following forest 
fires prior to stand origin. In the white spruce plots, slope was positively 


correlated with better soil profile drainage. It is interesting to note 
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that the correlation between Slope and drainage is not Significant for 
the lodgepole pine plots (R = -0.06) suggesting that the productivity 
increases on steeper Slopes in the lodgepole pine forests may not be due 
to better drainage as was suggested by Dumanski et al. (1973) for west- 
central Alberta. In fact, steeply sloping sites, especially in the lower 
Slope position may actually be more moist than their more gently sloping 
counterparts due to groundwater discharge or seepage within the site. 
Unfortunately, the effect of Slope position on productivity was not 
evaluated in this study. A more extensive consideration of drainage is 
included in the separate discussion of drainage. 

In the black spruce plots, slope gradient was positively correlated 
with thickness of the A horizon (usually an Aeg horizon) and profile 
drainage. Though this correlation is opposite to that of the pine plots 
(© =+.50), it reflects the negative correlation of Slope with the thickness 
of the organic layer and the tendency for peat accumulation, and restricted 


Ae development on level black spruce sites. 


1.3 Slope aspect (X4) 


Slope aspect was a significant independent variable in six growth 
prediction equations: lodgepole pine MAI and PAI per tree; white spruce 
MAI and age of maximum PAI (vegetation variables deleted); black spruce 
Site index (vegetation variables deleted); and aspen mean annual basal 
area increment (Tables 7, 8, 10, 12, 13). In order to be entered as an 
independent variable in a stepwise multiple regression, aspect must be 
coded in a form other than as azimuth degrees due to the circular nature 
of the compass and the obvious difficulty of coding northerly aspects and 
of averaging aspects. This problem was overcome by using a method adopted 


by Myers and Van Deusen (1960) and several authors since, who expressed 
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aspect as (sine of azimuth clockwise from southeast)+1, where the south- 
east exposure was least productive. This method gives a range of trans- 
formed aspect from 0 to 2 and gives equal ratings to 0° and 360° azimuth. 
In this study, the most productive aspect was found to be between north 
and north-northeast (lifecie and the least productive were southeast to 
south-southeast aspects, when MAI values were stratified according to 
Slope aspect. Thus northerly aspect was then assigned a maximum value of 
2, (sine 90° +1). The difference between 11.3° and 90.0° is the correction 
used. The aspect transformation thus became aspect = sine (azimuth+78.0)+l. 
Thus, positive correlations with Slope aspect reflect a tendency towards 
maximum productivity on north to northeast aspects. In the equations 
mentioned above for lodgepole pine and black Spruce, the significant 
correlations of aspect with the dependent variable were positive (Tables 
rok). 

Lloyd and Lemmon (1970) found that the optimum rectification for oak 
site index in the Northern Appalacians was 45-60° depending on soil series. 
Site index was lowest in the west and highest in the northeast. They 
decided to use Sa. aS an average correction as opposed to 45° used by many 
authors, but emphasized that too much reliance should not be placed on the 
appropriate rectification for azimuth for any one soil series due to one 
Or more of several factors: (1) perhaps not enough plots were available 
to adequately sample the population in question, (2) the data may not have 
properly sampled all sectors of the azimuth circle, and (3) some undefined 
interactions may have affected the relationships. 

Several other authors have discovered preferences of various tree 
species for particular slope aspects. Duffy and England (1967) observed 
that sites on steep, north-facing slopes were more productive for a given 


surface material than those on south-facing slopes in the Kananaskis 
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Experiment Station, Alberta. Duffy and Knight (1967) noticed the same 
trends and also that the best forest growth occurred on upper north-facing 
slope positions, but upper and mid south-facing slopes and imperfectly 
drained hilltops were less productive in the Whitecourt,’ Alberta forest area, 
Dumanski et al. (1973) found that aspect had minimal effect on medium and 
fine-textured materials but the north aspects tended to outproduce south 
aspects on sandy and gravelly materials in the Hinton - Edson area, Alberta. 
Other studies in the eastern United States have also indicated better 
growth on northerly aspects (Carmean 1965; Graney and Ferguson 1971, 1972; 
Lee and Sypolt 1974). The more favorable moisture regime on northerly 
aspects is often given as the explanation. Lee and Sypolt (1974) from 
observations on adjacent north and south-facing slopes in West Virginia 
suggested that in areas of high summer precipitation, marked differences 
in growth are not associated with soil moisture differences. At midday on 
south-facing slopes it appears that: (1) super-optimum temperatures for 
net assimilation occur frequently in forest canopies (2) soil temperatures 
lower than canopy temperatures limit absorption of available water and (3) 
higher canopy temperatures increase transpiration demand, create leaf water- 
deficits earlier in the day and magnify the midday depression of net 
assimilation (Lee and Sypolt 1974). In contrast to the findings of the 
present study, Lowry (1975) found that black spruce site index in eastern 
Canada was higher on south aspects and on mineral soil, indicating that 
the cooler north aspects were limiting to black spruce growth. 

West-facing slopes are warmer than east-facing slopes. The maximum 
solar radiation on east-facing walls was observed two to three hours after 
dawn, when soils and air are cool. The maximum radiation on west-facing 


walls was observed two to three hours before sunset. Thus, the same 


radiation intensity would be expected to heat a west-facing slope Cova 
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higher temperature than an east-facing slope (Cunniff 1958 cited in 
Taylor and Ashcroft 1972 p. Sly a 

It has been suggested by some authors (Squillace and Bingham 1958 and 
Hermann and Lavender 1968 cited in Spurr and Barnes 1973 p. 40) that 
differences in aspect and thus soil moisture may lead to ecotypes found on 


certain aspects. 


1.4 _Aspect and slope interaction (x5) 

Because slope gradient and aspect act together to influence tree 
growth, it seemed desirable to include an interaction or product of 
aspect and slope as an independent variable within the regression. To- 
gether slope and aspect determine the amount of insolation at the ground 
surface which in turn determines air and soil temperatures, influences 
Precipitaticn distribution ard soil moisture coréitions. 

Surprisingly, the slope X aspect interacticn was a Significant indepen- 
dent variable in only two equations: black Spruce PAI per tree and mean 
annual basal area increment (Table 9). Both correlations were positive 
(Table 19) suggesting that black spruce basal area and PAI per tree are 
favored by steep north aspects. Dukovnikov (1975) used an integrated 
slope, aspect, relief and altitude index for Pinus sylvestris, Picea abies 


and Abies alba in Bulgaria, that effectively integrated these variables. 


2. Internal Site (Soil) Factors 


Included in internal site factors are the various soil physical 
Properties which influence or reflect differences in forest growth includ- 


ing soil texture, horizon thickness, drainage, color, structure and 


consistence. 
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2.1 Texture and. horizon thickness 

Several measures of soil texture used in the regression analysis 
including silt plus clay contents of the A (usually Ae) and B (usually Bt) 
horizons, a ratio of clay in the A horizon to clay in the B horizon and 
stone content of the profile are discussed. The thickness of the A horizon 
and interactions of thickness of the A and B horizons with silt plus clay 
content of the respective horizons are also discussed here. 

Silt plus clay content of the A horizon (Var. X20) was a Significant 
independent variable in seven equations; lodgepole pine age of maximum PAI 
(vegetation variables deleted); black spruce MAI, site index, age of maximum 
PAI (vegetation variables included and deleted) and mean annual basal area 
increment; and aspen age of maximum PAI (vegetation variables deleted) (Tables 
72, 20, 11,-13, 14). Lodgepole pine age of maximum PAI, black spruce 
Site index and age of maximum PAI were positively correlated with silt plus 
clay content (Tables 17, 19) while black spruce mean annual basal area 
increment was negatively correlated (Table 19). 

It appears that increases in silt plus clay content of the A horizon 
indicate more impeded drainage in the already poorly drained black spruce 
forests. The positive correlation with age of maximum PAI indicates slower 
maturation of black spruce on the fine textured soils where profile 
drainage is likely to be more impeded. 

Elsewhere, other authors have found correlations of A horizon silt 
plus clay content with forest productivity. Mader (1963) found that five- 
year volume increment for red pine was positively correlated with moisture 
holding capacity of the A and B horizons. Mogren and Dolph (1972) found 
for north-central Colorado and south-central Wyoming that lodgepole pine 
site index was negatively correlated with percent soil particles greater 


than 0.25 mm diameter in the A horizon, however there was no significant 
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positive correlation with silt plus clay contents in the A and B horizons. 
They believed that correlations with site index were accounted for by water 
availability. Mader (1976) discovered that higher silt plus clay contents 
in the A horizon were associated with better white pine sites in 
Massachusetts, probably as a result of better water availability and 
fertility. Lutz and Chandler (1946) emphasized that excessively fine or 
coarse soils are not productive. 

The interaction of A horizon silt plus clay content and A horizon 
thickness (Var. X22), expressed as a product of the two was significant in 
only one equation: aspen site index (Table 10). The simple correlation of 
this interaction with aspen site index is positive (Table 20), indicating 
aspen's apparent preference for deep, fine textured A horizons. 

Thickness of the A horizon (Var. X21) considered alone was a signifi- 
cant independent variable in only two equations: black spruce mean annual 
basal area increment and aspen PAI/tree (Tables 9, 10). The only significant 
positive correlation was with aspen PAI per tree (Table 20) and is probably 
a reflection of more soil profile leaching and better drainage exerting a 
favorable influence upon PAI. 

Several other authors have determined relationships between A horizon 
thickness and forest productivity. Duffy (1964a) found for the Rocky 
Mountain House area, Alberta that lodgepole pine dominant height and 
average height were positively correlated with the thickness of the Ae 
horizon, while basal area per acre, total volume per acre and merchantable 
volume per acre were negatively correlated. Holmes and Tackle (1962) 
found that lodgepole pine height growth in Montana was positively 
correlated to the depth of abundant rooting which usually coincides with 
the depth of the Ae horizon. Hansen and McComb (1958), in southeastern 


Iowa, Observed that height and diameter growth of green ash and bleck 
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walnut were inversely correlated with the degree of erosion of the Ae 
horizon. Della-Bianca and Olson (1961) observed for the Piedmont hardwood 
and pine-hardwood upland forests, that height growth was positively correlat- 
ed with thickness of the A horizon and with percent organic matter in the 
A horizon X horizon thickness. Carmean (1965) described black oak site 
index in Ohio as a function of the depth of the Al+A2 horizons and several 
other factors for medium-textured, well drained soils and for fine-textured 
soils with restricted internal drainage. Mogren and Dolph (1972) had 
Similar results and observed that lodgepole pine site index was positively 
correlated with solum depth and believed that site index was accounted for 
by water availability. 

Silt plus clay content of the B (usually Bt) horizon (Var. X23) was 
a Significant independent variable in only one equation, aspen density, 
(Table 10) and the correlation was not significant (Table ZO)r 

Elsewhere, other authors have related B horizon silt plus clay content 
to forest productivity. Fralish and Loucks (1975) reported that aspen site 
quality in Wisconsin was negatively correlated with sand content of the Bt 
and positively correlated with available water-holding capacity. Duffy and 
Knight (1967) found the best forest growth in the Whitecourt Forest area 
of Alberta on moderately well drained sandy clay loam till sites on upper, 
north-facing slopes. Some sites had alluvial sandy loam overlays. Page 
(1976) observed that white spruce and balsam fir site index in Newfoundland 
were negatively correlated with silt plus clay content in the lower solum. 
Gessel and Lloyd (1950) observed Douglas-fir site index increases in 
north-western Washington when soil textures ranged from coarse to 
medium textures. Douglas-fir site index also increased on soils under- 
lain by an impeding layer when precipitation increased up to a total of 


150 cm per year. Hannah and Zahner (1970) discovered that the site index 
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of jack pine and bigtooth aspen in Michigan was significantly greater on 
SOils developed on outwash sands with textural bands of unsorted till-like 
material. Red pine site index in Minnesota was 2.5 m greater on plots with 
textural B horizons or textural banding and was more closely related to the 
surface soil (top 25 cm) than to subsoil when textural bands were absent 
(Alban 1974). Mader (1976) discovered that increased -silt plus clay content 
in the B horizon was associated with reduced site quality in natural stands 
of white pine in Massachusetts, perhaps reflecting reduced aeration and 
poorer rooting. Cation exchange capacity is also greater in fine textured 
soils. It is thus apparent that no broad generalizations can be made upon 
the influence of soil texture alone on site quality. A coarse-textured soil 
may support excellent productivity if the moisture regime is favorable. 
Influence of soil texture must be discussed in terms of other soil properties 
if it is to be meaningful, and generalizations cannot be made for large 
areas. For example, a fine-textured soil which retains optimum soil 
moisture in one area may retain too much in another area where precipitation 
is much higher, with resulting decreased yields due to flooding or lack of 
oxygen for tree roots. 

The interaction of silt plus clay content in the B horizon with horizon 
thickness, expressed as a product of the two (Var. X26), has often proved to 
be a useful measure for evaluating site quality. In this study, this 
interaction was a significant dependent variable in three equations: 
lodgepole pine age of maximum PAI; white spruce age of maximum PAI and 
mean annual basal area increment (Tables 7, 8). The significant simple 
correlation of the interaction with white spruce age of maximum PAI was 
negative (Table 18) indicating more rapid maturity of white spruce on 
soils with a likely high available water-holding capacity. The negative 


correlation of the interaction with white spruce mean annual basal area 
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increment (Table 18) is likely indicative of impeded drainage and 
restricted oxygen availability for white spruce roots on sites with thick, 
fine-textured Bt horizons. Other authors have shown relationships between 
Site productivity and the interaction of B horizon silt plus clay content 
and thickness. Duffy (1964a) demonstrated that percent sijtiplusrecriayein 
the Bt X effective horizon thickness /10 was negatively correlated with 
lodgepole pine dominant height, average height, total volume and 
merchantable volume on one soil series in the Rocky Mountain House area, 
Alberta. The correlation with basal area per acre was positive. 

The ratio of clay in the A horizon/clay in the B (Var. X27) was used 
aS an independent variable to express the degree of leaching or profile 
development and to determine whether this process might be related to 
forest productivity. This ratio was a significant independent variable 
in six equations: lodgepole pine PAI per tree; white Spruce age of maximum 
PAI (vegetation variables included and deleted); black spruce site index 
and age of maximum PAI (vegetation variables deleted); and aspen density 
and site index (Tables 7-10, 12, 13). The only significant Simple correla- 
tion was with lodgepole pine PAI per tree and was positive (Table The 
suggesting that lodgepole pine PAI per tree is greatest where the Ae 
horizon is well developed or where a sandy aeolian overlay may be present 
with a consequent favorable rooting environment in these soils. No 
examples were found in the literature where clay A/clay B horizon ratios 
were used as independent variables in growth prediction equations as just 
described. 

Stone content (Var. X28) based upon a visual estimate of coarse 
fragment® in the solum greater than 1 cm diameter was a significant 
independent variable in three equations: white spruce MAI; black spruce age 


of maximum PAI (vegetation variables deleted) and aspen age of maximum PAI 
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mliantes 8, 10, 13). Only the negative correlation of stone content with 
black spruce age of maximum PAI was Significant (Table 19). It is generally 
regarded that stones in the solum or the subsoil do not have an adverse 
effect upon forest productivity unless there are so many that rooting 
volume and/or water-holding capacity is reduced. Lutz and Chandler (1946) 
suggested that root space becomes constricted when stone content is 
greater than 20%, but stone content makes fine-textured soils more porous 
with root penetrability facilitated. Also the presence of stones on and 
in the soil decreases evapotranspiration losses. In the Wapiti map area, 
stone content within the profile seldom exceeds 202% and averages from 73 
in the pine plots to less than 1% in the aspen plots (Tables 3-6). Black 
Spruce age of maximum PAI is negatively correlated with stone content 
(Table 19) but is positively correlated with silt plus clay content in the 
A horizon (Table 19) indicating that stone content is a reflection of 
coarse texture. Mogren and Dolph (1972) observed that lodgepole pine site 
index in north-central Colorado and south-central Wyoming was negatively 
correlated with percent surface stoniness but there was no Significant 
correlation with percent stone in the solum. They attributed the 


relationship with site index to be a reflection of water availability. 


2.2 Soil Profile Drainage 

Soil profile drainage and site moisture regime are determined by 
topographic position, the physical nature of the parent material, depth and 
position of water table, soil profile features and precipitation (Farrar 
1960). Three criteria were used as independent variables to characterize 
internal soil profile drainage; internal soil drainage class on a five 


point scale (Appendix 2), depth to distinct mottling (canada Soil Survey 


Committee 1970), and hydraulic conductivity as inferred from the texture of 
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the least permeable horizon (Appendix 2). Only one profile drainage- 
related variable was allowed to enter any one regression equation, namely 
the variable that accounted for the greatest increase in the R? value. 

Internal soil drainage (Var. X29) where 1 = poor and 5 = rapid, was 
@ significant independent variable in eight equations: lodgepole pine site 
index and mean annual basal area increment; white spruce site index; black 
spruce MAI (vegetation variables deleted), density, age of maximum PAI 
and site index (vegetation variables deleted) and aspen MAI (vegetation 
variables deleted) (Tables 7, 8, 97 137714). Wone of the Simple correla- 
tions of internal soil drainage with the dependent variables were signifi- 
cant. 

Jeglum (1974) reported increased black spruce growth on northern 
Ontario peatlands as soil aeration improved and that moisture-aeration and 
nutrient regimes explained large proportions of the variation in minor 
vegetation. Similar results were observed by Lowry (1974) for black spruce 
in eastern Canada. 

Depth to distinct mottling (Var. X30), a morphological criterion 
useful in assessment of soil internal drainage, was a significant indepen- 
dent variable in seven equations: lodgepole pine MAI and PAI per tree; 
black spruce MAI (vegetation variables deleted), density, site index 
(vegetation variables deleted) and age of maximum PAI (vegetation variables 
deleted); and aspen PAI per tree (Tables 7, 9, 10, 11). None of the simple 
correlations of depth to distinct mottling with the dependent variables 
were Significant. 

Hydraulic conductivity (Var. X31) or its log (Var. X32), depending 
upon which gave the better correlation with the dependent variable, was 
a significant independent variable in seven equations: white spruce MAI 


(vegetation variables included and deleted) and site index (vegetation 
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variables deleted); black spruce MAI (vegetation variables deleted) and 

PAI per tree; and aspen age of maximum PAI (Tables 7149, 207 elT ) Tale, 
Hydraulic conductivity was Significantly correlated only with white spruce 
density and PAI per tree (both negative; Table 18). The inferred hy@raulic 
conductivity is very highly correlated with the silt Plus clay content of 
tke B hcrizon since hydraulic conductivity is inferred from texture 
(Appendix 2). 

Other authors have found relationships between soil Crainage and 
moisture availability with forest productivity. Griffith (196C) foune 
that available soil moisture in the B horizons during the growing season 
(May to September) was the most important single variable affecting 
coastal Douglas-fir growth, accounting for 47% of the total variability. 
Duffy (1964a) discovered that lodgepole pine total volume per acre was 
negatively correlated with available moisture on a weight basis in the 
surface 91 cm (36 in) on one soil series in the Rocky Mountain House 
area of Alberta. Similar results were found by Holmes and Tackle (1962) 
in Montana. Spurr and Barnes (1973 p. 95) pointed out that the factors 
that contribute towards making a soil well drained, especially convex 
surfaces, are apt to ensure good drainage of cold air and indicate a site 
with lower maximum and higher minimum temperatures than would be indicated 
by the regional climatic average. They further noted that well drained 
soils support plants of southerly distribution. These trends were observed 
in the Wapiti map area, and are especially evident when the vegetation of 
the poorly drained black spruce - Labrador tea - cloudberry type (Type 7), 
which is comprised of species of predominantly northern distribution, is 
compared with the rapidly drained lodgepole pine - blueberry - lichen type 


(Type 14) which is comprised of species that have ranges extending much 
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farther south than the species of the black spruce bog forests (Type 7). 
Productivity is limited in the black spruce bog forests (Type:7) due in 
part to excess water and is limited in the lodgepole pine - blueberry - 


lichen type by insufficient soil moisture. 


2.3 Soil Structure (X33) and Consistence (X34) 

Structure and consistence of the B (usually Bt) horizon were entered 
as independent variables into the regressions in order to determine if 
these properties had utility in predicting the various dependent variables. 
The scales used to code these variables are in Appendix 2. Structure of 
the B horizon (1 = structureless, 7 = platy), implying a gradient from 
1 to 7 of increased ease of root penetration, was a Significant independent 
variable in nine equations: lodgepole pine density; white spruce age of 
maximum PAI (vegetation variables included and deleted) and mean annual 
basal area increment; black spruce MAI; and aspen MAI, density, mean annual 
basal area increment and age of maximum PAI (vegetation variables deleted) 
(Tables 7, 8, 9, 10, 12, 14). None of the simple correlations of structure 
of the B horizon with the dependent variables were significant. 

Soil porosity in part, is a result of structure that allows gaseous 
exchange with tree roots as well as providing rapid water infiltration. 
Root penetrability is best into soils with aggregates in the 2-6 mm diameter 
range (Taylor and Ashcroft 1972). The type of soil aggregation is 
important to soil permeability. A well-aggregated fine-textured soil may 
function as a coarse textured soil. Structure is most important in soils 
high in silt and clay and favorable structure is promoted by the presence 
of organic matter (Spurr and Barnes 1973). Tree roots enhance soil 
structure by compressing soil particles into aggregates and decaying plant 


roots tend to bind the aggregates. Favorable soil structure is best 
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maintained by healthy forest stands of Species well adapted to their 
environment (Lutz and Chandler 1946). Soil structure may be destroyed 

by the removal of litter and by leaving the soil surface exposed to the 
pounding of rain that causes compaction and leaching (Spurr and Barnes 1973, 
Lutz and Chandler 9 46)%. 

Soil consistence (Var. X34) was also considered in the regression 
analysis, and was expressed on a five point scale (Appendix 2) for the B 
(usually Bt) horizon ranging from loose to extremely firm. Soil 
consisténce is strongly influenced by soil moisture content that in turn 
is affected by the time of year, with soils usually being more moist in 
the early part of the growing season, than later in the summer. For this 
reason, soil consistence as used here might not be as strong an indicator 
of site productivity as some of the more permanent soil properties such 
as texture. Nevertheless, consistence of the B horizon was a significant 
independent variable in five equations: lodgepole pine MAI and Site index; 
and aspen age of maximum periodic annual increment (vegetation variables 
included and deleted) and site index (vegetation variables deleted) (Tables 
7, 10, 14). Only aspen site index (vegetation variables deleted) was 
Significantly correlated with consistence of the B horizon (Table 20). 

The negative correlation was interpreted as being a result of increased 
resistance to root penetration as consistance ranges from loose to 
extremely firm. 

This is the only study the author is aware of where soil structure 
and consistence were introduced as independent variables in a stepwise 
multiple regression approach to evaluating site quality. The early 


indications are that these variables might prove useful. 
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2.4 Soil Color 

Color is probably the most easily determined of soil characteristics. 
Several inferences about other soil characteristics may be drawn from 
soil color. Dark colors (low value and chroma) may be indicative of 
undecomposed organic materials or in soils low in organic matter, may 
indicate compounds of iron and humus, elemental carbon, or compounds of 
manganese, and magnetite (Soil Survey Staff 1951). Red colors in soil are 
generally related to unhydrated iron oxide which is relatively unstable 
under moist conditions and thus indicates good drainage and aeration. The 
light gray of Ae horizons may indicate a very low content of organic matter 
and iron (Soil Survey Staff 1951). Yellow colors are largely due to the 
presence of hydrated iron oxides frequently found in imperfectly drained 
soils. Dark gray colors may result from poor aeration. The gray color of 
gleyed horizons of poorly drained soils may be due to the presence of 
ferrous compounds. Soil color is inherited from the parent material and 
reflects it's composition and lithology (Soil Survey Staff 1951, Lutz and 
Chandler 1946). 

Four attributes of soil color from Munsell notation were entered as 
independent variables into the regression analysis: hue of the B horizon, 
value of the B horizon, chroma of the A horizon and chroma of the B horizon. 

Hue is the dominant spectral color of the soil. Within the range for 
the hue YR (yellow-red), the hue becomes more yellow and less red as the 
hue numbers increase. Within the range for the hue Y (yellow), the hue 
becomes more yellow and less green as the hue numbers increase (Soil 
Survey Staff 1951). For the Wapiti area, within the YR hues, a decrease 
in hue number (increasing redness) generally corresponds to a change in 
internal drainage class from imperfectly to rapidly, increased aeration 


and increased oxidation of iron compounds to the ferric form, rather than 
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with a change in color associated with different parent materials. 

Hue of the B horizon (Var. X36) was a Significant independent vari- 
able in only five equations: white spruce MAI (vegetation variables deleted); 
black spruce mean annual basal area increment and age of maximum PAI 
(vegetation variables deleted) and aspen age of maximum PAI (vegetation vari- 
ables included and deleted) (Tables 97-107 27 13°, 8L4). Only the negative 
Simple correlation with aspen age of maximum PAI was significant =(Table 20). 
In the soils of the sampled plots, the hues of the B (Bt, Br, Bg eta, Btgj, 
Btj) horizons are almost exclusively within 7.5 YR7vLO YR*=and =25:5-Y charts, 
the 2.5 Y colors indicating more restricted drainage and less oxidation of 
iron compounds than does the 10 YR chart. The negative correlation with 
aspen age of maximum PAI indicates more rapid maturity on the soil with 10 
YR hues (better drainage). Since hue of the B horizon seems to be largely 
a function of drainage for the soils of the forest stands sampled here, 
the author believes that soil internal drainage class is probably a better 
representation of soil moisture regime. Drainage class is assigned partly 
on the basis of soil hue, and the drainage class scale used here (Appendix 
2) permits evaluation of a much wider range of soil drainage conditions 
than does hue alone. 

Color value of the B horizon (Var. X37), an expression of the relative 
lightness of the color, was a significant independent variable in eight 
equations: lodgepole pine MAI and mean annual basal area increment; white 
spruce density and age of maximum PAI (vegetation variables deleted); aspen 
MAI, PAI per tree and mean annual basal area increment; and black spruce 
MAI (vegetation variables deleted) and site index (vegetation variables 
deleted) (Tables 7, 8, 10, 12, 13). The only significant correlation of 
color value of the B horizon was with black spruce PAI per tree (vegetation 


variables deleted) and was positive (Table 19). Increasing value or lightness 
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is usually a reflection of more intense weathering or leaching, although 
the mineralogic composition of the soil is also an important determinant. 
Soil color value darkens by one-half to three steps between dry and moist 
(Soil Survey Staff 1951). Because all of the soil colors were determined 
in the field at close to field Capacity moisture conditions, the variation 
among hue, value and chroma due to differences-in soil moisture should be 
minimal. Soils under lodgepole pine forests often have well developed Ae 
horizons of high color value, often to the extent where a secondary Ae 
horizon develops within the first,giving the "bleached" condition 

(Dumanski et al. 1973). High soil color value of the Bt horizon may be an 
indicatton of degraGing lodgepole pine site quality. Dumanski et al. 

(1973) found Bisequa Gray Luvisols (now Brunisolic Gray Luvisols), "Bleached" 
Orthic Gray Luvisols and Orthic Humo-Ferric Podzols (i.e. soils with high 
color value in the Ae horizon(s)) to have "high@productivaity for lodgepole 
pine in the Hinton-Edson area of Alberta. However, it does not necessarily 
follow that high color value in the Ae horizon coincides with high color 
value in the Bt. It must be conceded that in this study the contribution 
of organic matter or soil mineralogical composition to soil color value 
could not be evaluated, and may be important. 

Chroma of the B horizon (Var. X38) was a significant independent 
variable in six equations: white spruce density and PAI per tree; black 
Spruce density, PAI per tree and age of maximum PAI; and aspen site index 
(Tables 8, 9, 10). Only the positive correlation with black spruce density 
was Significant. This correlation may indicate the greater likelihood of 
high black spruce density on soils with better aeration and consequent high 
color values through increased oxidation of iron compounds. Again, the 


effect of mineralogical composition and organic matter cannot be assessed 


here. 
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Chroma of the A horizon (Var. X35) was a significant independent vari- 
able in seven equations: lodgepole pine MAI (vegetation variables deleted) ; 
black spruce MAI (vegetation variables deleted) ,PAI per tree and age of 
maximum PAI (with aA without vegetation variables); and aspen MAI and 
Gensity (Tables-9, 10, 11, 13). Only the positive correlation with lodgepole 
Pine MAI (vegetation variables deleted) was Significant (Table 17). Presum- 
ably, the differences in productivity are accountable in the same manner as 
those for chroma of the B horizon. 

It seems likely that further investigation into the factors determin- 
ing soil color would be worthwhile, and that these factors as well as the 
various expressions of soil color might serve as useful indicators of 


forest site potential productivity. 


3. Stand Factors 

Stand factors are characteristics of the forest stand that are functions 
of the physical site or internal site (soil) factors and of the site 
productivity. Stand attributes used as independent variables in the 
regression equations include forest canopy cover, tree stem density, litter 


cover, deadfall cover and thickness of the organic horizon. 


3.1 Cancpy Cover (X60) and Stem Density (X64) 

Canopy cover or extent of forest crown closure as determined from 
ocular estimates, was a significant independent variable in nine equations: 
lodgepole pine MAI, PAI per tree, age of maximum PAI and mean annual basal 
area increment; white spruce MAI and mean annual basal area increment; and 
black spruce density, site index and PAI per tree, (Tables 7, 8, 9). The 
only significant correlations occurred between canopy cover and lodgepole 


pine MAI, PAI per tree and age of maximum PAI (Table 17), and seem to 
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indicate that Canopy cover is a good indicator of site productivity in 
mature lodgepole pine forests. 

Stem density (Var. x64) expressed as stems per hectare was a signifi- 
cant independent variable in only three equations: lodgepole pine age of 
maximum PAI; black spruce age of maximum PAI (vegetation variables deleted), 
and aspen age of maximum PAI (vegetation variables deleted) (Table sts) 4); 
However, none of the simple correlations between stem density and the 
dependent variables were significant. Stem density is a function of several 
factors including fire history, insect or disease outbreaks, site moisture 
regime and seed availability at the time of stand origin. No attempt was 
made to determine the causes of stand density; only the effect of density 


upon productivity was evaluated. 


3.2 Organic Matter and Deadfall (x63) 

Organic matter here refers to two of the independent variables used 
in the regression analysis: litter cover and thickness of Organic horizons. 
The organic horizon includes the LFH and/or other Organic horizons. It 
seemed more appropriate to include thickness of Organic horizon in this 
section rather than in the discussion of internal site (soil) variables 
because it is often an indication of conditions at the ground surface, 
and is largely determined by the nature of the forest the site supports. 

It is thus closely related to the other stand factors discussed here, namely 
litter and deadfall. 

Litter refers to the small dead plant material on the grcund surface 
and includes conifer needles, deciduous tree leaves, small twigs and branches 
and understory plant material. Litter yield is dependent upon several 
factors including climate and latitude (yield greater in tropics), altitude 


and exposure (yield generally greater at low to mid-elevations and was found 


ab Vl eleatedg ‘Ole Fe ee ee vii 
re. rN Be 1. he be 3 

sotual oon Mamie ill dscns Cae 

tae al ele gee aaate, hat ni abe 

ere ase 0 stove ibe he ous iia 

eee ee ist wlan 

44 ba Why whi Nene tae aehie te a6 sdigete saat 
porno (ete Sib ideale S36¥ |. 


oe f 


7842 


> pap gab bart saicas Jotdiet # ea 
ito bate Se ahs ore eee 
PeATggee teiade ies ehevor ett! si 


, Maia ewe? 


tse) Lswtan tot’ 


Net) Seiberpphyhh” mag ae ae oF c ster Sesti ere 
: 
ia cewltheht pea) deme tetas japayeeeh 


7 


dnneks sco) dexeak eet See ehbetSal nonesct g 


Mi 


den | Sty Sho GP “one earl a rotiaes 
lida ob). (hte Tiereetge Sivitanenths wis wt: ae 
) Sencsp ont Ty SOR Pe edsenlenk ne Heat 6 28 

in ot) Orage), a oe ere: dty yo Denletersh: ‘i i 

jes den i Naateraing Selb, eames tata Od tredalex pisvoty “audd | 

“it sé Laan daaty feeb Biase’ ath of sever satini | Pa 
wc Bab cy le) Otpage Meee ae Greggs eeloeen ste aatieton 
Lathes nt ity Seng, «Leg aaoge gnelg yaodetebaw) a 

sfudatic i<v~nout Ai aetony Bal) Beige et Gee wii co eabeetane 


baud} &8% nc = tle_wedeshe be veal ve “TSAI Be vi (ssenep biolyt’ 


147 


to be higher on NE aspects than on SW aspects in German coniferous forests), 
fertility (greater on fertile sites) and soil moisture (yield decreases 
from mesic to dry) (Bray and Gorham 1964). Because the relationship 
between the dependent variables and litter cover appeared nonlinear on 

the scatter diagrams, the litter cover values were transformed as Iflog 
litter cover. Thus a negative correlation of a dependent variable with the 
1/log litter cover variable reflects a positive correlation with litter 
cover. 

The 1/log litter cover variable (Var. X62) was a significant indepen- 
dent variable in seven equations: lodgepole pine MAI, site index, age of 
maximum PAI and basal area annual increment; white spruce site index and 
age of maximum PAI; and black spruce age of maximum PAI (Tables 7, Bip DO ).s 
The 1/log litter cover variable was negatively correlated with lodgepole 
pine MAI and site index as well as black Spruce site index (Tables 17, 19). 
The other correlations were positive (Table 18). The negative correlations 
of l/log litter with the MAI and site index terms are actually positive 
with actual litter cover, indicating that sites that are productive in 
terms of volume and height growth also have high foliage production which 
in turn is reflected by high needle fall and subsequent accumulation. 

Bray and Gorham (1964) have also indicated that leaf litter may often be 

an index to net production in coniferous forests. Litter tends to 
accumulate in cold, acid conditions of boreal forests (Spurr and Barnes 
1973 p. 182). Yuasa et al. (1974) observed that Japanese red pine needle 
fall was greater on better sites. It is interesting to note that the 1l/log 
litter cover variable was not significant in any of the aspen equations, 
perhaps reflecting the more rapid decomposition and lower accumulation 
rates of the deciduous leaf litter. Gymnosperms have been shown to 


yield up to one-sixth more total litter annually than angiosperms (Bray 
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and Gorham 1964). The l/log litter cover variable was positively 
correlated with age of maximum PAI in the lodgepole pine, white spruce and 
black spruce plots, reflecting more rapid maturity (maximum volume increment 
at an earlier age) on sites with lesser litter cover. The correlation of 
the 1/log litter cover variable with lodgepole pine mean annual basal area 
increment is positive (negative with actual litter cover), opposite to 
that with lodgepole pine MAI. indicating that the increments in pine volume 
growth are reflected in terms of height growth more so than stem diameter 
growth. The author is not familiar with other studies where litter cover 
has been used as an independent variable in growth prediction equations. 
Its use in this study indicates it might be a reliable indicator of produc- 
tivity in site studies of coniferous species. 

The thickness of the organic horizons (Var. X39) above mineral soil, 
Or the log of its thickness (Var. X40) in the cases of two lodgepole pine 
equations (Table 7) and two aspen equations (Table 9) was a Significant 
independent variable in nine equations: lodgepole Pine density, age of 
maximum PAI (vegetation variables deleted), mean annual basal area increment; 
white spruce PAI per tree and age of maximum PAI; black spruce mean annual 
basal area increment; and aspen MAI, site index and PAI per tree (Tables 7-11). 
The correlations of organic horizon thickness with lodgepole pine age of 
maximum PAI (vegetation variables deleted), white spruce age of maximum PAI 
and aspen MAI were positive (Tables 17, 18, 20), that with black spruce 
mean annual basal area increment was negative (Table 19). The remainder of 
the simple correlations were nonsignificant. 

In the lodgepole pine forests,the organic horizons are primarily pine 
needles and feathermosses and the positive correlation between organic 


horizon thickness and age of maximum PAI (vegetation variables deleted) 
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indicates later maturity in stands with thicker LFH layers. In the white 
spruce forests, the organic horizon often has a different composition than 
that in the lodgepole pine forests and is frequently developed under differ- 
ent conditions, especially where restricted drainage or seepage occurs. The 
Organic horizon is thickest in the white spruce forests on Gleysolic soils 
with humified or very well-decomposed organic matter and on soils with 
peaty, fibric, organic horizons. Conditions of restricted drainage and 
limited aeration may account for the negative correlations of organic 
horizon thickness with black spruce mean annual basal area increment and 

the positive correlation with white spruce age of maximum PAI, indicating 
later maturity on poorly drained soils. Page (1976) observed that white 
spruce and balsam fir site index in Newfoundland were negatively correlated 
with humus depth. Lowry (1975) reported that black spruce site index was 
higher on eastern Canada sites where the thickness of the FH layer is less 
than 20 centimeters. The organic horizon in the aspen forests is composed 
primarily of deciduous leaf litter and is usually less than 10 cm thick 

with an average of 7.5 cm (Table 6). Aspen MAI is positively correlated 
with organic horizon thickness, indicating deep litter accumulation on highly 
productive sites. 

Deadfall (Var. X63), or dead trees lying on the ground as a result of 
natural mortality through insect damage, disease, windthrow or fire within 
the current forest generation, expressed as percent cover of the ground 
surface, was a significant independent variable in seven equations: lodgepole 
pine MAI; white spruce site index, PAI per tree and age of maximum PAI; black 
Spruce site index; and aspen site index and age of maximum PAI (Tables 7-10). 
Only the negative correlation of deadfall with white spruce site index and 
the positive correlation with white spruce age of maximum PAI (Table 18), 


are significant. The negative correlation with white spruce site index 
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indicates some event in the stand's history which has resulted in mortality 
and subsequent windthrow. Such natural stand thinning may "release" 
residual trees. The positive correlation between deadfall cover and white 
Spruce age of maximum PAI, may also be an indicaticn of past events which 


resulted in a suppression of the growth rate of the white spruce. 


4. Vegetation 

The discussion of vegetation as independent variables in the 
regression is separated into four sections: tree regeneration, shrubs, herbs, 
and mosses and lichens. Vegetation variables accounted for a large percent- 
age of the R? value in many of tke equations in Tables 7 to 10. The influence 
cf the vegetation variables on the R? value is apparent when Tables 7 to 
10 are compared to Tables 11 to 14. The R? value for the lodgepole pine 
MAI equation for example Ccrcps from 0.€6 tc C.24 when the vegetation variables 
are deleted. This reduction in R? value was much more evident on the lodge- 
pole pine and white spruce equations than for the black spruce or aspen 
equations. The drop in F ratio for the black spruce and aspen equations 
was more significant, to the point that the F ratios fcr black spruce MAI, 
aspen MAI and aspen site index were not significant at the p = 0.05 level. 


A discussion of each of the vegetation variables follows. 


4.1 Tree Regeneration Density (X80) 


Tree regeneration includes all tree species less than 1.5 m tall and xs 
expressed as stems per hectare. Summation of tree regeneration for all 
species was enterec fcr each plot in the regression analysis. Regenera- 
tion density was a significant independent variable in five equations: 
lodgepole pine tree density, site index and PAI per tree; white spruce MAI; 


and aspen tree density (Tables 7, 8, 10). The correlation with lodgepole 
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pine tree density was negative (Table 17), indicating conditions 

favorable to seedling establishment in Open stands. The positive correla- 
tion of regeneration density with lodgepole pine PAI/tree (Table 17) 
probably indicates greater PAI growth in more open stands, also a function 


of tree density. 


4.2 Shrub Cover 

The cover values of two shrubs, Ledurm groenlandicum and Rosa acicularis 
were included as independent variables in the regression analysis. Ledum 
groenlandicum is a common understory species in moist to wet boreal forests, 
usually where the soil surface is acidic. Heringa and Cormack (1953) found 
Ledum most abundant on till soils in west-central Alberta. Ledum has its 
greatest cover in the black spruce plots (5.9%: Table 5) followed by lodge- 
pole pine (4.1%: Table 3), white spruce (0.4%: Table 4) and aspen stands 
(0.1%: Table 6). Ledum groenlandicum cover (Var. X82) or the log of Ledum 
cover (Var. X83) in the case of lodgepole pine stem density, was a 
Significant independent variable in eight equations: lodgepole pine density; 
white spruce MAI, site index, PAI per tree, age of maximum PAI and mean 
annual basal area increment; black spruce site index; and aspen age of 
maximum PAI (Tables 7-10). The significant correlations with pine density 
and white spruce age of maximum PAI and mean annual basal area increment 
were positive, while that with black spruce site index was negative 
(Tables 17-19). The positive correlation with white spruce age of 
maximum PAI indicates later maturity on sites with high Ledum groenlandicum 
cover. The positive correlation with mean annual basal area increment 
indicates that white spruce trees have more taper on sites with higher 


Ledum groenlandicum cover. 


Rosa acicularis is another common shrub species and is most abundant 
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in open areas, especially on disturbed sites at low elevations, usually 
favoring the drier site. Rosa cover (Var. X81) is 3% in the aspen plots 
(Table 6) but is 1% or less for the other tree species (Tables 3 to Ds). 

Rosa acicularis cover was a Significant independent variable in seven 
equations: lodgepole pine density and mean annual basal area increment; 
white spruce MAI, site index and mean annual basal area increment; and 
black spruce MAI and age of maximum PAI (Tables 7-9). The only significant 
correlation with Rosa acicularis cover was with black spruce MAI (Table Loy 
and is likely an indication of better growth and earlier maturity with 


improved soil drainage as reflected by higher Rosa acicularis cover. 


4.3 Herb Cover 

Herbs include all vascular plants with soft, non-woody aerial portions 
which die kack to the rootstOck each fall. They include forbs (broad 
leafed) and graminoids (grasses and sedges). Five common species of 
herbaceous plants were included in the regression as independent variables: 
Cornus canadensis, Rubus pubescens, Calamagrostis canadensis, Epilobium 
angustifolium and Rubus pedatus. 

Cornus canadensis or Canada bunchberry is a common understory species 
in the boreal forest. It has its highest cover in the white spruce forests 
(8.6%: Table 4) followed by the lodgepole pine forests (8.33%: Table EB 
black spruce forests (5.8%: Table 5) and aspen forests (5.6%: Table 6). 
Cornus canadensis cover (Var. X86) was a significant independent variable 
in ten equations: lodgepole pine MAI, site index and mean annual basal area 
increment; white spruce density and site index; black spruce site index 
and mean annual basal area increment; and aspen density, PAI per tree and 
age of maximum PAI (Tables 7-10). The significant simple correlaticns 


were between Cornus and lodgepole pine MAI, lodgepole pine site index, 
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black spruce site index and mean annual basal area increment and were 
positive (Tables 17, 19). In general the Presence of Cornus canadensis 
seems to indicate favorable site conditions for the four tree species 
Sampled. Illingworth and Arlidge (1960) found that Corrus Canadensis 
characterized the best lodgepole pine sites in central British Columbia. 

Rubus pubescens is less common than Cornus canadensis in the boreal 
forest, being most abundant at lower elevations on moist Sites. It has its 
highest cover in the aspen forests (2.2%: Table 6) follcwed by the white 
Spruce forests (1.7%: Table 4), lodgegcle pine forests (0.4%: Table 3) and 
black spruce forests (0.3%: Table 5). Rubus pubescens cover (Var. X87) was 
a Significant independent variable in five equations: lodgepole pine site 
index; white spruce PAI per tree; black spruce site index; and aspen MAT 
and mean annual basal area increment (Tables 7-10). The simple correlations 
between Rubus pubescens and the dependent variables were positive and 
Significant (Tables 17, 18, 20) with the exceptions of aspen MAI that was 
negative (Table 20) and black spruce site index (non-significant) , 
indicating that even though Rubus pubescens is most abundant in the aspen 
forests it is not necessarily a good indicator of favorable aspen site 
conditions. 

Calamagrostis canadensis or marsh reed grass is common on moist sites 
especially in open areas and is capable of rapid increases in cover follow- 
ing drastic disturbances such as clearcutting (Corns and LaRoi 1976). 

Marsh reed grass has its highest cover in the aspen forests (0.4%: Table 6) 
followed by the black spruce forests (0.2%: Table 5), lodgepole pine and 
white spruce forests (0.1%: Tables 3, 4). Calamagrostis canadensis cover 
(Var. X88) was a significant independent variable in five equations: white 
Spruce MAI, PAI per tree and mean annual basal area increment; black spruce 


PAI per tree; and aspen PAI per tree (Tables 8-10). All correlations 
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between marsh reed grass cover and the dependent variables were non- 
Significant. Tllingworth and Arlidge (1960) found for central British 
Columbia that the lodgepole pine - Calamagrostis type was second only to 
the lodgepole pine - Cornus - moss type with respect to productivity. 

Epilobium angustifolium or fireweed is common in open areas in the 
boreal forest, particularly on recently disturbed or cleared sites after 
fire or clearcutting. It has its highest cover in the aspen forests 
(3.1%: Table 6), followed by the lodgepole pine forests (leles "Table 3). 
white spruce forests (0.3%: Table 4), and black spruce forests (OL 23: 
Table 5). Fireweed cover (Var. X85) was a significant independent variable 
in six equations: white spruce PAI per tree; black spruce MAI, site index 
and mean annual basal area increment; and aspen density and site index 
(Tables 8-10). Only the Simple correlations between fireweed cover and 
black spruce MAI, site index, and mean annual basal area increment and 
aspen PAI per tree were significant and all were positive (Tables 19, 20). 
The positive correlations with the black spruce dependent variables 
probably reflects the better black spruce growth on more well drained sites 
where fireweed is likely to be more abundant. The positive correlation 
with aspen site index is presumably an indication of greater height growth 
in more open stands. 

Rubus pedatus is a species that is often abundant in lodgepole pine- 
black spruce forests at high elevations in the foothills, usually in 


association with Vaccinium membranaceum and appears to be an indicator of 
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cool, moist conditions. Rubus pedatus has its highest cover in the lodgepole 


Pine forests (1.1%: Table 3) followed by the black spruce and white spruce 
forests (0.1%: Tables 5, 4) and is absent in the aspen forests. Rubus 
pedatus cover (Var. X84) was a significant independent variable in only 


two equations: black spruce MAI and PAI per tree (Table 9). Neither of 
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the simple correlations was significant. 

Total vascular plant cover or the sum of the cover values of all 
understory shrubs and herbs was entered into the regressions in addition to 
the cover values of the individual plant species. Vascular plant cover 
was highest in the lodgepole pine forests (41.3%: Table 3) followed by 
the aspen forests (40.3%: Table 6), white spruce forests (34.1%: Table 4) 
and black spruce forests (31.7%: Table 5). The abundance of vascular plant 
cover seems to be related to the degree of crown closure and to site moisture 
regime, with generally higher vascular plant cover on the open moist sites. 
Heringa and Cormack (1953) related plant cover to texture and observed 
that most species in lodgepole pine stands of central Alberta were most 
abundant on silt loam textures. This effect may be a reflection of 
drainage, as drainage and texture are closely related as discussed earlier. 

Vascular plant cover (Var. X89) or its log (Var. X90), depending on 
which gave the closet linear relationship to the dependent variable in 
question, was a significant independent variable in nine equations: lodge- 
pole pine density, PAI per tree and mean annual basal area increment; white 
spruce density and PAI per tree; and aspen density, site index, PAI per 
tree and mean annual basal area increment (Tables 7, 8, 10). The simple 
correlations between vascular plant cover and pine mean annual basal area 
increment, white spruce density and aspen density were negative, the re- 
mainder positive (Tables 17, 18, 20). Only the correlations with pine 
mean annual basal area increment and aspen site index were non-significant. 
These observations can be accounted for by differences in canopy closure 


and site moisture regime as discussed above. 


4.4 Mosses and Lichens 


Three common mosses, total moss cover and total lichen cover were 
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entered into the regression as independent variables. The mosses were 
Hylocomium splendens, Pleurozium schreberi and Polytrichum commune. 

Hylocomium splendens is a common feathermoss in boreal and subalpine 
forests. Its cover was greatest in the white spruce forests (21.9%: Table 4), 
followed by the black spruce forests (20.9%: Table 5), aspen forests 
(5.1%: Table 6) and lodgepole pine forests (5.0%: Table 3). Hylocomium 
splendens cover (Var. X91) was a significant independent variable in four 
equations: black spruce density, age of maximum PAI and mean annual basal 
area increment; and aspen mean annual basal area increment (Tables OAL O)i. 
The only significant correlation was positive between Hylocomium splendens 
and black spruce mean annual basal area increment (Table 9), indicating 
greater black spruce basal area growth on well drained sites where Hylocomium 
is likely to be a dominant moss as opposed to poorly drained sites where 
Sphagnum spp. are likely to be dominant. 

Pleurozium schreberi is also a common feathermoss of the boreal forest. 
Its cover was greatest in the lodgepole pine forests (31.9%: Table 3), 
followed by the black spruce forests (18.6%: Table 5), white spruce forests 
(10.6%: Table 4) and aspen forests (3.0%: Table 6). Pleurozium schreberi 
cover (Var. X92) or its log (Var. X93) depending on which gave the better 
linear relationship with the dependent variable in question, was a signifi- 
cant independent variable in three equations: aspen density, age of maximum 
PAI and mean annual basal area increment (Table 10). All the simple 
correlations between Pleurozium schreberi cover and the dependent variables 
were non-significant. 

Polytrichum commune, another boreal forest moss is locally abundant 
especially on mineral soil after fire, logging or other surface disturbances. 
This moss species is not abundant in the sampled forest populations having 


its greatest cover in the white spruce forests (0.6%: Table 4) and less 


i ale Pen ae 
¥ i a 


sty. S55 555 eet Jenks i et 
aigrncas apes Uae wale 
ire Laisa ehh iain nila 


Wines soars apa an ‘ied 


s. 


6 abba see 0a) sauenon se 


ry i yee ‘ hae 
gical’ sod) atvena? ‘sake sivindto 


favs atts Shee ans b eee Prue 


we ak obde tse 
- aka. Ran Ra eee Pe ey teint: ohana 
LAO), Heel tae Penh Sewer inne asaen sages De 

Ge. yawwdeed pea) alae ndisaretnee’ } 

Pye, eee Ae J rveeates Wak | HHS ee atic ni 1K 


i698 Dba tasd lier Ag Asunve: aoe levee 


ee] aw 


= 
Cc 


ihacag vhs nee a! wom Jeu beak 
stdin? wh a ddd § oni 

cy innit a Orta eS» decisti % inane! at ‘contend oi 
(Laat: yee, EY eeaeaes waa deat ads nk deed 
(2 elawh \odtig BA nieusod. shunye Goatd edd) | 
ove nitteon gel ce Bae ae eteatek yYisqeb han 
Tr Stale) pet ak) 160 (eR aa 
aw poten Oy ilatier taaterages: ort rbd ke 

& saath Dewiaaianger owen nk aiaelase itt 


1 i eee) . iG “J LENT 


Sie Ws Si 


j She aN » aij phere) PY ie ee #9 %8 isesa Laws 
sav dienes’ sad ae Rae puetiondne ne leo exgvcnnsd toed | 


ca>>) 2b gate eee Gieeed sodzone .anynnos mutsheg ie 


‘ 


Va Sree 


shennds30.t spade wees galerie aetet 26774 fton-igteria ‘rm yhisieg 
tatoos jose’ Secieogh kd 4. snehoude son 8b /eRsaeee ‘ 


gitver aawigelyy 


et2eeot, eeguye et itdw eit ni 32) 


f= i ye 1. stat -2d40) 


iS? 


than 0.5% in the black spruce, pine and aspen forests (Tables pte 0) 
Polytrichum commune cover (Var. X94) was a significant independent variable 
in only one equation: white spruce density (Table 8), with which it is 
positively correlated, perhaps as a reflection of greater seedling density 
On mineral seedbeds in the early stages of forest development. 

Total moss cover or the cumulative cover of the three species dis- 
cussed above plus others,was greatest in the black spruce forests (63.3: 
Table 5) followed by the lodgepole pine forests (45.0%: Table 3), white 
spruce forests (40.7%: Table 4) and aspen forests (9.6%: Table Gi... socal 
MOSS cover (Var. X95) or its log (Var. X96), depending upon which gave the 
better linear relationship with the dependent variable,was a significant 
independent variable in five equations: lodgepole pine PAI per tree and 
mean annual basal area increment; white spruce mean annual basal area 
increment; and black spruce MAI and PAI per tree (Tables 7-9). Only the 
Simple correlations between total moss cover and pine and black spruce 
PAI per tree were significant and both were negative (Tables k Ry coat bs) Se 
perhaps indicating suppressed growth in old dense stands. 

Total lichen cover, comprised primarily of Cladonia and Cladina Spp. 
plus Peltigera aphthosa was greatest in the lodgepole Pine forests (1.63: 
Table 3), followed by the black spruce forests (0.8%: Table 5), white spruce 
forests (0.4%: Table 4) and aspen forests (0.1%: Table 6). Lichen abundance 
was greatest in the lodgepole pine - blueberry - lichen (Type 14) forests 
on well to rapidly drained loamy sand to sand textures. Lichen cover 
(Var. X97) was a significant independent variable in four equations: lodge- 
pole pine MAI, site index, and PAI per tree (Table 7) that were negatively 
correlated with lichen cover (Table 17), and black spruce density (Table 9), 
that showed no significant correlation. The negative correlation of pine 


Productivity with lichen cover is probably a reflection of poorer pine 
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growth on the dry sites where lichen cover is likely to be greatest, 
Illingworth and Arlidge (1960) found that abundant Cladonia spp. were 


indicative of poor lodgepole pine site conditions in central British 


Columbia. 
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J. SQURCES OF ERROR IN REGRESSION ANALYSIS 


Several sources of error are often inherent in forest - site studies 
(Duffy 1964a, Ralston 1964, and others). These are discussed with reference 


to the present study. 


(1) One of the factors investigated may have been truly unrelated to tree 
growth but remained correlated due to chance. In the present study, 
attempts were made to derive plausible explanations for significant 
correlations between independent and dependent variables. Definite 
causal relationships were however, difficult to determine. The 
validity of certain factors may often be checked by an analysis of 


new data. 


(2) One must be careful not to disproportionately emphasize soil or 
vegetation properties, bearing in mind that vegetation response to 
change is more rapid than that of soils (Grigal and Arneman 1971). 
In the present study, consideration was given to soil, site and 
vegetation properties and the contribution made to the regression by 


the vegetation variables was determined. 


(3) The scale of site factor classification is important. Broadening or 
narrowing of the scale may yield closer relationships with the 
dependent variable (Grigal and Arneman 1971). The effect of 


narrowing and broadening variable scales was not evaluated in the 


present study. 


(4) The method of measuring a given independent variable may not have 
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been sufficiently refined to have shown a correlation with the 
dependent variable. For example, quantitative determinations in the 
laboratory of soil texture, for example, may have yielded closer 
relationships between soil texture properties and the dependent 
variables than those attained using textural classes as determined 


imasthe field: 


Large differences in local climate may not be considered. Due to 

the very limited availability of climatic data pertinent to the pre- 
sent study, the effect of local climate could not be evaluated except 
in a very broad sense as related to differences in elevation, for 
example, Ralston (1964) suggested that meteorological differences 

can be partly resolved by the subdivision of large geographic areas 
into relatively homogeneous climatic provinces. The climatic 
divisions made by McIver et al. (1972) for Central Alberta were of 
limited utility in the present study due to the very small mapping 


scale. 


The correlation between dependent and independent variables may be 
nonlinear. This source of error was minimized in the present study 
by making scattergrams of dependent versus independent variables 

as described in the methods section, and applying the appropriate 
transformation to the independent variable to approximate a linear 
relationship. 

Stand history effects could not be ascertained. Sites with 
abnormalities in growth caused by pathogens, insects and competing 
vegetation could often be avoided but effects associated with genetic 


origin are unknown and must be regarded as sources of error. 
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The growth factors or "independent" variables may be effects rather 
than causes of tree growth (Broadfoot 1969). Independent variables 
used in the present Study that could be interpreted as effects of tree 
growth include litter cover, deadfall cover, thickness of the organic 
horizon, the percent cover of the various plant species and possibly 
soil structure and consistence. However, all the properties mentioned 


above are undoubtedly influenced by relatively permanent soil properties. 


The independent variables should be sampled uniformly throughout 

their range, and a uniform distribution should be used if linearity 
is uncertain (Demaerschalk and Kozak 1974). An attempt was made in 
the present study to sample stands over wide range of environmental 
variation. Linearity of independent variable response to dependent 


variables was handled as in (6). 
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V SUMMARY AND CONCLUSIONS 


Fifteen forest vegetation types were identified from data from 
137 plots in the Wapiti map area, Alberta (NTS 83L) using 
association table analysis and a two-dimensional ordination and 
were described in terms of floristic composition, forest 
productivity and soils within Mixedwood, Lower Foothills, Upper 
Foothills and East lope Rockies (Rowe 1972) Forest Sections. The 
fifteen forest types with mean productivity (expressed as 
Baer) are as follows: (1) Populus tremuloides/Viburnum 
edule/Rubus pubescens (4.3 ), (2) Populus balsamifera/Rosa 
acicularis/Thalictrum venulosum (5.0 ), (3) Picea glauca/Rubus 
pubescens-Maianthemum canadense (4.7 ), (4) Picea glauca/Equisetum 
arvense/Hylocomium splendens (3.6 ), (5) Picea glauca/Hylocomi um 
splendens (3.6 }, (6) Picea mariana/Ledum groenlandicum/ 
Equisetum arvense (2.1), (7) Picea mariana/Ledum groenlandicum/ 
Rubus chamaemorus (0.4), (8) Pinus contorta/Viburnum edule/ 
Rubus pubescens (4.2 ), (9) Pinus contorta/Spiraea lucida/ 

Cornus canadensis (4.0 ), (10) Pinus contorta/Alnus crispa/Cornus 
Peeedenaic (3.7 ), (11) Pinus contorta/Ledum groenlandicum/ 
Pleurozium schreberi (3.1 .), (12) Pinus contorta/Picea mariana/ 
Ledum groenlandicum-Vaccinium membranaceum (3.1 ), (13) Pinus 
contorta/Menziesia glabella/Rubus pedatus (2.7 ), (14) Pinus 
contorta/Vaccinium myrtilloides/Cladonia spp. (1.3 ), (15) 


Picea engelmannii-Abies lasiocarpa/Menziesia glabella (2.4). 


Successional relationships among the 15 vegetation types were 


generalized from the plot data. Most of the 137 forested plots 
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are of fire origin. Two Major successional trends are apparent: 
lodgepole pine, aspen and balsam peplar to Picea glauca/Rubus 
pubescens -Maianthemum canadense; and lodgepole pine and upland black 
spruce to Picea spp. -Abies lasiocarpa/Hylocomium splendens forests. 
Logging and pulp-cutting activities will prevent climax vegetation 
from becoming established over large areas and will probably 

result in successional series somewhat different from those after 


fire. 


The two-dimensional ordination graphically demonstrated the 
similarities among the 137 plots and 15 vegetation types and 
showed the distribution of several individual plant species and 


the soil units within the plots and vegetation types. 


Separate multiple linear regression equations were calculated for 
83 lodgepole pine-, 30 white spruce-, 15 black spruce- and 17 
aspen-dominated plots, using mean annual total volume increment, 
stem density, site index, periodic annual total volume increment 
per tree, age of maximum periodic annual total volume increment 
and mean annual basal area increment as dependent variables. 
Independent variables included external site factors including 
elevation, slope and aspect; internal site (soil) factors including 
texture, horizon thickness, internal drainage, structure, consistence 
and soil color; stand factors including canopy cover, stem density, 
organic matter and deadfall; and vegetation factors including 

tree regeneration, individual shrub and herb species cover and 

moss and lichen cover. When both vegetation and physical site 


factors were allowed to enter the regressions, all equations were 
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Significant at the p=.01 level with large amounts of the variabili- 
ty in the dependent variable being accounted for by the equations: 
pine (44-66%), white spruce (70-94%), black spruce (98-99%) and 
aspen (94-99%). Testing the equations with an independent sample 
reinforced their validity. Important implications are several: 
(1) Several indices of site productivity including MAI, site index, 
PAI per tree and basal area annual increment can be reliably 
estimated for all four species groups within the greater part of 
the Wapiti map area. (2) Tree stem density in natural stands can 
be reliably estimated for the four tree species utilizing only site 
and vegetation factors, indicating factors other than stand history 
may be responsible for determining stem density. (3) Age of 
maximum PAI can be reliably estimated in a manner similar to that 
of stem density for the four tree species indicating that management 
decisions such as rotation age might be made by an appraisal of the 
Site whether or not a mature forest cover was present. With 
vegetation variables deleted from the regression equations, the 
equations lost much of their predictive ability with smaller amounts 
of variability being accounted for by the equations: pine (24-49%), 


white spruce (53-58%), black spruce (88-97%) and aspen (74-94%) . 


Each of the 49 independent variables examined appeared in at least 
one of the prediction equations. Those that appeared with the 
highest frequency were those related to soil moisture and to 
climate and microclimate. Several plant species, indicative of 
particular habitat conditions also entered into many of the growth 
prediction equations, greatly increasing the predictive ability of 


several equations especially those for lodgepole pine and white 
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spruce. 
Productivity of the soil mapping units could not be reliably 

compared due to within map unit variability in mean annual total 
volume increment being greater than that between map units. This 


may in part be due to a small sample size for some map units. 


It is concluded that the Wapiti map area includes several well 
organized systems which can be defined by vegetative, edaphic 

and other environmental criteria. The forest ecosystems and forest 
productivity-environment relationships described here for the 
Wapiti map area may be expected to be similar in adjacent areas 

of western and north-western Alberta foothills where climate, 
vegetation and soils are similar, though caution is advised 

where use of growth prediction equations outside the study area 


is contemplated. 
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APPENDIX 3. Locations of the 137 sampled plots 


Latitude 

54° 7.8'N 
54°11.8'N 
54°14.3'N 
54°15.4'N 
54°15.0'N 
54°21.8'N 


54°25.9'N 


54°29.4'N 
54°98.1'N 
54°19.2'N 
54°19.8'N 
54°27.1'N 
54°27.2'N 
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54°56.0'N 
54°56.0'N 
54°54.2'N 
54°38.3'N 
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APPENDIX 3 cont. 


Plot Plot 

No. Latitude Longitude No. Latitude Longitude 
51 54°08.8'N 118° 36.0'w 76 54°28.3'N 118° 59.2'w 
52 54°02.8'N 118° 43.8'W Jay 54033.22N 9 Liew iso ok 
53 54°05.2'N 118° 28.1'wW 78 54°19.8'N 118° 55.7'W 
54 54°01.3'N 118° 36.9'w 79. Salo lin /oligmesy ee Ew 
55 54°00.7'N 118° 03.0'W 80  54°24.8'N 119° 16.6'W 
56 54°00.8'N 118° 09.4'W 81. 54°24.2IN 119" 12-5tW 
57 54°33.0'N 118° 31.5'W G2) Bato 9IN iow aiediw 
58 Baaal cin J1s> 33.3'W 83, “BaeaviOIN) 9 1LoemeO oon 
59 BA8g3.8'N 118° 27.9%" 84 54°59.3'N 119° 58.7'W 
60 BAe36.0!N 218° 31.5! 85 54°58.5'N 119° 45.0'W 
61 BAga5.7!N 118 33.3'W 86,  §S5900.0'N 119s 59.4uu 
62 Baae.8IN /1lg> 42.2% 87 54°59.8'N 119° 59.2'wW 
63 B4e22.2'N 118° 36.0!wW 88 54°50.3'N 119° 24.7'W 
64 54°26.5'N 118° 44.3'W 89 54°42.5'N 118° 53.8'W 
65 54°27.21N 118° 36.8!W 90 454-42-7'N. lle 954-8 1y 
66 54°25.3'N 118° 48.8'W 91 54°43.5'N 118° 50.7'wW 
67 Eamog.4in 11g 47.5! 92 54°43.3'N 118° 56.5'W 
68 54°31.2'N 118° 59.2'W 93. 54°39.5'N 118° 59.9'W 
69 54°27.4'N 118° 46.1'W 94 54°45.9'N 118° 55.1'W 
70 54°29.5'N 118° 45.6'W 95 54°44.7'N 118° 58.5'W 
71 54°40.1'N 118° 45.0'W 96 54°44.9'N 118° 59.0'W 
72 54°25.6°N 118° 56.2'W 97 54°40.5'N 119° 01.1'wW 
73 Sioa seh 1g 59.27 98  54°40.4'N 119° 00.8'W 
74 54°23.8'N 118° 58.4'W 99 54°50.1'N 119° 25.9'W 


15 oan se 58-610) 100 54°51.5'N 119° 28.2'W 
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APPENDIX 3 cont. 


Plot Plot 
NO. Latitude Longitude No. Latitude Longitude 
101 ea51et'N 119° 2eto'w 126 54°19 e0!No) ptlos sos0 lw 
102 Bu Ae. 2IN | 119° 27-0hW 127 “54° i5e6'Nu ito oho w 
103 B4amaye3 (Ne 11O° -30.8'W 128 254° lets Neel Ome De or 
104 BA 45.6'N. «119° 35s2%W 129. -54°10,49N | 119-840. 3" W 
105 B4-46.9'N. 119° 23 /8"w 130° °$4°31.6°s. 1197 aca w 
106 Bane] 7 Nao” O57 w 131° S4223.0'N alo tad w 
107 64°49.8'N 119° 28.3'W 132). "54024. 53Ne8 0 ea5roaw 
108 64247.6'N 119° 30.1'W 1330 54023-88619 45 78'W 
109 BA e259 Ne 119° «4572 'W 134%" 54°96. 4°N | ios ase 
110 Ba-35.0°N = 119 46.1'W 135. 954-2652") blo esa 

° ° 364 J54-2776EN 1 Oamsaes tg 
VT BA 27f0'N. 119° 52,4'W 1 ; ; 

e) oO ’ Oo ' o ' 
112 54043.3'N?) 119° 49801'w 137° 54°26-1ON 119 seo 5 


3 S4°42.2'N 119 47.7'W 
Pra) «6 454°39.6'N 119° 48.6'W 
115 -454°41.3'N 119° 46.4'W 
116 54°43.1'N 119° 41.0'W 
Mig 54°41.4'N 119° 31.4'W 
Bree -54°43.8'N °° 119° 31.28'W 
ee 54°43.8'N 119° 39.4'W 
feo 954°44.3'N 119° 38.2'W 
ei) 54 43.4'N 119° 31.5'W 
122 paeie tN: 119. 36. W 
fea) 454°46.0'N 119° 32.0'W 
124 54°44.9'N 119 35.9'W 


125 54°14.4'N 119° 50.8'W 
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Parent 


a f 
No. material 
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(oho Gest ial 


Sahel 
Br 


Cn_till 
Br. 
(ere aephalpe 

Av 

(Ed Ao RI 

Av Gv 

Cn till 
Br 

Gry Gish 

Gv 


Ow 
Gm tide 


Org 
org/till 
OwS 
Cretal 
Av 

Lac till 
Lac till 
SsBr 

Br 

Mag geal l 


Bac tril 


Soil 
unit 


MBN6/St 
to Sky4 


MBN4/TOR1 
SHP1 
Iosl 
SHP1 
JIRV4 
MBN5/TOR2 
SHP1L 
JRV4 
MLB8 
TRV4A 
SHew 
KNZ 
SKY1/P 
BKM6 
MLB8 
Iosl 
DON1L 
S1P1/P 
TORL 
TOR 
DON] 


DONL 
COP1 


Classification 


Br.GL 


O.GL 


GLSZ.GL 


O.GL 


Internal 
drainage 
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Summarized site data for the 137 sampled plots 


Elevation 
(ft above 
sea level) 
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4200 
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APPENDIX 4 cont. 


Elevation 
Plot Parent > Soil Internal _ (ft above 


No. material Classification drainage sea level) Slope % Aspect 


25 Lac till DON1 O.GL MW 2500 4 N 

26 acy Sil) LDG6 Sz.Ch MW 2400 3 N 

27 Lac till LDG8 GLE.EB I 2750 0-2 N 

28 ES ART1 O.HFP W 2750 ll NNW 
29 Org “KNZ TY .M 2 2750 0-1 level 
30 ES  ARTL E.EB R 2750 9 SSE 
31 Org KNZ T.M vP 2750 3 N 

32 Cn tilt MBN6/s BR.GL MW 5050 0 level 
33 Av/till LDG3 BR.GL MW 4400 2 N 

MBNG 
34 en till MBN6 BR.GL Mw 3750 8 E 
TORL 

35 Cn til MBN6 BR.GL MW 4100 3 N 

36 Cn till MBN6 BR.GL MW 4100 2 E 

37 ere bial REB7 ©.GL MW 4100 30 E 

38 Cr “€i11 RBB7 E.EB MW 4100 43 sw 
39 ericill RBB7 BR.GL MW 4700 11 NE 
40 Gn ei 21 MBN6 BR.GL MW 4250 0 level 
41 Cn till MBN6 BR.GL MW 4450 8 N 

42 Av/Cn till MBN6 BR.GL MW 4400 2 N 

43 Cn tilt MBN6 BR.GL W 4150 18 E 

44 Lac S1P1 O.LG P 2250 level level 
45 Av/Lac LDG7 GL.GL I 2250 1 N 
46 Dan S1P1 GL.GL I 2150 level level 
47 Av-Lac DVS2 O.LG Mw-I 2200 1 N 

48 Till EDS4 GL.GL I 2800 level level 


49 Br Sapte 0.G 
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APPENDIX 4 cont. 


Elevation 
Plot Parent Soil. Internal (ft above 
No. material unit Classification draina e- sea level Slope % Aspect 
50 (Siow we ikae MBN5 GL.GL a 4750 5 N 
5a) (Gioy venlatah MBN4 O.GL MW 4250 5 E 
52 Cnet MBN6/t BR.GL MW 4800 abs! SSW 
Br COPZOr Tor 
53 Br TORL BR.GL 3g 4550 level level 
54 aus ea val MBN6 BR.GL MW Ss0 7 NNW 
5 Br COP1 O.GL Mw-I 3900 4 NW 
56 Avs BKM8 E.EB W 4400 5 NW 
Si (soy UeslnAl Yee S1P1 HU.LG/P P 2800 ~ 8 W 
hac tl, 
58 Lac till S1P1 HU.LG aL 2950 2 W 
59 erow acaba SKY1 0.G/P P 2950 3 N 
60 Race eae S1P1 O.LG 1g 3150 3) E 
61 Cn) eile or SKY1l O.LG/P I-P 3050 5 NE 
Lacwta. 
62 Crue cect EDS4 BR.GL I 3050 level level 
63 Org KNZ T.M P 2750 Ih nds) N 
64 Crivptee el EDS4 BR.GL MW 3500 Pe NW 
65 Org KNZ T.M P 3300 5 W 
66 Cnet EDS4/t GL.GL I 3600 Phas) E 
Br TORL 
67 Org KNZ T.M P 3050 level level 
68 Grey venlalal EDS4 0.GL MW 3000 5 N 
69 (Gre) eaglel EDS4 BR.GL Mw-I 3200 10 W 
70 Cn tals EDS4 BR.GL MW 3100 10 N 
yak Cnr tiled: EDS4 GL.GL MWw-I 3050 2 N 
AD) Cn till EDS4 BR.GL MW 3500 10 N 
73 Cn e222 SKY1 HU.LG P 3900 2 NE 
74 Cn till EDS4 BR.GL MW 3900 10 s 
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Elevation 
Plot. Parent Soil. Internal, (ft above 
No. material unit Classification drainage sea level) Slope $% Aspect 
75 Cheeta it EDS4 BR.GL MW 3850 10 Ss 
Br 
76 Crit EDS4 BR.GL MW 3500 3 s 
77 Cnr cial EDS4 OvGL MW 3250 3 E 
78 Cnyce dsl SKY1/P R.G P 3850 7 E 
79 En tril EDS4 GLBR.GL MW-I 3950 1 E 
80 Cn till MBN6 BR.GL MW 4600 2 N 
81 Gny tak MBN6 BR.GL MW-I 4550 3 SE 
82 Cnepla EDS4 GL.GL iD 4100 level level 
Br TORL 

83 Br TORI GL.GL I 3500 als s 
84 CnPtli Yor EDS4- 

Lac stat DONL O.GL 1 2850 3 N 
85 Cymer I DON1 O.GL MwW-I 2550 0 level 
86 Lac till DON1 GL.GL I 2800 0 level 
87 Lac till DON1 GL.GL i 2850 0 level 
88 ES HRT3 BR.GL MW 2450 3 W 
89 Br TORL 0.GL MW 2650 10 W 
90 Br TORL O.GL MW-I 2700 VS N 
91 Br TORI 0O.GL MW 2550 8 N 
92 Br TORL GL.GL bs 3100 3 N 
93 Br TORL GL.GL I 3250 0 level 
94 entail EDS4 0.GL MW-I 2750 9 NE 
95 Av Gv BKM8- O.EB MW 2700 in) NW 

shia 2600 8 NW 

96 hacudedel S1P1 0.G I 
97 Enger il EDS4 BR.GL MW 3150 15 W 
a8 Lac till DONL GL.GL I 3250 0 level 
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APPENDIX 4 cont. 


Elevation 
Plot Parent soil. Internal (£t above 
No. material unit Classification drainage sea level Slope $ Aspect 
99 ES HRT1 E.EB R 2550 15 N 
100 ES LDG8 BR.GL W 2500 0 level 
AvLac 
101 Avs GUN1 O.LG P 2550 5 NW 
102 Avs LDG8 
Br COP1 BR.GL MW 3100 5 W 
-W 
103 Avs BKM6 BR.GL W 2850 6 NW 
104 Lace ta. DON1L GLSZ.GL I 2700 4 N 
105 Br COPL O.EB W 3400 22 NNW 
106 ES LDG6 O.GL MW 3150 4 NW 
Cita EDS4 
107 ES LDG8 E.EB MW 2850 20 NW 
Br COP] -W 
108 Br COPL O.EB W 2950 ns NW 
1093 Creer? MLB8 BR.GL W 4200 18 W 
110 (axes Geniahit MLB8 BR.GL MW 4300 7 NNE 
ita bil Gre ceded MLB8 O.GL Mw-I 3300 URS) NE 
ilk . (eke {entail MLB8 BR.GL MW 3500 15 SW 
stabs) Avs MLB8 BR.GL MW 3600 5 N 
Creci ld. 
114 AvS LDG8 GLBR.GL I 3600 6 NNE 
Corer MLB8 
11s Crecead MLB8 BR.GL MW 3600 9 NW 
116 iQ gealilal EDS4 BR.GL MW 3100 22 NW 
his Avs S1P1 0.G I- 2900 10 NNE 
Lac till P 
118 LacecLi. DON1 O.GL MW 2700 15 W 
119 Cn till EDS4 0.GL MW 3000 7 NE 
120 (eyo tentulak (she EDS4 O.GL MW 2900 5 N 
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AeA Sah BS a) 


Elevation 
Plot Parent Soil. Internal (ft above 
No. material unit Classification drainage sea level Slove % Aspect 
kak Bac eae S1P1/P 0.G P 2750 4 NW 
a2 Av LDG6 or GL.GL I 3000 6 NE 
Lac till DON1 
23 ES BKM8 E.EB WwW 2750 3 NW 
124 Av LDG5 GL.GL ig 2700 3 NNE 
Lac tiLy 
25 Gracia RBB7/t O.GL W 4200 36 NE 
Br 
126 Avs IOSl CU.R W 4200 5 NW 
+CV 
127 Ow JRV4 O.GL W 3900 0 level 
128 Cr ci bw RBB7 BR.GL W 3900 17 WNW 
129 Creede): RBB7 BR.GL W 4000 SD ENE 
130 Creer: MLBS8 O.EB ie 4000 9 SE 
Br TORL 
131 Crete RBB7 E.DYB W 3800 40 SW 
is xd gestalt MLB8/t BR.GL W 4100 23 sw 
Br COP1 
JESS Av JRV4 BR.GL W 3700 10 E 
ese sootnbit 
134 Cre ta MLB8/t BR.GL MW 4000 9 NW 
Br COP1 
135 Org KNZ T.M P 4000 re} NW 
136 Gre ishe GUNL 0.G P . 4000 plat W 
AVEhy/ eet get IESE MLB8 BR.GL MW 4200 5 Sw 
Br COP1L 
Footnotes 
2 Plots 1 - 43 sampled during 1972 


Plots 44 - 83 sampled during 1973 
Plots 84 - 137 sampled during 1974 


aterials derived from the continental glaciation, s-stony. t-thin 
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_APPENIX 4 cont. 


Lacwevan. 
Lac - 
Br - 
Av - 
Gv - 
Ow - 


ES - 
Org = 


Soil Units 


lacustro - till of glacio-lacustrine origin 

stratified fine textured lacustrine deposits of glacio-lacustrine origin 
weathered sandstone (Ss) and shale unless '‘lithic' is indicated 

alluvium of stream deposition 

gravels of stream deposition 

outwash deposits of flowing water from a nearby glacier, usually gravels 
or sand (S) 

eolian sand 

organic soils. P-peaty phase 


are described in Twardy and Lindsay (1978) 


Soil classification abbreviations follow CSSC (1976) 


Soil internal drainage scale 


jek 
oy 
MW - 
Tce 
ys 
VP - 


rapid 

well 
moderately well 
imperfect 

poor 

very poor 


Elevation and tree height in meters = feet +3.28 
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Footnotes 


a 


b 


c 


m bancyees = fe neces. x 0.070 


Volume in a = fe X 0.028 
uh 


m°ha = ee x 0.230 


trees et = trees Bee ee ial 


Tree species abbreviations: 
Pl - lodgepole pine 


Sw - white spruce 


Sb =- black spruce 
Se - Engelmann spruce 
Fa - subalpine fir 


A - Aspen 


Pb balsam poplar 
B - white birch 
L - tamarack 


W - tree willow (Salix scouleriana) 
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APPENDIX 7 Vascular plants identified in forested plots within 
“me the Wapiti map area, arranged by family. 


POLY PODIACEAE 


Dryopteris spinulosa (Muell.) Watt 7/T,28/T,33/T,34/T,41/T,62/T,63/T,83/T, 
92/T,93/T,98/T,105/T,108/T,120/T1 


Gymnocarpium dryopteris (L.) Newm. 
EQUISETACEAE 
Equisetum arvense L. 
Equisetum scirpoides Michx. 
Equisetum sylvaticum L. 
Equisetum variegatum Schleich. 7/T,10/T,88/T 
LYCOPODIACEAE 
Lycopodium annotinum L. 
Lycopodium clavatum L. 83/1, 92/1, 202/77, LL4st. 
Lycopodium complanatum L. 
PINACEAE 
Abies lasiocarpa (Hook.) Nutt. 
Larix laricina (DuRoi) K. Koch 
Picea engelmannii Parry 
Picea glauca (Moench) Voss 
Picea mariana (Mill.) BSP. 
Pinus contorta Dougl. var. latifolia Engelm. 
GRAMINEAE 
Agropyron repens (L.) Beauv. NOYT 2 OL 
Calamagrostis canadensis (Michx.) Beauv. 
Elymus innovatus Beal 
Oryzopsis asperifolia Michx. 46/T, 90/1, 118/T 
CYPERACEAE 


Carex aenea Fern. 
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APPENDIX 7 cont. 


Carex Capillaris L. 

Eriophorum vaginatum L. ssp. spissum (Fern.) Hulten 14/T 
LILIACEAE 

Disporum trachycarpum (S. Wats.) B+tH. 94/1 

Maianthemum canadense Desf. var. interius Fern. 


Smilacina racemosa (L.) Desf. var. amplexicaulis (Nutt.) S. Wats. 25/T,27/T, 
M7 191s 14 94/1, 95/ phd /T,119/T,123/T, Ler 


Streptopus amplexifolius (L.) DC. 
Veracrum escholtzii A. Gray 33/T,36/T,40/1,41/5,42/T,43/3 
ORCHIDACEAE 


Calypso bulbosa (L.) Oakes 47/T,49/T,52/T,68/T,84/T,86/T,88/T,92/T,96/T, 
S771, 98/T;1L06/T 


Habenaria hyperborea (L.) R. Br. 


Habenaria obtusata (Pursh) ‘Torr. 44/T,59/T,61/T,6/7/T, /8/7T,85/T,97/T,98/T, 
POZA LOU Tp Lk s/LplZuse 


Listera cordata (L.) R. Br. 1/T,9/T,17/T,33/T,64/T,65/T,74/T,79/T,82/T, 
me oy it, t2e7 rT, 130/T 


SALICACEAE 
Populus balsamifera L. 
Populus tremuloides Michx. 
Salix bebbiana Sarg. 
Salix Spp. 
Salix scouleriana Barratt 
BETULACEAE 
Alnus crispa (Ait.) Pursh 
Betula papyrifera Marsh 


Betula pumila L. var. glandulifera Regel 1/T,15/T,29/T,31/T,35/T,37/T, 


SO/sE 111 /T 


Corylus cornuta Marsh 45/T 
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APPENDIX 7 cont. 


SANTALACEAE 


Geocaulon lividum (Richards.) Fern. 28/14. 29/0, 45/7 60/1, 10/4,S8/2 89/1) 
94/1, 96/T, 104/T, 114/77, 121/T 


RANUNCULACEAE 
Actaea rubra (Ait.) Willd. 
Anemone parviflora Michx. 2/T 
Gtemeatisiverticellaris DC. var. columbiana (Nutt.) A. Gray 94/T 
Petpnznzumagiaucumis. Wats. 1/T;°5/T, 15/1, 18/T,31/T,50/T,55/T,61/T,118/T,135/7T 
Thalictrum Venudosumerrel. 3/2, 19/T,460/T,1l23/1 
SAXIFRAGACEAE 
Mitella nuda L. 
Ribes glandulosum Grauer 69/T 
Ribes lacustre (Pers) Poir. 
Wareineteriroliata L. 34/T, 36/T, 41/1, 43/T, 93/10, 106/T 
ROSACEAE 
Amelanchier alnifolia Nutt. 
Fragaria virginiana Duchesne var. glauca S. Wats. 
Geum irivale “L.015/2;49/2,;63/1T,65/T,L2l/l 
Prunus pensylvanica L.f. 79/T 
Prunus virginiana L. 46/T 
Rosa acicularis Lindl. 
Rubus acaulis Michx. W/L 5 71 10/70 647 tel tee 
Pipe eiehenaenorus UL. 14/27,29/3131/1,67/4,139/7 
Rubus parviflorus Nutt. 25/1,34/1,68/T,92/1,103/1,105/T,107/1,123/1 
Rubus pedatus J.E. Smith 
Rubus pubescens Raf. 
Rubus strigosus Michx. 
Sorbus scopulirna Greene 


Spiraea lucida Dougl. 


Fey a’ 
vy te & 


uit) etal oye aay! MA ‘ae 
tel 
pott \iket the ru ae 


Vs f ry Fae “es Nat det 


ae tedin 
1h mei ar 


m4 PEONts€ Es e > é if wie ONE. a, TVAS veh a 


ed oe. +e Wir NS By, la ahpsdous enh 
Nap bemvee veatneey Nek. Pe 
Tet ed, BON Aus 


Teh 42 Bae 


ih 


ah ante 


qed POEL PNG TABT thete be ated 
ae Lat\O8. rte, Exes CENM | 4 
PAPO WeGhad NEO AART ASB INHE INE tid gun 


ad te’ or 


LEGUMINOSAE 


Lathyrus ochroleucus Hook. 


Vicia americana Muhl. 10/15, 12/2 /18/1, 19/1, 85/1), 87/7, 90/T, 96/7, L00/T 06/1, 


LOU t vee), b20/7T,L27/T 
EMPETRACEAE 
Empetrum nigrum L. 
: VIOLACEAE 
Pee2aragurca, J.E. Smith 91/7,117/T,122/T 
Viola renifolia A. Gray 
ELAEAGNACEAE 
Shepherdia canadensis (L.) Nutt. 
ONAGRACEAE 
Epilobium angustifolium lL. 
ARALIACEAE 
Aralia nudicaulis lL. 
Opplopanax horridum (Sm.) Mig. QA/T SOT 1 L/an ore 


UMBELLIFERAE 


Heracleum lanatum Michx. 1/T,5/T,41/1,61/T, 78/T,126/T,127/T,130/T 


Osmorhiza depauperata Philippi 
CORNACEAE 

Cornus canadensis L. 

Cornus stolonifera Michx. 
PYROLACEAE 


Moneses uniflora (L.) A. 
P24 TDI OAT 


Pyrola asarifolia Michx. 
Pyrola secunda L. 


Pyrola virens Schweigg. 


Gray 2/T,4/T 6/T 1 7/T, 22/7 59/T,61/T,81/T,99/T114/T, 


13/7, 112/T,124/T,127/T ,132/T,133/T 
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APPENDIX 7 cont. 
ERICACEAE 


Arctostaphylas rubra (Rehder and Wills.) Fern. 135/T 


Arctostaphylos uva-ursi (L.) Spreng. 10/T,16/T,19/T,23/T,30/T,44/T,87/T, 
SI oa s7 
Gaultheria hispidula (L.) Bigel. 
Ledum groenlandicum Oeder 
Menziesia glabella A. Gray 
@x<Gyeoccus microcarpus Turez. 14/2,29/1,31/T,67/T 
Phyllodoce empetriformis (Smith) D. Don 2/T,11/T 
Rhododendron albiflorum Hook 
Vaccinium caespitosum Michx. 
Vaccinium membranaceum Dougl. 
Vaccinium myrtilloides Michx. 
Vaccinium vitis-idaea L. var. minus Lodd 
BORAGINACEAE 
Mertensia paniculata (Ait.) S. Don 
SCROPHULARIACEAE 
Castilleja miniata Dougl. 19/T,83/T,85/T,87/T,100/T,106/3,122/1 
Pedicularis bracteosa Benth. 2/T,6/T,9/T,15/T,16/T,127/T,130/T,135/T 
RUBIACEAE 
Galium boreale L. 
Galium triflorum Michx. 
CAPRIFOLIACEAE 
Linnaea borealis L. var. americana (Forbes) Rehd. 


Fonicera dioica Ls var. glaucescens (Rydb.) Butters 45/T,47/T,84/T,85/T, 


S071, 9170, 95/1, 121/T 


Lonicera involucrata (Richards.) Banks 


Symphoricarpos albus (L.) Blake 25/T, 26/7 ,44/T,45/T, 101/T,118/T,121/1 


Viburnum edule (Michx.) Raf. 
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APPENDIX 7 cont. 


COMPOSITAE 
Achillea millefolium L. 
Antennaria neglecta Greene 9/T 
Arnica cordifolia Hook. 
Arnica latifolia Bong. 3/1 807 1781/4,125/T 
Aster Ciliolatus Lindl. 
Aster conspicuus Lindl. 
Petasites palmatus (Art.) A. Gray 
Petasites sagittatus (Pursh) A. Gray 18/T 


Solidago decumbens Greene 26/1,86/T 


1 Plot number/percent cover follows species with >20% plot presence not 
Shown in table 1. T = less than 0.5% cover 


2 Plot occurrence of tree species obtainable from Appendix 5. 
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APPENDIX 8 Non-vascular plants indentified in forested plots 
within the Wapiti map area, arranged by family. 


AULACOMNIACEAE 
Aulacomnium palustre (Hedw.) Schwaegr. 31/T. 
BRACHYTHECIACEAE 
Tomenthypnum nitens (Hedw.) Loeske 15/15 
DICRANACEAE 
Dicranum spp. includes predominantly D. polysetum Sw. and D. undulatum Brid. 
DITRICHACEAE 
Ceratodon purpureus (Hedw.) Brid. 15/16,86/T 
ENTODONTACEAE 
Pleurozium schreberi (Brid) Mitt. 
HYLOCOMIACEAE 
Hylocomium splendens (Hedw.) B.S.G. 
HY PNACEAE 
Ptilium crista-castrensis (Hedw.) DeNot. 
MNIACEAE 
Mnium spp. includes predominantly M. cuspidatum Hedw. and M. spinulosum B.S.G. 
POLYTRICHACEAE 
Polytrichum commune Hedvw. 
SPHAGNACEAE 


Sphagnum spp. includes predominantly S. fuscum (Schimp.) Klinggr. and S. 


Warnstorifil Russ 
MARCHANT IACEAE 


Marchantia polymorpha L. (liverwort) 4/2,105/2 


LICHENS 


PELT IGERACEAE 


; ‘ i 
Peltigera spp. includes predominantly P. aphthosa (L.) Willd. with lesser amoun 


of P. canina (L.) Willd. 
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"Other lichen" includes predominantly species of the Cladoniaceae, 
especially Cladina mitis (Sandst.) Hale and W. Culb., C. 
rangiferina (L.) Wigg., plus other Cladonia species. Stereocaulon 
tomentosum Fr. also often contributes significantly to the generally 


very small "other lichen" cover category in Table l. 


* Plot number/percent cover follows species with <20% plot presence 


not shown in table 1. T = less than 0.5% cover. 
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Canada Land Inventory Mean Annual 
Total Volume Increment Classes 


MAI Ge acver* yr +) MAI (m> ie yr +) 
greater than lil greater than 7.8 
2 rms Is 8) 6.4 — 97.8 
71 - 90 5.0 = 6.3 
Sio= 70 | 3.60.7 4.9 
31 = 50 Pe Pa ees: 
EL 230 Ones 2ek 
less than 10 less than 0.8 
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